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Executive  Summary 

This  technical  report  (quarterly)  details  the  work  for  Office  of  Naval  Research  (ONR)  by  Tennessee  Tech.  The  goal 
of  this  project — jointly  funded  by  ONR,  NSF,  and  ARO — is  to  build  a  general  purpose  testbed  with  time  reversal 
capability  at  the  transmitter  side.  The  envisioned  application  is  for  UWB  sensors  and  tactical  communications  in 
RF  harsh  environments  where  multipath  is  rich  and  can  be  exploited  through  the  use  of  time  reversal.  The  report 
summarizes  the  results  for  each  of  two  major  tasks.  Specifically,  the  project  can  be  broken  into 


•  Task  I -Theoretical  Research 

•  Task  2-Experimental  Testbed 

In  the  part  of  theoretical  research,  the  central  result  is  to  enable  non-fading  transmission  through  the  use  of  lime  re¬ 
versal.  Next,  the  optimum  pulse  shaping  is  also  considered,  to  extend  the  range.  Finally,  the  chirp  pulse  waveform — 
practical  and  useful — is  considered.  This  chirp  UWB  system,  combined  with  time  reversal,  will  be  promising  for 
long  range  sensor  networking. 

In  the  experimental  testbed  part,  progress  in  the  second  generation  of  the  testbed — with  time  reversal — is  described. 
At  this  moment  of  writing,  the  first  generation  of  the  testbed — without  time  reversal — is  working  at  the  Pi’s  Lab, 
over  the  air  with  a  bandwidth  of  more  than  500  MHz.  The  second  generation  will  be  built  upon  the  first  generation. 
The  primary  modification  is  the  arbitrary  modulation  waveform  generation  at  the  transmitter. 
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Chapter  1 


Introduction 


Civilian  and  naval  vessels  have  long  been  potential  targets  for  criminal  and  terrorism  activities.  There  are  concerns 
that  terrorists  can  ship  various  types  of  weapons  of  mass  destruction  (WMD)  to  international  ports  using  commercial 
ships  and  their  cargoes,  One  of  the  major  steps  in  preparedness  for  such  danger  is  reliable  wireless  comm  unications 
both  between  the  boarding  party  as  well  as  effective  ship-lo-ship  and  ship-to-shore  communications  once  the  first 
sign  of  threat  is  detected.  This  need  is  made  possible,  only  recently  [1],  with  the  advent  of  revolutionary  ultra- 
wideband  (L’WB)  technology  [2,  3,  4,  5,  6,  7,  8,  9],  This  proposed  research  is  motivated  for  RF  challenged  envi¬ 
ronments  such  as  hospital,  coal-mine,  intra-ship,  intra- vehicle,  intra-engine,  manufacturing  plants,  assembly  lines, 
nuclear  plants,  body  area  network  sensors  sunounded  by  vehicles  and  tanks,  etc. 

In  the  part  of  theoretical  research,  the  central  result  of  this  report  is  to  enable  non-fading  transmission  through  the 
use  of  time  reversal.  Next,  the  optimum  pulse  shaping  is  also  considered,  to  extend  the  range.  Finally,  the  chirp 
pulse  waveform — practical  and  useful — is  considered.  This  chirp  UWB  system,  combined  with  time  reversal,  will 
be  promising  for  long  range  sensor  networking. 

In  the  experimental  testbed  part,  progress  in  the  second  generation  of  the  testbed — with  time  reversal — is  described. 
At  this  moment  of  writing,  the  first  generation  of  the  testbed — without  time  reversal — is  working  at  the  Pi’s  Lab, 
over  the  air  with  a  bandwidth  of  more  than  500  MHz.  The  second  generation  will  be  built  upon  the  first  generation. 
The  primary  modification  is  the  arbitrary  modulation  waveform  generation  at  the  transmitter.  It  is  expected  that  the 
receiver  part  of  the  second  generation  mainly  reuses  that  of  the  first  generation. 

A  key  understanding  is  achieved  so  far:  to  further  extend  the  range,  the  design  of  modulation  waveforms  must 
be  carefully  selected,  according  to  some  criteria.  A  good  criterion  is  to  maximize  the  signal  to  noise  ratio  (SNR) 
measured  at  the  point  after  the  matched  filter.  The  frequency  selectivity  caused  by  the  huge  bandwidth  (or  pulse 
waveform  distortion  in  the  absence  of  multipath)  will  lead  to  the  unique  problem  of  pulse  waveform  optimization. 

Another  key  understanding  is  that  the  effective  channel  enabled  by  time  reversal  can  be  treated  "Non-Fading” — a 
dream  for  wireless  engineers  from  half  an  century.  The  system  design  will  greatly  simplified  under  this  non-fading 
framework.  Although  this  theoretical  breakthrough  is  reached,  the  design  principle  and  system  parameters  to  exploit 
this  new  phenomenon  are  far  from  clear. 

One  lesson  learned  from  this  research  is  that  time  reversal  is  basic,  for  a  dense  multipath  environment — a  continuous- 
time  impulse  response  of  extremely  long  delay  spread. 
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The  key  UWB  system  design  consists  of  two  parts:  (1)  time  reversal  to  deal  with  dense  multipath;  (2)  the  selection 
of  UWB  modulation  pulse  waveform — software  defined  or  cognitive.  From  a  information-theoretical  viewpoint,  the 
two  parts  as  a  whole  form  so-called  “pre-coding”. 
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Chapter  2 


Non-Fading  Transmission 


We  propose  a  time  reversal  system  paradigm  for  localization,  communication,  and  networking.  UWB  transmission 
may  lead  to  (almost  deterministic)  no-fading  communication — a  dream  for  wireless  industry  over  half  a  century — to 
enable  some  applications,  for  which  narrowband  techniques  may  fail. 

The  transceiver  structure  has  been  the  central  topic  in  the  UWB  community.  The  Pi’s  research  roughly  follows 
this  trend.  First,  the  RAKE  structure  [10,  11,  12,  13,  14,  15,  16,  17]  is  used,  and  found  too  costly — after  a  long 
painful  search,  due  to  the  large  number  of  paths  (fingers)  and  timing  requirement  [6,  7,  8,  9].  In  particular,  pulse 
waveform  distortion  will  cause  a  design  problem,  first  recognized  (one  decade  ago)  by  the  PI  based  on  his  channel 
model  [10,  II,  12,  13,  14,  15,  16,  17],  which  was  accepted  into  IEEE  standard  [18],  The  orthogonal  frequency 
division  multiplexing  (OFDM)  system  [19,  20]  is  slightly  more  complicated  than  the  RAKE  system.  Secondly,  the 
transmitted  reference  (TR)  structure  [21,  22,  23,  24,  25,  26,  27,  28,  29,  30,  31,  32]  is  used  to  reduce  the  timing 
requirement,  and  avoid  the  need  of  channel  estimation.  Thirdly,  for  sensor  networks  of  low-cost  and  low  power, 
non-coherent  energy  detection  (analog  frontend)  [33,  34,  35,  36,  37,  38,  39,  40,  41,  42,  43]  is  preferred  to  the  TR, 
although  both  systems  are  adopted  in  IEEE  802.15.4a.  Inter-symbol-interfcrencc  (ISI)  limits  the  achievable  data 
rates  of  the  TR  and  energy-detection  systems. 

Fourth,  for  data  rate  higher  than  the  TR  and  energy-detector,  time  reversal  at  the  transmitter — perhaps  viewed  as 
continuous-time,  channel-matching  filter,  or  generalized  RAKE  [17]  at  the  transmitter — is  combined  with  transmit¬ 
ted  reference  structure  [44,  45,  46,  47],  and  non-coherent  energy  detection  [48,  49,  50].  It  is  even  extended  to  the 
UWB  multiple  input  and  multiple  output  (MIMO)  system  [8,  52,  53,  54,  55,  56,  57,  58].  Time  reversal  has  been 
used  in  UWB  radio  [59,  60,  61,  62,  63,  64,  65,  66,  67,  68,  69,  70,  71,  72,  73].  It  is  used  extensively  in  acoustics 
by  Fink’s  team  [74,  75, 76, 77,  78,  79,  80,  81,  82, 83,  84,  85,  86],  and  later  to  underwater  acoustic  communications 
mainly  by  Kuperman’s  team  [87,  88,  88,  89,  90,  9! ,  92,  93,  94,  95,  96,  97,  98,  99,  100,  101,  102,  103]  to  deal  with 
rich  multipath  phenomenon  which  is  similar  to  UWB.  It  has  been  first  used  in  wire-line  communication  by  Bell 
Labs  [104,  105]  and  IBM  [106].  In  spirit,  time  reversal  is  similar  to  the  classical  phase  conjugation  concept  in 
optics  [107,  108,  109],  and  retro-directive  array  in  microwaves  [110,  111,  112,  1 13,  62,  114,  115,  116,  117,  118, 
1 19,  120,  121,  122,  123,  124,  125,  126].  The  objective  of  all  the  above  schemes  is  to  exploit  the  channel  itself  as 
part  of  channel-matched  transceiver.  This  principle  has  the  advantage  of  real-time,  since  the  operation  is  finished 
in  a  continuous-time,  analog  form,  in  contrast  to  the  algorithm  convergence  of  the  discrete-time  form  used  by  such 
systems  as  equalizers  [127,  128,  129,  130,  131,  132,  133],  A  key  observation  [53,  52,  134,  135,  136]  has  been 
recognized  that,  if  lime  reversal  is  always  used  at  the  transmitter  side,  no  fading  needs  be  explicitly  addressed  in  a 
UWB  impulse  radio.  This  effect  derives  from  the  (multipath  and  frequency)  diversity  gain  offered  by  the  huge  signal 
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bandwidth.  The  proposed  research  is  based  on  this  key  effect — to  exploit  multipath  and  signal  bandwidth  (or  pulse 
duration). 

With  time  reversal — of  the  channel  impulse  response  (CIR)  illustrated  in  Fig.  4.3 — at  the  transmitter  side,  the 
effective  channel  between  a  pair  of  a  modulation  waveform  and  a  demodulation  waveform  is  the  autocorrelation  of 
the  CIR — illustrated  in  Fig.  2.1(c).  This  new  channel  is  much  more  stationary  than  its  old  counterpart  [134,  135], 
An  analogy  is  as  follows.  If  the  CIR  is  a  realization  of  a  random  process — the  first  order  statistics  is  changing 
rapidly  from  one  realization  to  another  realization,  the  second  order  statistics  is  much  more  reliable  and  stable.  Its 
time-varying  property  highly  relics  on  the  number  of  multipath  which,  in  return,  is  determined  by  the  pulse  signal 
bandwidth,  as  a  function  of  RF  environments.  The  significance  of  this  proposal  is  to  convert  the  conventional  CIR 
problem  into  a  new,  autocorrelation  problem,  and  systematically  establish  a  theoretical  framework  based  on  this 
principle.  This  understanding  is  in  parallel  with  the  historical  development  of  the  white  Gaussian  noise  theory  in 
early  1950s  [137,  138,  139,  140,  141,  142,  143,  144,  145,  146,  147,  148,  149,  150,  151,  152]:  what  really  matters  is 
not  the  CIR's  first  order  statistics,  but  its  second  order  statistics! 

Once  the  key  conceptual  breakthrough  is  made,  the  formulation  of  this  idea  in  a  mathematical  framework  is  rela¬ 
tively  straightforward:  The  system  model  consists  of  a  linear-time-invariant  (LTI),  dispersive  channel  with  additive 
(Gaussian  or  non -Gaussian)  noise  and  interference  [20].  We  have  converted  a  rime-varying ,  dispersive  (wireless) 
channel  into  a  time-invariant,  dispersive  channel — that  is  much  easier  to  handle.  An  exhaustive  research  has  been 
done  regarding  this  classical  system  model  [153,  154,  19,  155,  156,  133,  157,  158,  159,  160],  since  it  is  the  model 
used  for  wire-line  [  153, 1 56,  133, 1 57, 19] — DSL  and  VDSL  [19] — for  which  the  Bell  systems  have  the  deep-packet 
resources  [127,  128,  129,  130,  131,  132,  133].  The  well-understood  mode!  is  an  excellent  tool  for  the  proposed 
research,  and  paves  the  way  for  our  next  engineering  breakthrough.  As  one  example,  the  modal  modulation  used 
by  Holsinger  (1964)  [153]  is  a  starting  point  for  our  new  adventure.  Time  reversal  is  used  at  the  transmitter.  The 
so-called  lime  reversal  modal  modulation  is  proposed  here.  Our  system  model  consists  of  an  effective  LTI,  dispersive 
channel1  with  additive  (Gaussian  or  non-Gaussian)  noise  and  interference.  The  simplified  mathematical  framework 
will  lead  to  a  better  design,  since  the  dimensions  of  bandwidth  and  multipath  can  be  focused  on,  through  the  effective 
LTI  channel.  The  unified  theory  links  up  all  the  fundamental  system  parameters  such  as  bandwidth,  multipath,  and 
space,  and  box  them  into  an  effective  CIR.  In  practice,  all  we  only  need  to  do  is  to  replace  the  LTI  model  of  the  clas¬ 
sical  theory  with  our  effective  LTI  model:  it  is  that  simple!  Now,  the  open  question  is  to  spell  out  all  the  necessary 
conditions — such  as  pulse  duration  and  waveform  as  a  function  of  RF  environments — for  this  new  framework  to  be 
valid. 

Built  upon  the  above  new  understanding,  theory  and  prototype  testbed  are  proposed  to  test  and  validate  this  new 
system  design  paradigm.  A  5  GHz  UWB  radio  testbed  has  a  special  meaning  in  the  context  of  the  new  paradigm, 
since,  for  indoor  channels,  2  GHz  has  been  found  to  be  the  critical  bandwidth,  beyond  which  the  effective  channel 
can  be  considered  as  “deterministic  (non-fading)  in  practice,  provided  receivers  exploiting  the  full  channel  energy 
are  employed”  quoted  from  [161 J2.  Mixed  signal  processing  is  the  bottleneck.  To  overcome  this  problem,  a  new 
patented  design  [162,  163]  from  radar  community  can  be  used. 

Similar  to  short  pulse  radar  [164,  165,  166,  167,  168,  169],  waveform  optimization  for  detection  is  a  basic  problem. 
For  UWB  MIMO,  we  have  derived  a  formulation  based  on  the  eigenfunctions  of  the  integral  equation.  This  problem 
may  be  related  to  the  multiuser  MIMO  problem  that  can  be  effectively  solved  via  a  convex  optimization  [170, 
171,  172,  173,  174,  175],  Bear  in  mind  the  fact  that,  since  a  practical  UWB  system  works  in  a  very  low  SNR 
range,  optimization  for  detection  (maximizing  SNR)  is  the  primary  problem — not  the  one  to  maximize  the  mutual 

'if  the  CIR  of  the  LTI,  dispersive  channel  is  h(r),  then  the  CIR  of  the  effective  LTI  channel  is  h(— r)  *  A(t )  where  is  the  linear 
convolution. 

2  In  the  long  process  of  studying  time  reversal,  the  Pi's  team  recognized  this  result,  independently  of  [161]. 
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information. 


2.1  Significance 


Time-varying  fading  is  the  hallmark  of  narrowband  wireless  communications.  It  seems  safe  to  say  that  almost  all  the 
major  progress  made  in  wireless  industry,  such  as  MIMO  and  OFDM,  has  been  motivated  to  exploit  fading  through 
(multipath  and  spatial)  diversity.  Fading  caused  by  multipath,  however,  may  be  negligible  for  UWB  communication 
using  a  sufficiently  large  signal  bandwidth,  e.g.  2  GHz  for  indoor  [161].  The  proposed  research  is  motivated  to 
explore  the  dimensions  of  bandwidth  and  multipath  in  real-life  environments — an  issue  unique  to  UWB  impulsive 
radio.  The  number  of  paths  is  roughly  proportional  to  the  bandwidth,  when  it  is  below  2  GHz  for  indoors.  In 
particular,  an  extremely  large  number  of  (time-invariant)  resolvable  multiple  paths — treated  as  continuous-time, 
dense  waveform  over  a  duration  of  800  ns,  see  Fig.  4.3 — will  be  encountered  in  real-life  applications  [1,  176,  177, 
178,  1 79],  and  will  be  investigated  as  an  example.  When  fading  can  be  ignored  in  practice,  the  mathematical  tools  are 
simplified  to  a  level  accessible  by  engineers  with  only  undergraduate  education  [20,  180],  The  body  of  knowledge 
based  on  fading  communication  mathematics  [181,  182] — generally  very  difficult  for  engineers  and  students,  even 
at  a  graduate  level — can  be  circumvented,  perhaps  for  the  first  time,  in  the  wireless  communication  curriculum.  The 
new  mathematical  framework  will  pave  the  way  for  a  novel  system  paradigm  and  a  new  understanding  of  the  mystic 
multipath. 


2.2  UWB  Short  Pulse  Propagation 

The  non-fading  (deterministic)  treatment  of  a  wireless  channel  is  based  on  the  experimental  study  of  the  UWB 
propagation  channel.  There  is  a  class  of  extremely  rich  multipath  environments.  This  class  of  channel  has  been 
recently  investigated  by  the  Pi’s  team  [53,  52,  134,  56,  57,  58,  135].  It  is  found  that  the  number  of  paths  is  much 
larger  than  that  of  the  office  environment — compare  Fig.  4.3  and  Fig.  2.1(d).  Communications  and  networking  in 
confined  metal  environments  [134,  135,  57]  in  this  class — Fig.2.1(a) — is  motivated  for  applications  such  as  intra¬ 
ship  [1],  intra-vehicle  [176],  intra-engine  [52],  manufacturing  plants,  assembly  lines,  nuclear  plants,  body  area 
network  sensors  surrounded  by  vehicles  and  tanks,  etc. 


2.2.1  Related  Work 


One  difficult  problem  is  communications  over  such  RF  challenged  environments  as  a  confined  metal  box.  Narrow 
band  wireless  technologies  have  proved  ineffective  in  these  environments,  due  to  resonance  caused  by  the  metal 
walls  [177,  178,  176,  179],  Chirp-like  signals  are  observed  in  a  jet  engine  [178].  Only  recently  [1],  this  problem  has 
the  promise  to  be  solved  with  the  advent  of  UWB  technology.  Lawrence  Livermore  National  Laboratory  (LLNL) 
just  demonstrated  UWB  transmission  in  a  ship  in  2006  [1  ].  The  Pi's  team  demonstrated  this  in  an  engine  [52],  and 
an  metal  cavity  [134,  135,  57].  Strict  physics  foT  the  latter  case  is  given  by  L.  B.  Fclsen  [177]  who  introduced  the 
UWB  problem  to  the  PI  in  1993,  and  W.  Geyi  [179]  who  was  introduced  to  this  problem  by  the  PI  in  2007  IEEE 
APS  conference.  In  general,  this  problem  is  open,  especially  in  the  context  of  non- fading  transmission.  The  mobility 
of  antennas  needs  to  be  considered. 
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Figure  2.1:  UWB  short  pulse  propagates  in  RF  harsh  environments,  (a)  Rectangular  metal  cavity:  (b)  CIR  in  a  rectangular 
metal  cavity.  The  pulse  has  its  spectrum  from  3  GHz  to  10  GHz.  (c)  autocorrelation  of  CIR  in  (b);  (d)  CIR  in  an  office 
environment.  The  experimental  setup  is  identical  to  that  of  (b). 


2.2.2  Research  Plan 

The  central  task  of  the  proposed  research  in  this  context  is  to  validate  the  system  design  assumption,  and  parameterize 
the  effective  channel  h(-r)  *  h(r).  When  the  testbed  is  available,  the  field  test  will  be  carried  out,  to  demonstrate 
the  system  principle  in  different  RF  harsh  environments.  A  lot  of  new  experiments  are  needed,  in  the  process  of  the 
proposed  research.  Fortunately,  the  PI,  a  veteran  in  channel  modeling,  has  the  expertise  and  the  required  equipment 
to  perform  these  tasks. 


Preliminary  Work 

The  PI  has  accumulated  more  than  one  decade’s  experience  in  this  channel  modeling  area.  His  Lab  has  two  thc- 
statc-of-the-art  channel  sounding  systems  for  UWB  communication.  The  time  domain  system  Fig.  2.3(b)  is  based 
on  digital  sampling  oscilloscope  (DSO),  Tcxtronix  TDS  7000E3.  The  time  domain  sounder  has  the  capability  to 
handle  the  MIMO  channel  sounding — not  using  virtual  array  approach.  The  frequency  domain  system  uses  a  vector 
network  analyzer  (VNA),  Agilent  N5230A  (300kHz- 1 3. 5GHz).  These  sounders  helped  to  produce  two  PhD  degrees, 
five  MS  degrees,  and  five  undergraduates  used  them  through  NSF’s  REU  (Research  for  Undergraduates)  program. 

Fig.  2.1  illustrates  a  new  phenomenon — an  extremely  long  train  of  multipath  echoes.  Fig.  4.3  conlains  much  more 
paths  than  that  of  Fig.  2.1(d).  The  measurement  is  performed  using  frequency  domain  technique  to  analyze  the 
characteristics  of  the  UWB  channel  in  the  confined  metal  environment.  The  rectangular  confined  metal  cavity  in 
Fig.  2.1(a)  measures  16  feet  by  8  feet  by  8  feet.  The  materials  of  the  walls  are  aluminum.  The  sounding  is  carried 
out  by  sweeping  a  set  of  narrowband  sinusoid  signals  (tones)  through  a  wide  frequency  band.  The  data  of  channel 
frequency  response  is  collected  using  the  VNA  based  channel  sounder  in  the  Pi’s  lab.  The  VNA  sweeps  from  3  GHz 
to  10  GHz  using  7001  points  with  frequency  step  of  1  MHz.  The  power  of  each  tone  is  10 dBm. 

2.3  Time  Reversal  Modal  Modulation 

2.3.1  Related  Work 

When  time  reversal  is  used,  the  effective  channel  is  regarded  as  deterministic  (no-fading)  channel.  The  proposed 
time  reversal  modal  modulation  is  different  from  OFDM  in  several  ways.  First,  OFDM  is  a  special  limit  of  this 
modulation,  when  the  pulse  duration  becomes  infinite.  Most  UWB  pulses  arc  short  in  duration,  similarly  to  UWB 
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II 


(a)  ft>) 


Figure  2.2:  Time  reversal  modal  modulation  is  proposed  to  simplify  the  transceiver  complexity.  The  channel  matched  filter 
h('l'h  - 1)  is  moved  to  the  transmitter  side.  The  self-reproducing  4>„{t)  can  he  defined. 

radar.  Secondly,  OFDM  can  be  viewed  as  multiple  tones  (multiple  carriers),  implying  each  tone  assumes  a  single 
frequency  value.  This  is  not  true  for  the  proposed  modulation.  The  effective  channel  can  be  diagonalized  by  its 
eigenfunctions  (eigen- waveforms) — changing  with  RF  environments,  similar  to  the  radar  problem  [165,  166,  167, 
168,  169].  These  eigen-waveforms  arc  used  for  modulation,  and  may  be  wideband  (say  100  MHz).  Thirdly,  the 
pros  and  cons  of  this  modulation  relative  to  OFDM  are  unclear,  in  the  context  of  system  implementation. 

23.2  Research  Plan 

The  proposed  system  formulation  for  optimum  modulation  waveform  is  based  on  two  assumptions: 

•  The  channel  is  linear,  lime-invariant  (LTI)  for  every  symbol  transmission. 

•  The  spectra  of  noise  and  interference  signals  are  available  at  both  the  transmitter  and  receiver. 

The  first  assumption  has  been  established  recently  by  the  Pi's  group,  through  experiments  [53,  52,  134,  56,  57,  58, 
135],  in  the  framework  of  time  reversal.  Knopp’s  team  reaches  the  same  result  [  161  ],  from  a  channel  modeling  view. 
The  second  assumption  is  the  practice  in  current  industry. 

The  mathematical  formulation  for  the  so-called  modal  modulation  traces  back  to  the  well  known  work  of  Holsingcr 
(1964)  at  MIT  [153,  19],  originally  motivated  for  telephone  line — widely  regarded  as  the  first  work  for  OFDM — 
regarding  optimum  modulation  for  a  deterministic,  dispersive,  LTI  channel  with  memory.  It  is  found  that  OFDM  is 
an  asymptotic  case  of  Holsinger’s  modal  modulation,  when  the  duration  of  the  modulation  pulse  goes  infinite.  For 
a  narrowband  technique,  the  modulation  pulse  is  a  sinusoidal  pulse  of  infinite  duration.  Since  modal  modulation 
is  optimum — to  be  made  more  cleaT  later,  OFDM  is  accepted  widely  without  question.  With  the  advent  of  UWB 
devices,  we  must  reexamine  the  practice  for  narrowband  techniques.  With  the  aid  of  the  effective  channel  introduced 
in  Section  2,  the  mathematical  framework  of  Holsinger’s  is  sufficient:  Our  task  is  to  replace  his  LTI  channel  model 
with  the  effective  channel  model — comparing  Figs.  4.6  and  4.7. 

It  is  found  recently  by  the  Pi’s  team  that  modal  modulation  is  not  equivalent  to  OFDM  for  UWB  modulation  pulses 
of  short  finite  duration.  It  appears  that  this  finding  requires  a  paradigm  shift  in  that  UWB  transmission  does  nol  have 
to  use  a  sinusoidal  pulse! 

Single  Input  and  Single  Output  (SISO)  Let  us  fix  some  notations  for  more  precise  discussions  about  Holsinger’s 
modal  modulation,  illustrated  in  Fig.  4.6,  for  colored  noise  and  an  arbitrary,  finite  observation  interval.  If  x(t)  is 
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the  input  to  the  LT1  channel — based  on  assumption  1 — with  channel  pulse  response  h(t).  The  symbol  duration  of 
ac(t)  is  [0,'/s],  and  the  arbitrary,  finite  waveform  duration  of  x{t)  is  [0, 7'j — which  implies  that  x(t.)  is  equal  to  zero 
if  t  is  between  T  and  The  output  waveform  r  (?)  =  fnr  x  (r)h  (t  -  t)  dr,  0  <  t  <  T\,  for  an  arbitrary,  finite 
observation  interval  [0,7i].  A  colored  noise  with  a  spectral  density  N(f) — based  cm  assumption  2 — defines  a  kernel 
function  K\[t,  s)  by 


r 


&n(t  —  a)K\{(jys)dcT  —  —  s),  0  <  i, s  <  T\ 


(2.1) 


(2.2) 


(23) 


where  Hn('r)  =  defines  a  function  K(t ,  s)  by 

fTl  fTx  ,  s 

K(t,s)=  I  I  h(a  —  t)K\(aJs)h(p  -  t)dadp  0  <t,s<T 
J  0  Jo 

and  define  a  possibly  infinite-size  set  of  eigenfunctions*  #*(0*  through  a  integral  equation  [183,  184] 

An  ■  I pn(t)  =  f  K(t,T)<pn(T)dT,  0  <t<T 

Jo 

That  is  a  function  </?„(0  that  when  convolved  with  over  the  interval  [0tT]  reproduces  itself,  sealed  by 

a  constant  An,  called  eigenvalue.  The  ifn(l)  n  —  l,2,...,oo  form  a  set  of  complete*  orlhonormal  basis  func¬ 
tions,  so  x(t) — an  arbitrary  modulation  waveform  of  finite  duration  T— can  be  expanded  in  terms  of  x(f)  - 

oo 

^2  xntpn  (0i  0  <  t  <  T  and  xn  —  JQ  x  (£)  <pn  (t)dt,  n  =  1, 2, ...  ,00,  so  the  energy  of  the  transmitted  sig- 

n=0 

00 

nal  waveform,  Et  -  f*  )x  (£)|2di,  can  be  expressed  as  ft*  =  When  white  Gaussian  noise  is  assumed,  a 

00 

simple  matched  filter  can  be  used:  The  output  signal  to  noise  ratio  (SNR)  is  2 E/JV>,  where  E  =  ^  and 

n- 1 

No/ 2  is  the  double  sided  noise  spectral  density.  So  the  probability  of  error  for  the  detector  is  given  by 


1  i 

Pe  ~  "AT  /  / _ exp(--i2)ttc 

v27T  J^/2E/N„  1 


(2.4) 


In  order  to  minimize  lJe,  SNR  should  be  maximized.  Since,  by  convention,  Aj  >  A2  >  Ay  >  •  •  ■ ,  if  the  energy  of 
x(?)  is  fixed,  maximizing  E  means  a;(t)  =  xnp\ (?) — the  optimal  modulation  waveform.  So,  the  maximum  SNR 

can  be  expressed  as  SNRma*  —  with  the  maximum  eigenvalue  Xmaz  —  A|  for  our  above  notation.  This  is 

ideal  for  range  extended  UWB  sensors — the  primary  goal  of  this  project. 


From  another  (capacity)  point  of  view,  the  modal  modulation  system  have  many  parallel  orthogonal  channels,  and 
each  channel  corresponds  to  one  eigenfunction  We  can  use  the  classical  water  filling  scheme  to  calculation 

coefficients  xn  n  =  1 , 2, . . . ,  00,  based  on  the  energy  constraint  of  transmitted  signal  x(t). 

A  new  time  reversal  modal  modulation,  illustrated  in  Fig.  4.7,  is  proposed  to  simplify  the  transceiver  complexity. 
The  h(Tft  -?),  time  reversal  filter,  is  moved  to  the  transmitter  side.  Without  loss  of  generality,  we  use  white  Gaussian 
noise  as  an  example  to  illustrate  the  concept.  Similarly  to  the  modal  modulation,  the  self-reproducing  eigenfunctions 
4>n(t)  with  eigenvalues  illustrated  in  Fig.  4.7  can  be  defined  as 

fin  ■  <i>n{t )  =  f  Ehh(t  -  T)<j>„(T)dT,  0  <  t  <  T 

Jo 

where  Rhh{t)  —  MO  *  M— 0*  are  four  proposed  tasks  in  the  S1SO  system: 


(2.5) 
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Figure  2.3:  UWB  MIMO  System,  (a)  MIMO  System;  (b)  Channel  Sounding. 

•  Concept  proof  The  no-fading  mathematic  framework  needs  to  be  proved,  and  made  more  precise,  by  spelling 
out  system  parameters  such  as  bandwidth  and  waveform.  Real-world  experimental  data  will  be  used  for  this 
purpose. 

•  Optimization  of  modulation  waveforms  over  different  RF  environments  This  requires  modulation  wave¬ 
form  optimization,  as  a  function  of  waveform  duration  T,  defined  in  Eq.  (4.5). 

•  Interference  spectrum  estimation  and  its  impact  Cognitive  radio  must  have  the  function  of  “detect  and 
avoid”  other  narrowband  radios  such  as  Wi-Fi  and  WiMax.  This  is  a  challenging  problem,  even  for  a  narrow- 
band  radio.  One  simple  method  is  to  estimate  the  background  interference  over  the  ultra-wide  bandwidth — say 
5  GHz — and  optimize  the  waveform,  based  on  the  key  equation  (4.5). 

Multiple  Input  and  Multiple  Output  (MIMO)  Let  us  consider  the  system  block  diagram  of  the  time-reversed 
impulse  MIMO  in  Fig.  2.3.  Let  us  denote  hnm{t)  the  channel  impulse  response  (CIR)  relating  the  m-th  element  at 
the  transmitter  to  the  n-th  element  at  the  receiver.  The  set  of  impulse  responses  is  called  the  propagation  operator. 
If  pulse  waveform  vector  a(t)  =  [o]  {*),  02(f),...,  of  finite  duration  T  is  sent  from  the  transmitter,  then  the 

vector  b(£)  —  [fci(O) &2(f)t —,Ajv{f)]t  is  received  at  the  receiver,  for  an  arbitrary,  finite  observation  interval  [0,'4'j. 
The  superscript  “f”  represents  the  conjugate  transpose  (Hermitian)  of  a  vector  or  matrix.  It  follows  that 

b(0  =  H(0  *  ■((),  0  <  t  <  T2  (2.6) 

where  H(f)  is  the  matrix  of  jV  x  M  with  elements  of  hnm{t).  The  notation  of”*”  represents  elemcnt-by-elcment 
convolution.  On  the  other  hand,  due  to  the  spatial  reciprocity,  if  one  sends  a  signal  %,n  =  1 , ...,  N  from  the  receive 

array,  the  signal  m  =  1 . M  is  obtained  in  the  transmit  array.  Defining  column  vectors  e  and  f  as  a  and  b, 

respectively,  it  follows  that 

f(t)  =  Hf(i)  *e(t),  0  <t<T2  (2.7) 

As  a  consequence,  H(f)  permits  one  to  calculate  the  forward  propagation  of  impulses  from  the  transmit  array  to  the 
receive  array,  while  tf(I)  the  backward  propagation  from  the  receive  array  to  the  transmit  array.  If  the  time  reversed 
matrix  H*(—  t)  is  used  as  part  of  the  precoding  matrix — in  analogy  with  Fig.  4.7,  the  equivalent  matrix  channel  is 


RfiH(f)=H(f)*Ht(-«)  (2.8) 

Eq.  (2.8)  will  contain  the  temporal-spatial  focusing — unique  to  UWB  pulse  signals.  One  objective  of  UWB  MIMO 
is  to  enable  non-fading  transmission.  Cooperative  distributed  antennas  can  be  used.  MIMO  multiuser  networking  is 
our  future  research  focus,  in  the  framework  of  time  reversal. 


in  analogy  with  (4.14),  the  self-reproducing  eigenfunction  vector  3,(f)  can  be  defined  as 


Rhh(*  -  T)On{r)dT,  0  <  t  <  T 


(2.9) 


Once  Eq.  (2.9)  is  solved,  the  rest  of  signal  processing  is  similar  to  the  SISO  system  above.  An  alternative  approach  to 
solve  for  optimum  MIMO  waveform  may  use  the  recent  result  of  information  theory,  based  on  convex  optimization 
[170, 171,  172,  173,  174,  175],  There  are  three  proposed  tasks  in  the  MIMO  system: 
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(c) 


W) 


Figure  2.4:  Comparison  of  modulation  waveforms  for  optimum  modal  modulation  for  different  environments,  (a)  wedge 
(street  comer);  (b)  rectangular  plate  (building);  (c)  waveform  comparisons:  A  passhand  channel  has  an  ideal  passhand  filter  as 
its  frequency  response;  (d)  optimum  modal  modulation  waveform  calculated  using  CIR  of  Fig.  4.3 


•  Concept  proof  Similar  to  the  SISO  system,  experimental  and  theoretical  work  is  required  to  extend  the  frame¬ 
work  to  the  MIMO  system.  The  focus  is  on  the  validity  of  the  framework,  by  considering  multiple  antennas. 
The  Pi's  lab  has  two  channel  sounding  systems  (in  both  time — in  Fig.  2.3(b) — and  frequency  domain)  for  this 
purpose. 

•  Spatio-temporal  optimization  Optimization  of  waveform  is  required,  by  taking  into  account  the  space  di¬ 
mension,  through  Eq.  (2.9). 

•  Spectrum  estimation  Sensing  the  spectrum  is  a  critical  part  of  cognitive  radio.  Algorithms  for  spectrum 
estimation  need  be  investigated. 


2.3.3  Preliminary  Work 

Compared  with  the  benchmark  of  the  time  reversal  chirp  system,  there  is  a  big  room  for  optimization  (Table  ??) — 
as  a  function  of  waveform  duration — by  using  time  reversal  modal  modulation  system.  This  has  been  justified 
from  two  different  metrics.  Different  propagation  environments  define  different  optimum  modulation  waveforms,  as 
illustrated  in  Fig.  2.4. 

From  Detection  Point  of  Hew  The  SNR  after  matched  filter  in  the  receiver  side  is  the  metric  to  determine  the  system 
performance. 

From  Data  Rate  Point  of  View  In  the  time  reversal  modal  modulation  system,  there  are  many  parallel  orthogonal 
channels,  and  each  channel  corresponds  to  the  orthonormat  eigenfunction  (f)  and  eigenvalue  from  Eq.  4.14. 
One  waveform  design  example  of  the  fa  (f)  is  given  in  Fig.  2.4.  We  assume  «i  >  «2  >  *3  >  ■  ■  ■ .  Meanwhile, 
the  bits  transmitted  in  these  different  channels  are  assumed  to  be  totally  independent.  If  L-level  pulse  amplitude 
modulation  (PAM)  is  considered,  there  are  L  equally  likely  amplitude  values.  In  time  reversal  system,  the  transmitter 
waveform  filter  is  GT(t),  white  the  time  reversal  filter  is  h(-t).  In  the  time  reversal  chirp  system,  3r(f)  is  chirp 
waveform.  When  BER  p,  is  set  to  10-4  and  the  same  PAM  constellation  is  considered,  e.g.,  for  L  =  8,  the  data  rate 
gap  (in  dB)  between  time  reversal  modal  modulation  system  and  time  reversal  chirp  system  is  shown  to  7.5,  1 1.1, 
13.9,  for  the  waveform  duration  ofT  (in  ns)  equal  to  I,  10,  100. 
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Chapter  3 


Optimum  Transmission  Waveform  in 
Canonical  Channels 


The  purpose  of  this  chapter  is  to  introduce  readers  the  time  domain  modeling  for  the  UWB  signals.  For  wide¬ 
band  signals  such  as  UWB,  by  principle,  time  domain  frame  work  appears  to  be  a  better  frame  work  to  study  its 
propagation.  The  time  domain  impulse  responses  for  a  series  of  scattering  environments  have  been  derived  for  the 
first  time.  Specifically,  several  electromagnetic  communication  fundamental  channels,  including  perfect  conduct 
electromagnetic  wedge,  large  rectangular  plates,  and  parallel  planes  waveguide,  are  studied.  The  receiving  signal, 
calculated  through  convolution  in  the  time  domain  impulse  response,  can  provide  waveform  information,  as  well 
as  the  amplitude  information  (Path  loss).  The  pulse  waveforms  passing  through  these  channels  are  compared  with 
their  original  waveforms.  The  comparison  show  that  for  wideband  signal,  pulse  could  be  extremely  distorted  by 
diffraction  mechanism.  The  results  converted  from  frequency  domain  by  applying  a  inverse  fast  Fourier  Transform 
(IFFT)  serves  as  sanity  checks  for  our  time  domain  derivation.  Given  the  time  domain  channel  impulse  response,  we 
also  study  the  optimum  transmit  waveform  such  as  the  signal  to  noise  ratio  (SNR)  at  the  receiver  side  is  maximized, 
under  the  constraint  of  fixed  transmitted  power. 


3.1  Introduction 


The  purpose  of  this  chapter  is  to  derive  the  channel  impulse  responses  for  a  series  of  environments  and  analyze 
the  corresponding  optimum  transmit  pulse  waveforms.  The  UWB  system  arc  generally  based  on  the  transmission 
of  many  consecutive  pulses.  Therefore,  by  principle,  time  domain  framework  appears  to  be  a  better  framework  to 
study  the  radiation,  propagation  and  reception  of  UWB  signals.  Moreover,  comparing  with  traditional  frequency 
domain  framework,  time  domain  framework  provides  more  physical  insight.  For  example,  except  predicting  path 
loss,  time  domain  convolution  with  TD-UTD  can  also  predict  the  receiving  pulse  waveform  and  show  the  possible 
pulse  distortion  caused  by  diffraction,  which  has  been  proved  to  be  a  significant  issue  for  wideband  signals  such  as 
UWB  [l]-[3]. 

The  essential  benefits  of  the  time  domain  framework  motivate  this  work.  However,  the  traditional  time  domain 
methods  such  as  Finite  Difference  Time  Domain  (FDTD)  and  time  domain  integral  equation  become  intractable 
when  the  width  of  the  incident  pulse  is  very  narrow  compared  with  the  geometric  dimensions  of  the  scattering 
objects,  which  is  the  case  when  we  study  the  propagation  of  UWB  signals.  In  this  chapter,  we  will  try  to  attack  the 
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propagation  problems  for  wideband  signals  directly  in  the  time  domain,  with  a  closed  form  solution.  We  derived 
the  time  domain  impulse  response  for  a  series  of  scattering  environments.  The  pulse  waveform  after  diffraction 
is  compared  with  the  incident  waveform.  The  results  converted  from  frequency  domain  by  applying  a  inverse  fast 
Fourier  Transform  (IFFT)  will  serve  as  sanity  checks. 

Since  waveform  is  distorted  by  the  channel,  we  need  to  design  transmit  waveform  in  accordance  with  the  different 
channels.  Another  contribution  of  the  chapter  is  that  we  analyze  the  optimum  transmit  waveform  based  on  derived 
time  domain  impulse  response. 


3.2  Pulse  Distortion  and  Its  Impact  on  UWB  System  Design 

3.2.1  Pulse  Distortion  Caused  by  Wide  Bandwidth 

Traditional  communication  theory  are  based  on  the  assumption  of  distortionless  propagation,  which  suggests  that  the 
receiving  signal  consist  of  a  series  of  the  replica  of  the  transmitting  signals.  This  assumption  is  valid  when  the  signal 
bandwidth  is  not  wide.  However,  for  wideband  signals,  especially  pulse  based  signals  like  UWB,  this  assumption 
might  not  hold  anymore.  As  a  simple  example.  Fig.  3.2  shows  the  waveform  distortions  for  pulses  with  different 
band  widths  after  they  pass  through  a  90-degree  wedge  (diffraction)  channel.  The  incident  pulse  in  Fig.  3.2(a)  is 
a  second  order  derivative  of  Gaussian  waveform  with  a  pulse  width  of  0.2ns,  corresponding  a  5.86  GHz  10-dB 
bandwidth  in  the  frequency  domain,  and  the  pulse  in  Fig.  3.2(a)  is  about  10ns  wide,  corresponding  a  128  MHz 
1 0-dB  bandwidth.  For  the  convenience  of  waveform  comparison,  the  amplitude  of  both  incident  waveform  and  the 
diffracted  waveform  have  been  normalized.  It  is  evident  that  the  shorter  pulse  (wider  bandwidth)  experienced  more 
distortion  after  they  passing  through  the  same  diffraction  channel.  When  the  signal  bandwidth  is  not  very  wide  (c.g., 
the  128  MHz  waveform  in  Fig.  3.2(a)),  the  distortion  caused  by  channel  is  negligible.  However,  when  the  signal 
width  is  increasing,  the  distortion  becomes  noticeable.  The  details  of  the  calculation  of  the  diffraction  waveform  for 
a  wedge  diffraction  will  be  given  in  the  Section  3.5.1. 
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Figure  3.2;  Pulse  distortion  for  pulses  with  different  bandwidths.  The  incident  pulse  for  (a)  is  about  10.2  ns  wide, 
corresponding  to  a  10  dB  bandwidth  of  128  MHz,  and  for  (b)  is  about  0.2  ns  wide,  corresponding  to  a  10  dB 
bandwidth  of  5.86  GHz.  Amplitudes  were  scaled  to  the  same  level  to  better  compare  the  waveforms 


3.2.2  The  Impact  of  Pulse  Distortion  on  System  Design 

Consider  a  communication  system  illustrated  in  Fig.  3.3.  Here,  p{f)  is  the  transmitted  pulse  waveform,  and  h(f ) 
denotes  the  CIR.  The  received  signal  before  q  (/)  can  be  written  as 


r(t)  =  p(t)  *  h(t)  +  n(f),  (3.1) 

where  n(t)  is  Additive  White  Gaussian  Noise  (AWGN)  with  a  two-sided  power  spectral  density  of  Hi/2.  Given 
s(f),  a  matched  filter  matched  to  s(-f)  provides  the  maximum  signal -to-noise  power  (SNR)  ratio  at  its  output.  The 
matched  filter  here  can  be  virtually  decomposed  into  two  filters,  with  the  first  filter  matched  to  the  CIR  h(t)  and 
the  second  matched  to  the  pulse  waveform  p(f).  As  shown  in  Fig.  3.3,  this  decomposition  can  be  done  by  setting 
C|(()  —  h(-t)  and  c2(()  —  p(— t).  conventionally  channel  is  modeled  as  a  tapped  delay  line. 

L 

h(i)=£aia(*-7i)  (3.2) 

t=i 


where  L  is  the  number  of  multipath  components  and  aj  and  rt  are  its  individual  amplitudes  and  delays.  The  received 
signals  are  just  a  series  of  replicas  of  the  transmitted  signals,  with  different  attenuations  and  time  delays. 

In  practice,  the  estimation  of  h(t)  is  generally  a  very  difficult  task,  especially  when  there  is  per-pulse  distortion 
for  the  wide  band  signals.  The  receiver  suffers  significant  performance  loss  (as  big  as  4  dB),  if  the  system  is 
designed  without  the  necessary  compensation  for  pulse  distortion  [10],  Along  with  system  performance  loss,  pulse 
distortion  can  cause  timing  and  synchronization  error  as  well.  It  has  been  shown  in  ( 10]  that  the  timing  error  caused 
by  distortion  due  to  cylinder  diffraction  can  be  much  larger  than  the  Cramer-Rao  lower  bound,  thus  is  another 
fundamental  issue  that  should  be  considered  in  system  design. 
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Figure  3.3:  General  model  for  dala  transmission 


A  new  channel  model  that  takes  into  account  per-path  distortion  is  then  proposed  in  [15] 

L 

h(t }  =  £  h,(t)  *  S(t  -  r,)  (3.3) 

(=i 

where  distortion  factor  ht(t)  includes  both  attenuation  and  distortion  of  the  signal  when  pass  through  the  channel. 

Note  that  waveform  information  can  only  be  provided  by  time  domain  framework,  which  justifies  our  effort  for  the 
time  domain  analysis. 


3.3  Time  Domain  Channel  Modeling 

In  this  section  we  study  several  simple  models  that  is  mathematically  tractable.  For  these  models,  the  analytical 
results  can  often  be  obtained  through  exploring  the  physical  mechanism  behind  it  These  analytical  results  provide 
us  unique  insight  into  the  signal  propagation  in  different  practical  environments.  In  this  section,  we  are  aimed  to 
derive  the  time  domain  channel  impulse  response  for  these  models.  Wc  will  firstly  start  with  its  frequency  domain 
model  and  then  transform  it  to  time  domain  via  Laplace  transform. 


3 .3.1  Wedge  Channel 

Rather  than  seeking  the  exact  solution  for  Maxwell’s  equations,  Geometric  Theory  of  Diffraction  (GTD)  uses  gen¬ 
eralization  of  Fermat’s  principle  to  find  high-frequency  harmonic  Maxwell  equation  solutions  for  many  complicated 
objects.  The  Uniform  Theory  of  Diffraction  (UTD)  was  proposed  to  complement  GTD  in  the  situations  where  GTD 
was  invalid  [4],  The  original  GTD/UTD  was  proposed  in  frequency  domain  and  then  was  widely  used  in  narrow- 
band  radio  propagation,  e.g.,  for  path  loss  prediction.  In  1990,  time  domain  UTD  was  firstly  derived  by  Veruttipong 
and  Kouyoumjian  in  [5]  by  applying  an  inverse  Laplace  transform  to  the  con-esponding  frequency  domain  results. 
TD-UTD  diffraction  coefficient  for  a  straight  Perfect  Electric  Conduct  (PEC)  wedge  can  be  expressed  as  [5] 


D*'h(t) 


-1 

2n\/27rsm/Jo 


4 

E  K*F{Xmtt) 

TTl—l 


(3.4) 
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where 


Kl'h  =  ±coS(i±f±^} 


(3.5) 


K Z'h  =  ±cot(*^£l) 


and 


(3.6) 


The  value  of  Xm  in  (3.6)  is  given  by 


A, 

*2 

*3 

*4 


2Lcos2(2n^-+j^--^) 


2Lras2(2n^Vj±^J) 


(3.7) 


2  Lcos2 ( 2wr  -  ) 


where  is  the  nearest  integer  solution  of  equation  2irnN±  —  0±—  ±7r,  and  L  is  the  distance  parameter  determined 
by  the  incident  wavefront.  For  plane  wave  incidence,  we  have  L  =  r  sir?  fa,  where  r  is  the  distance  along  the 
diffraction  ray  from  the  edge  to  the  observation  point.  The  rest  of  parameters  used  in  (3.4)-(3.7)  are  described  as 
follows: 

us,  h ”  soft  and  hard  boundary  conditions 


n  the  parameter  describing  wedge  angle;  wedge  angle  do  =  (2  —  n)n 

fa  the  angle  between  the  edge  and  the  incident  ray 

4>  the  angle  between  the  wedge  surface  and  the  diffraction  ray 

<(>'  the  angle  between  the  wedge  surface  and  the  incident  ray 


These  parameters  are  also  illustrated  in  Fig.  3.1,  which  shows  a  special  case  when  $  =  jr/2. 

if  the  observation  angle  is  not  in  the  transition  zone,  we  have  the  time  domain  diffraction  coefficient  for  GTD  (Keller) 
solution  [5] 


(3.8) 


It  should  be  noted  that  Veruttipong’s  TD  wedge  diffraction  coefficient  given  in  (3.4)  is  a  special  case  of  the  diffraction 
coefficient  for  curved  wedge  in  [6],  where  an  analytical  signal  representation  for  the  transient  fields  has  been  used. 
The  details  for  analytical  TD-UTD  will  be  given  in  the  Section  3.5.1. 


3.3.2  Rectangular  Building 

The  first  propagation  model  of  interest  is  rectangular  model.  Rectangular  is  one  of  the  typical  shapes  often  encoun¬ 
tered  in  a  practical  propagation  environment.  Moreover,  a  lot  of  obstacles  in  the  propagation  environment  can  be 
approximated  by  a  rectangular,  with  a  limited  degree  of  accuracy  [7].  This  model  has  been  widely  studied  in  the  past 


32 


CHAPTER  3 .  OPTIMUM  TRANSMISSION  WAVEFORM  IN  CANONICAL  CHANNELS 


z 


Figure  3,4:  UWB  pulse  diffracted  by  a  rectangular  building,  Tx:  Transmitter  antenna,  Rx:  Receiver  antenna.  4fdr 
denotes  the  distance  between  Tx/Rx  and  the  building.  The  size  of  the  building  is  2  a  x  b.  The  lower  half  plane  plotted 
by  dash  line  denotes  the  image  of  the  building  relative  to  the  ground. 


decades.  However,  most  of  research  is  focused  on  frequency  domain  such  as  path  loss  prediction  under  different 
frequency  band.  In  this  chapter,  we  investigate  its  time  domain  characteristics  such  as  pulse  distortion  introduced  by 
diffraction  in  the  edge  of  the  building.  We  also  analyze  the  optimum  transmit  waveform  based  on  this  time  domain 


channel  model. 

As  illustrated  in  Fig.  3.4,  we  consider  a  simple  UWB  channel  consists  of  a  transmitter,  a  receiver,  and  a  rectangular 
building  that  blocks  the  LOS  transmission.  In  Fig.  3.4,  dr  denote  the  distance  from  the  building  to  Tx,  Rx, 
respectively.  The  size  of  the  building  is  2  a  x  6.  The  lower  half  plane  plotted  by  dash  line  denotes  the  image  of  the 
building  relative  to  the  ground  and  is  considered  as  perfectly  absorbing. 

A  good  analysis  of  the  path  loss  prediction  in  frequency  domain  can  be  found  in  [7]  and  the  references  therein. 
Generally,  two  approaches,  Fresnel-Kirchhoff  scalar  approach  and  GTD  approach  considering  a  four-ray  propagation 
system,  are  applied  to  study  the  effects  of  rectangular  building  diffraction. 


Propagation  Model  Based  on  Fresnei-Kirchhof  Integral 

A I .  I  Frequency  Domain  Model 

Fresnel-Kirchhoff  is  one  of  the  traditional  approaches  used  to  calculate  the  propagation  loss  due  to  the  blocking  of 
a  rectangular  building.  We  Define  the  channel  transfer  function  as  the  ratio  of  the  strength  of  the  received  signal  to 
the  transmitted  signal.  Mathematically  transfer  function  can  be  expressed  as  H{ju)  —  *  Then  the  transfer 

function  describes  the  rectangular  diffraction  can  be  expressed  as  [7]: 


(3.9) 


where  A:  is  the  wavenumber  and  d  is  the  propagation  distance.  In  our  case  wc  have  d  —  4  +  dT  .In  (3.9),  variable  A 
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and  B  are  ihe  normalized  half  width  and  normalized  height  of  the  building,  given  by 


A  —  ai 


12  (dt+dr) 

A  did* 


h  W+dr) 


0.10} 


where  A  is  the  wavelength.  It  should  be  noted  that  the  validity  of  (3, 9)  is  based  on  the  condition  that  the  distances  to 
and  from  the  obstacle  must  be  greater  than  its  dimensions  and  the  wavelength,  i.e*T  4,  dr  »  aT  6,  A. 

A  1.2  Time  Domain  Model 

In  this  section,  we  will  derive  the  time  domain  channel  impulse  response  based  on  (3.9).  We  firstly  work  on  one  of 
the  integrals  and  formulate  them  to  a  form  ready  for  Laplace  transform. 


exp(-j^y2)dy  =  exp  j~(\/^y)2j dy 


■  MT 


exp  (—z2)dz 

=  ^(l-j)(l-er/c(A')) 

where  the  complementary  error  function  is  defined  as 

2 

erfc(x)  —  -7=  /  exp (—t2)dt  , 

V*  Jx 

and 

A'  =  A-  =  y/srl  , 

Here  ,s  —  ju)y  the  propagation  delay  ra  —  a2(dt  +  dr)/(2dtdrc)1  and  c  is  the  speed  oflight. 
Similarly  we  have 


(3.11) 


exp {-j^y2)dy  =  “(1  -  i)(l  ~  er/c( v/srt))  . 


(3.12) 


where  Tb  =  (P(dt  +dr)/(2dtdrc). 

Substitute  (3, 1 2)  into  (3.9)  and  after  some  manipulations  we  can  reformulate  (3.9)  as 

H{s)  =  ^  [1  -  (1  -  erfci^sTa))  (1  -  erfc{yfsTb))\  e^*7” 
where  tq  —  dfc  is  the  common  time  delay. 

Recall  the  Laplace  transform  pair 

«r/c(^sa  'XT  -  r.)  . 


(3.13) 


7 T  tyjt  -  T„ 


(3.14) 


where  the  unit  step  function  u{£)  is  defined  as 


,,,  /  1,  t>0 

“<1)  =  {o,  (<0 
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By  applying  inverse  Laplace  transform  to  (3.13),  time  domain  impulse  response  can  be  obtained 

MO  =  mt-n).  [dj  -  r.)  +  -  r.))  -  (^=  *  •  <3I5> 


When  a  — +  oo,  ra  —*  oc,  thus  (3J5)  reduces  to 

MO  =  |5(t  -  r0)  *  £  yfr^u(t  -  n) , 


(3.16) 


which  is  exactly  the  time  domain  formula  for  the  well-known  expression  for  knife  edge  diffraction  [P.503,  [17]]. 
This  serves  as  our  sanity  cheek. 


Eq.  (3.15)  consists  of  three  parts.  The  physical  explanations  for  the  three  parts  arc  strailforward.  They  represent  the 
response  from  horizontal  edge,  vertical  edge,  and  the  mutual  interaction  between  these  two  edges,  respectively. 


Propagation  Model  Based  on  GTD 


Ray  concept  is  involved  when  we  attack  the  diffraction  problem  with  GTD/UTD.  The  total  received  field  strength  is 
calculated  by  adding  contributions  from  the  various  rays,  which  is  also  referred  to  as  multipath  components  in  the 
wireless  communication.  [7]  provided  a  good  GTD-based  model  in  the  frequency  domain 


//  (jM 


(3.17) 


where  the  frequency  domain  GTD  diffraction  coefficient  for  half  plane  (zero-angle  wedge)DH,A(_fw,  t//,  t/0  can  6e 
expressed  as 


bs’h(j<j,  n',  if)  - 


e'^V A 
4tt 


1 


fO'-O'  ^ 


1 


)  ■  cos(^) 


lCOs(  2  , 

The  corresponding  angle  and  the  distance  parameters  in  (3.17)  and  (3.18)  are  given  by 


(3.18) 


=  tan  *{£) 

4>'  = 

d\  —  y/  df  +  O2 
d2  = 


^  =  il/  +  n  +  tan~ *{£)  +  tan  ~ 1  ( £ ) 
<t>  =  +  7T  +  fan_1(^-j  +  tan_I(£) 

di  -  y/d*  +  a2 
^  =  y/df+¥ 


(3.19) 


The  key  term  to  transform  (3.12)  to  time  domain  is  the  GTD  diffraction  coefficient,  which  here  can  be  obtained  by 
setting  n  -  2  and  =  n/2,  from  the  GTD  solluion  (3.8)  for  wedge  diffraction  in  Section  3.3. 1 


D*'h{t,n\n) 


zl  fl 

2tt  V  2 1 


I 
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cx«(^) 


(3.20) 
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Figure  3.5:  UWB  pulse  diffracted  by  parallel  plates. 

3.3 3  Two  Parallel  Half  Plane  Model 


Diffraction  of  a  pulse  or  a  transient  wave  (acoustic  or  electromagnetic)  by  a  half  plane  is  one  of  the  early  diffraction 
problems.  Exact  time  domain  solution  of  the  CIR  for  such  a  channel  has  been  obtained.  In  this  paper,  another 
case  for  the  pulse  diffracted  by  two  parallel  half  planes  will  be  studied.  Again,  we  firstly  borrow  the  existing  exact 
frequency  domain  solution  and  then  transform  it  to  time  domain  by  Laplace  transform.  It  will  be  interesting  to 
compare  the  field  diffracted  by  two  parallel  half  planes  with  that  diffracted  by  a  single  half  plane. 


In  reality,  a  lot  of  environments  can  be  approximately  modeled  as  the  problem  of  diffraction  by  two  parallel  planes. 
Fig.  3.5  illustrates  an  example.  We  assume  the  transmitting  antenna  is  located  far  away  from  the  planes  and  thus  the 
incident  wave  could  be  viewed  as  a  plane  wave.  Moreover,  we  assume  the  observation  point  is  far  from  the  planes 
so  that  the  distance  parameter  p  — *  oo.  We  restrict  that  the  incident  angle  <k  is  in  the  range  of  0  <  <  x  and  the 

observation  angle  <p  ^  it  ±  <&,  which  means  <t>  is  not  in  the  vicinity  of  the  geometrical  optics  shadow  boundary  and 
the  reflection  boundary.  The  transfer  function  for  such  an  channel  would  be  [  1 6] 


T  cos  [Am(sin  <$>  +  sin  <£(>)]  + 

(ab  +  sfe )  J0  cos  [fcQ(sin  0  ~ sin  ^  ~ 


/4nka 

7 


'  m— U 

V  cos  <h  cos  dhi  1  /‘im)1/3 


(3.21) 


where  o  is  the  distance  between  two  planes.  Note  that  another  assumption  made  here  is  ka  >>  1,  which  can  be 
satisfied  cither  if  the  distance  a  is  large  or  the  frequency  is  high. 


Taking  a  Laplace  inversion  (by  setting  s  =  the  time  domain  impulse  response  can  be  obtained  as 


m  =  < 


oo 

Ai  E 


+ 


11  J^0  (am-t-t)*/1 
°°  i 

E  (am)3/* 


<5  (f  —  a/c  (sin  +  sin<fo)  -  p/c)  + 

8  (f  +  a/c  (sin  +  sin<fo)  -  p/c) 

S(t.  —  a/c  (sin  </>  —  sin  0o)  —  2(2m  +  1  )a/c  —  p/c) 
6  (t  +  a/c  (sin  <f>  -  sin  -  2(2m  +  l  )a/c  -  p/c] 
8(t  —  a/c  (sin  0  +  sin  <j>o)  —  2(2 m)a/c  —  p/c)  + 

<5  ((  +  a/c  (sin  <f>  +  sin0o)  -  2(2rn)a/c  -  p/c) 


•«(<)  * 
u(t)  * 


*  (3.22) 
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where 


V2  1  sin^siD^^o  I  sm^sin|^o  /  I  I  \  c 

n  ^/p/cCOS<(>  +  oob<Pq  1  y/2n  cos<p  +  cos  <k\  \  cos  <p  cos  )  ^fpa 


(3.23) 


Physical  meaning  for  (3.22)  is  as  follows.  The  first  term  in  (3.22)  consists  of  a  superposition  of  edge  waves  from  the 
two  half  planes,  with  the  condition  that  each  is  excited  by  the  incident  pulse  alone.  The  edge  wave  is  characterized 
by  a  kernel  ofu(f)/v^  and  can  be  explained  with  fraction  calculus  as  “semi-integral”. 

/(0*^  =  r(i/2)£>-,/2{/W} 

where  T(x)  is  the  Gamma  function  and  the  fractional  integral  of  order  7  is  defined  by 

s-7 {/cm  =  f^)/lCf_e)7  */(m 

The  second  and  the  third  terms  are  infinite  sums,  which  represent  the  successive  mutual  interactions  between  the 
two  half  planes.  These  two  terms  are  characterized  by  a  kernel  of  u(t)  and  we  have 

/(*)•!.(*)=  f  nm 

j —00 

The  coefficient  or  p^n/a  decreases  rapidly  as  the  value  of  m  is  increased.  Practically,  compulation  of 

infinite  sums  in  (3.22)  requires  but  few  terms  to  obtain  reasonable  accuracy.  This  will  be  verified  in  our  numerical 
results  in  the  Section  3.5. 


It  should  be  noted  that  if  both  the  transmitter  antenna  and  the  receiver  antenna  are  located  in  the  cannon  between 
two  half  planes,  then  there  is  no  diffraction  phenomena  and  multiple  reflections  dominate  the  received  signal.  Time 
domain  C1R  for  such  a  channel  has  been  studied  in  [  1 8]. 


3.4  Time  Domain  Optimum  Transmit  Waveform  Design 

Section  3.2.2  shows  the  optimum  receiver  should  be  a  matched  filter  matched  to  the  waveform  distorted  by  the 
channel,  not  the  transmitted  waveform  itself.  However,  from  system  point  of  view,  such  a  matched  filter  alone  is 
not  enough.  We  can  further  maximize  SNR  at  the  receiver  side  by  carefully  designing  the  transmitted  waveform 
[6].  The  primary  difference  between  the  problem  in  this  paper  with  the  standard  matched  filter  problem  is  the  pulse 
distortion  caused  by  channel.  We  will  show  later  that  if  channel  is  distortionless,  then  the  SNR  at  the  receiver  side  is 
independent  of  transmitted  waveform  and  then  there  is  no  optimum  transmitted  waveform. 


3.4.1  Optimum  Transmit  Waveform  for  Distorted  Channel 

The  transmitted  waveform  optimization  problem  can  be  stated  as  follows.  Given  the  channel  impulse  response  ) 
and  fixed  transmitted  power  Pt  (average  power),  we  wish  to  achieve  the  maximum  SNR  at  the  receiver  side  by  joint 
design  the  transmitted  waveform  and  a  good  receiver.  This  problem  has  been  discussed  in  [6]  for  radar  detection. 


3.4.  TIME  DOMAIN  OPTIMUM  TRANSMIT  WAVEFORM  DESIGN 
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Assuming  the  transmitted  pulse  p{t)  (to  be  optimized)  is  confined  to  the  symmetric  time  interval  [ -T/2,  T/2].  The 
energy  of  transmitted  pulse  is  then 


(3.24) 


It  follows  from  detection  theory  that  the  best  receiver  is  still  a  matched  filter  matched  to  the  signal  component  in 
r(t),  the  SNR  after  such  a  matched  filter  is  given  as. 

SNR  =  2  Ey/N0.  (3.25) 


where  Ey  =  /7.^2jp(t)  *  h(t)\2  dt  is  the  received  signal  energy.  The  problem  is  then  reduced  to  find  the  optimum 
p(t)  such  that  Ey  is  maximized,  under  the  constraint  of  fixed  Ep. 

It  has  been  shown  in  [13]  (p.  125)  and  [6]  that  the  optimum  p(t)  can  be  obtained  by  solving  the  following  homoge¬ 
neous  Fredholm  integral  equation 

fT/ 2 

Pn<t>n(t)  =  *(t-  T)<pn(r)  dT ,  (3.26) 

J-T/2 

and  let  p(/)  =  where  is  the  eigenfunctions  corresponding  to  the  maximum  eigenvalue  /{,  and  the  time 
reversal  operator  n(t)  is  defined  as  k(1)  =  h(t)  *  h(—t).  Without  loss  of  generality,  it  has  been  assumed  that 
Po  >  Pi  >  p2>  *  •  *  •  Pulse  waveform  reproduces  itself,  sealed  by  a  constant  po,  when  convolved  with  the  time 
reversal  operator  over  the  interval  [—T/2,  T/2] . 

With  optimum  p{t)  and  c(t),  we  achieve  maximum  SNR 

SNR  =  2pqEp/Nq.  (3.27) 

It  is  worth  noting  that  to  meet  the  requirement  of  optimum  information  transmission  might  result  in  another  transmit 
waveform.  Maximization  SNR  requires  the  transmitter  to  put  as  much  as  possible  energy  to  the  frequency  with  high 
transmission  capability  through  the  channel,  under  the  fixed  transmitted  energy  constraints,  while  ignore  the  other 
frequency  modes  with  small  transmission  capability.  It  is  possible  that  these  smaller  modes  contain  a  significant 
quantity  of  information  useful  for  the  system.  We  consider  an  extreme  case  here.  When  T  —*  oo,  it  has  been  shown 
in  [15]  (p.  360)  that  all  the  energy  should  be  placed  into  one  single  frequency  Snax  in  the  operation  band  for  which 
H(f)  is  the  greatest.  This  leads  to  a  transmitted  waveform  of  a  sine  wave  at  fmax-  However,  it  should  be  noted 
that  a  signal  of  arbitrarily  narrow  bandwidth  conveys  information  at  a  zero  rate.  Therefore,  trade-off  should  be  made 
between  SNR  and  information  transmission. 

To  illustrate  the  principle,  we  study  the  optimum  transmit  waveform  for  the  following  channels.  Consider  practically 
all  the  communication  systems  are  band  limited,  we  artificially  impose  a  band  pass  filter  in  the  system.  The  transfer 
function  of  this  filter  Hf{f)  is  defined  as 


A  <  I/I  <  h 

otherwise 


where  /i  and  f2  arc  the  low  and  high  frequency  of  the  filter,  respectively.  Its  time  domain  impulse  response  can  be 
expressed  as 


hp(t)  =  2f28inc(2f2t)  —  2fi$inc{2fit) 
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The  autocorrelation  of  the  impulse  response  of  this  filter  can  be  written  as 

hr(t)  *  hr(-t)  =  [2f2sinc(2f2t)  -  2fisinc(2fit)\  *  [2f2sinc{-2f2t)  —  2f\sinc(—2fit)] 

=  2  f2$inc(2f'2t)  +  2f\sinc(2  f\t)  —  4f\sinc(2f\t) 

=  hF(t)  (3.28) 

Apparently  autocorrelation  of  the  function  hp{t)  reproduces  itself 

It  should  be  noted  that  for  cognitive  radio  purpose,  the  frequency  response  of  this  filter  could  be  in  an  arbitrary 
passband  (c.g.,  fcth  element  subband)  located  under  the  FCC  spectral  mask  (3.1  GHz  -10,6  GHz  )  [19]. 

Assuming  the  CIR  of  a  specific  channel  hv(t)f  then  the  overall  of  the  CIR  will  be  E{£)  —  hp(t)  *  hc{t). 


Wedge  Channel 


Consider  a  case  that  light  of  sight  (LOS)  is  blocked  by  a  PEC  wedge,  as  illustrated  in  Fig.  3.1,  the  CIR  of  such  a 
channel  without  consideration  of  bandwidth  can  be  expressed  as 


^e(t)  —  \fift 


am  ± 
- 

mr  Bin 


srr^ 


coa  ^  -coe(  ^  ^  )  coo  ; 


■*(0 


(3.29) 


where 


C( 


I 


cos  ~  -cos(^)  cos^-cos(^) 


is  a  constant  whose  value  is  depending  on  the  environment. 

Note  that  is  an  even  function.  The  autocorrelation  of  the  above  CIR  is 

L  ff\  .  jl  (a\  _  **(0  n(0 

c( )  M  0 "  Vi  Vt 


Then  the  autocorrelation  of  the  overall  CIR  for  a  wedge  channel  is  written  as 

«(i)  =  hc(t)  *  hc(—t)  *  hf(t)  *  hp(~t) 

=  C2n  ■  f  \2f-28inc{2f2t) —  2f\sinc(2f\t)\dt 
J—OO 

=  C 2  ■  |Sr(27T/it)  -  5r(27r/2t)l  +  CV/2  (3.30) 


where  sine  integral  function  Si(x)  is  defined  as 


Note  that  Si(x)  can  be  implemented  with  a  Matlab  built-in  function,  substitute  (3.30)  into  (3.26)  and  solve  (3.26) 
for  4>(t)  numerically  we  can  obtain  the  optimum  waveform  for  wedge  channels. 
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Rectangular  Channel 

The  time  domain  OR  for  such  a  channel  has  been  derived  in  Section  33.2.  It  follows  from  Eq.  (3J5)  that  the 
OR  consists  of  three  parts.  Practically,  the  contribution  from  the  third  pari  (mutual  interaction  between  two  edge 
waves)  in  (3.15)  is  very  weak  and  could  be  ignored.  We  will  also  verify  the  validation  of  this  assumption  through 
our  numerical  results  subsequently.  Then  the  OR  can  be  approximated  by  the  first  two  parts 


(Trb^"(i)  * ' 4(1  - r“> +  '  (7 T^T‘(,)  * S(t  ~ Tb) 


(3.31) 


Note  in  (33 1 )  wc  have  ignored  the  common  time  delay  to,  for  simplicity. 

For  the  term  1  most  of  the  energy  is  focused  at  the  singularity  where  /  —  0,  thus  we  have  — — ^-7s  — 

V*t  To.frJ  V  *  U+Ta,b)v* 

1 


1  Then  (3.3 1)  can  be  further  approximated  as 


h(t)  »  Ca ^  *  tf(t  -  t.)  +  Cb^±  *  6{t  -  7i) 


(332) 


whCTe  =  2^  and  Cb  =  5F75- 

The  time  reversal  operator  k{I)  for  such  a  rectangular  channel  can  be  expressed  as 

«(0  »  (<2  +  Cl)  •  \Si(2-Tifft)  -  Si{2*ht)\ 

+CaCb  -  [Si(2itfit)  -  Si{2wf2t)\  *  [«((  -  ra  +  r„)  +  6{t  +  ra  -  rh)]  +  (C„  +  Cb)2n2/2  } 


Two  Half  Planes  Waveguide  Channel 

As  slated  in  Section  33.3,  the  half  planes  waveguide  model  is  a  little  more  complicated  than  the  other  two  models 
studied  above,  due  to  the  exist  of  multiple  replicas  (caused  by  waveguide  phenomena).  However,  we  can  simplify 
the  problem  by  focus  on  the  main  lobe  of  the  k(£)  instead  of  the  whole  k(£).  The  main  lobe  of  for  two  half 
planes  channel  can  be  approximated  as 

*{i)  £=  2 Al  -  [SiftTtfxt)  -  Si(2icf2t)\ 

(£m=i  ™-3)  ’  [b  («(2»At)  -  Si(2nht))  +  &  (SS&M  _  (334> 

For  the  case  of  high  frequency  where  /i  — *  oo  and  f2  — >  oo,  interestingly,  n(t)  is  reduced  to  a  form  of 

k(£)  ^  C*  ■  [(Si&wht)  —  Si( 2tt f2t))\ 

where  Cf  —  2  A  f}  F  -1  m~s)  is  a  constant  depending  on  the  environment  parameters  such  as  the  location 

of  the  receive  and  transmit  antennas(/?,^<p),  and  the  distance  between  the  two  half  planes  (a).  Therefore,  under 
these  assumptions,  the  optimum  transmit  waveform  for  the  two  half  planes  waveguide  channel  will  be  the  same  as  a 
wedge  channel,  which  has  been  solved  above. 
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3,4.2  Numerically  Computing  Eigenfunctions 

Except  for  some  classical  case  such  as  band  limited  signals  [13,  6],  the  eigenfunctions  and  eigenvalues  are  gener¬ 
ally  difficult  to  be  solved  in  a  closed  form.  In  this  section  we  are  interested  in  computing  the  eigenfunctions  and 
eigenvalues  numerically.  The  idea  is  to  discretize  the  function  and  borrow  the  powerful  matrix  theory  to  solve  the 
problem.  We  firstly  discretize  the  autocorrelation  function  «(f )  with  sufficient  low  sampling  rate  such  as  to  include 
all  the  significant  channel  frequencies. 

As  stated  in  the  previous  section,  the  optimal  pulse  waveform  (maximizing  SNR)  is  a  solution  to  the  integral  function 


I*n4>n(t)  = 


(3.35) 


Note  that  the  TR  operator  n(t)  is  an  even  function,  i.e.,  n(t  -  r)  =  «(r  - 1),  the  above  equation  can  be  reformulated 
as 


Pn4>n{t)  = 


r  -  f)^n(r)  dr 


(3.36) 


r-T/2 

Pn  ■  =  /  K(t  -  T)(f>n(t)  dr  n= 

JT/2 


<M-*) 

M-N  + 1) 

k{-2N) 
k(-2N  +  l) 

...  *(0) 

...  k(1) 

M-N) 
M-H  +  i) 

1  <M*) 

*(0) 

...  k(2N) 

M*) 

(3.37) 


(3.38) 


We  then  transform  the  integral  function  to  a  matrix  problem.  To  find  the  solution  of  the  integral  function  (3.35)  is 
then  equivalent  to  finding  the  eigen  function  of  matrix  n{t)(2N  +  1  x  2N  +  1}  in  (3.38).  Several  remark  notes  can 
be  drawn  from  (3.38) 


•  The  optimum  waveform  obtained  via  the  numerical  solution  introduced  above  is  not  necessarily  the  unique 
one;  though  the  maximum  eigenvalue  is  unique,  the  corresponding  eigenfunction  might  not. 

•  If  there  is  no  pulse  distortion  (e.g.,  free  space  transmission),  the  time  reversal  operator  n(t)  =  -  r0), 

where  ao  and  r0  are  the  constants  depending  on  the  distance  between  the  transmitter  and  the  receiver.  For  this 
special  case,  there  is  no  optimum  transmitted  waveform  since  any  waveform  will  reproduce  itself  after  being 
convoluted  with  the  time  reversal  operator. 

•  We  can  improve  computation  accuracy  by  decreasing  sampling  interval  At.  However,  the  computation  com¬ 
plexity  is  increased  correspondingly. 


3*5,  NUMERICAL  RESULTS 
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It  is  also  evident  that  the  availability  of  time  domain  OR  h(t)  (and  thus  its  time  reversal  operator  «(f))  is  critical 
to  computing  the  optimum  waveform  <&(£).  The  following  will  show  how  to  get  the  time  domain  OR  for  several 
canonical  environments. 


3.5  Numerical  Results 


3.5.1  Calculation  of  Singularity  in  the  Impulse  Response  of  Wedge  Diffraction 

A  singularity  appears  in  the  impulse  response  for  most  of  channels  where  diffraction  phenomena  is  involved,  e.g., 
in  Eq.  (3.6)  (wedge  diffraction),  when  t  —*  0,  then  — *  oo,  thus  in  Eq.  (3.4),  there  is  a  singularity  at 

t  ~  0.  How  to  numerically  calculate  this  singularity  is  not  a  trivial  issue.  Since  the  singularity  is  exactly  where 
the  energy  concentrates  and  also  defines  the  type  of  the  pulse  distortion,  it  can  not  be  directly  removed  using  the 
lime  windowing.  From  an  Electromagnetic  point  of  view,  this  singularity  at  the  diffracted  wavefront  is  essential. 
An  integral -differential  ( Int-Diff)  approach  has  been  proposed  and  proved  to  be  an  efficient  way  to  remove  this 
discontinuity  at  t  =  0  during  the  calculation  [3, 1 0],  (n  this  section,  we  will  further  verify  this  approach  by  comparing 
it  with  the  well  known  ATT  approach  [8]. 


I  nt-Diff  Approach 


Mathematically,  the  principle  of  this  Int-Diff  approach  can  be  described  as  follows. 


(339) 


hence,  we  have 


P(0  *  9(t)  = 


dt 


p(t) 


■■  [  9{t)dt 
Jn 


(3 ,40) 


where  g(t)  has  a  singularity  in  the  time  domain.  The  integration /J  g(t)dt  removes  the  singularity  and  can  often  be 
obtained  in  a  closed  form.  The  reason  why  we  call  this  approach  Int-Diff  approach  is  because  additional  integration 
and  differentiation  operation  arc  involved  to  get  the  convolution  result. 


Int-Diff  approach  provides  a  very  efficient  way  to  do  the  numerical  convolution  for  the  case  when  there  is  a  singularity 
in  the  functions  to  be  convolved.  This  approach  is  specially  important  when  we  deal  with  the  diffraction  problem  in 
the  time  domain  because  there  is  always  a  singularity  in  the  diffraction  coefficient  due  to  the  nature  of  diffraction. 

Again,  we  consider  wedge  diffraction  case  as  an  example.  Notice  that  the  singularity  in  (3.4)  is  solely  caused  by  term 
We  will  focus  on  this  term  first.  The  rest  of  work  is  just  linear  summation.  According  to  int-diff  approach, 
firstly  we  will  need  to  integrate  Happily,  we  have  a  closed  form  for  the  integration  of  F(A'mi  £),  which 

can  be  expressed  as 

F(Xm,t)  =  /0l  F{Xm,t)dt 

=  2\f^atan(\fifc) 


(3.41) 
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Let  D*’h{i )  denotes  the  integration  of  the  wedge  diffraction  coefficient  !7'h{t),  then  wc  have 

4 


£>*'h(f)  = 


(3.42) 


Note  that  singularity  in  lP-h{t)  has  been  removed  in  the  new  integration  expression  Ds'h{t).  The  latter  can  be  used 
to  do  the  convolution,  followed  by  a  differentiation  operation  as  a  compensation. 


Analytic  Time  Transform  (ATT) 

ATT  is  another  classic  technique  used  to  perform  an  efficient  convolution  with  a  broad  class  of  excitation  pulse 
functions  int-diff.  In  this  paper,  this  classic  approach  will  be  briefly  introduced  and  serve  as  a  sanity  check  for  the 
int-diff  approach. 

The  second  version  ATT  in  [8]  is  defined  as 

/  (f)  =  -  /°°  dr  ,  for  lm(t)  >  0  (3.43) 

KJ-oot-T 

The  “+”  sign  over  /  means  an  analytical  signal  in  a  complex-time  domain,  Let  p(f),  and  h(t)  denote  the  incident 
pulse  and  the  impulse  response  of  the  system,  respectively,  then  we  have. 


p(t)*h(t)  =  -  Re 


P  {!)*  h.  (t) 


(3.44) 


Notice  that  the  singularity  on  real  time  axis  (/m(()  =  0)  in/i  (£)  is  only  a  branch  point  singularity  and  not  a  pole. 

Therefore,  the  difficulty  in  calculation  the  singularity  in  h(t)  has  been  removed  by  ATT  in  ft  (i).  Moreover,  if  we 
choose  the  incident  pulse  such  that  in  the  frequency  domain  p(tu)  can  be  modeled  using  exponential  function  as 


N 


or  equivalently,  in  time  domain 


p(w)  =  ^  Ane  /or  w  >  0 


+  .  .  +  .  , 
P  (*)  =zlAn  6  (*  +  .?“«) 

Tl— 1 


then  Hie  convolution  can  be  evaluated  in  closed  form  as 


r  n 


p(t)  *  h(t)  —  Re 


^  ^  ^Ti  h  (i  H“ 


Ln=l 


where  the  shifting  properly  of  the  analytic  delta  function 


^  J  (t  +  jan)*h(t)  =h  ( t+jan ) 


(3.45) 


(3.46) 


(3.47) 


(3.48) 


has  been  used. 
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Figure  3.6:  Comparison  between  the  approachs  oflnt-Diff  and  ATT.  In  (a)  amplitudes  have  been  normalized  to  the 
incident  pulse  to  compare  the  power  of  the  incident  pulse  and  the  received  pulse.  In  (b),  the  incident  pulse  has  been 
scaled  to  the  same  amplitude  of  the  received  pulse  to  compare  their  waveforms 

To  verify  our  result  with  int-diff  approach,  we  use  the  same  incident  pulse  as  used  in  [8].  Mathematically,  the 
frequency  domain  of  the  incident  pulse  p(u>)  is  defi  ned  by 

PH  =  c0(l  -  e-uT)p'e~'jPlT  (3.49) 

where  T  =  j,  fc  represents  the  center  frequency,  which  can  be  adjusted  easily  according  to  differ¬ 
ent  requirements.  The  peak  of  the  p(u>)  is  normalized  by  choosing  Q  —  ^  )  •  Time  domain 

waveform  of  the  incident  pulse  can  be  conveniently  obtained  by 


p(f)  =  Re 

p{t) 

—  Re 

.  Pi 

iy  — 

7T  ^>t  + 


gQ  («) 


n— 0 


j(n  +  P2)T 


,Pi 


where  (Vi)  =  n!(ff,n?r 


(3.50) 


Let  D  (i)  denote  the  analytic  impulse  response  of  wedge  diffraction,  then  using  the  convolution  property  in  Table 
2. 1  in  [8],  the  receiving  pulse  will  be 


r(t) 


Re  '  \  P  (t)*  D  ( t ) 

Re  E£o<?0 (n)b(I+j(n+P2)T) 


(3.51) 


Instead  of  integration,  ATT  approach  remove  the  singularity  by  performing  an  analytic  transform.  The  analytic 
transform  of  an  time  domain  function  is  not  so  straightforward,  especially  when  the  function  itself  is  complicated. 
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According  to  [8],  the  analytic  wedge  diffraction  coefficient!)  (f)  can  be  expressed  as 


(3.52) 


where 


+ 

f*  i  0 


(3.53) 


and  the  rest  of  the  parameters  are  the  same  with  that  of  TD-UTD,  presented  in  section  II. 

Fig.  3.6  shows  the  received  waveforms,  calculated  with  ATT  and  Int-Diff  approach,  respectively.  The  incident 
waveform,  marked  as  red  slash  curve,  is  also  plotted,  as  a  reference.  Firstly,  it  can  be  easily  observed  that  the 
received  waveforms,  calculated  with  different  approaches,  are  in  excellent  agreement.  Secondly,  it  can  also  be 
observed  from  Fig.  3.6(a)  that  the  diffracted  pulse  is  much  weaker  than  the  incident  pulse,  due  to  the  huge  energy 
loss  caused  by  diffraction.  To  better  compare  the  waveform  between  incident  pulse  and  the  diffracted  pulse,  we 
manually  scaled  their  amplitudes  to  the  same  value  and  their  waveforms  are  shown  in  Fig.  3.6(b).  It  is  evident 
that  the  diffracted  pulse  waveform  is  much  different  from  the  incident  one,  indicating  a  pulse  distortion  caused  by 
diffraction.  Finally,  the  amplitude  information  (signal  strength)  can  be  predicted  by  traditional  frequency  domain 
UTD,  the  pulse  waveform  information  that  presented  in  this  paper,  however,  can  only  be  provided  by  time  domain 
formulations. 

The  parameters  used  in  our  simulation  are  as  follows:  For  the  incident  pulse,  /(  =  2,  —  1,  and  fc  —  2  GHz.  The 

parameters  for  the  wedge  are:  Po  =  rr/2,  n  —  3/2,  4>  =  5rr/6,  <A  —  tt/3,  and  the  soft  boundary  condition  has  been 
considered. 

3.5.2  Rectangular  Building 

For  all  the  figures  (unless  otherwise  specified)  in  this  paper ,  the  waveform  labeled  “TD”  is  obtained  by  convolution 
of  the  incident  pulse  p(t)  with  the  time  domain  CIRs  derived  in  this  paper,  and  the  waveform  labeled  “FD+IFFT”  is 
obtained  by  applying  inverse  Fast  Fourier  Transform  (FFT)  to  the  product  of  P(f)H (/),  where  P(f)  and  II  ( f )  are 
Fourier  Transform  of  p(f)  and  channel  transfer  function,  respectively.  We  have  verified  our  time  domain  derivation 
by  comparing  the  ‘TD”  result  and  “FD+IFFT”  result. 

In  the  following  numerical  results,  we  choose  the  second  order  derivative  of  Gaussian  pulse  as  the  incident  pulse, 
which  mathematically  is  defined  as 


(3.54) 


where  w  is  the  parameter  controlling  the  width  of  the  pulse  (and  therefore  the  frequency  bandwidth  of  the  transmit 
signal)  and  tc  is  time  shift  to  pul  the  pulse  in  the  middle  of  the  window. 

Fig,  3.7  shows  a  comparison  between  GTD  and  Fresnel-Kirchhof  models.  Here,  we  set  w  -  0. 1  ns  and  ^  =■  0.3 
ns.  The  environment  parameters  used  in  the  simulation  are  as  follows:  4  —  1000m,  dr  —  900m,  a  =  15m,  and 
b  —  10m. 

Several  observations  can  be  made  from  Fig.  3.7. 
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Figure  3.7:  Simulated  pulse  waveforms  after  a  second  order  derivative  of  Gaussian  pulse  pass  through  a  rectangular 
channel  described  in  Fig.  3.4 

•  For  both  GTD  model  and  Fresnel -Kirchhof  Model,  time  domain  results  show  very  good  agreement  with  the 
results  transformed  from  the  frequency  domain  via  IFFT  (labeled  “FD+IFFT”).  This  verifies  our  lime  domain 
derivation. 

•  The  diffracted  waveforms  based  on  GTD  and  Fresnel- Kirchhof  agree  very  well  to  each  other,  except  that 
the  Fresnel -Kirchhof  result  has  an  additional  small  waveform  response.  This  extra  “baby  response”  can  be 
identified  from  its  time  domain  formula  to  be  the  mutual  interaction  (coupling)  term  between  two  edges,  which 
has  not  been  taken  into  account  by  the  GTD-based  model. 

Again,  for  the  convenience  of  comparison,  common  propagation  delays  has  been  neglected. 


3.5J  Two  Parallel  Half  Planes  Model 

As  mentioned  in  Section  3.3.3,  when  a  UWB  pulse  is  incident  on  the  parallel  half  plane  channel,  the  received  signal 
consist  of  three  parts,  with  each  part  including  two  same  (but  separated  in  the  time)  signal  components.  The  relative 
signal  strengths  of  these  three  parts  depend  on  the  value  of  the  coefficient  4)  and  A  i,  which  in  turn  depend  on  the 
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Figure  3.8;  UWB  pulse  diffracted  by  a  parallel  plane  waveguide 


environment  and  the  location  of  the  antenna.  Fig,  3,8  shows  the  received  signal  when  a  second  order  derivative  of 
Gaussian  pulse  is  incident  on  the  parallel  half  plane  channel.  The  incident  pulse  width  parameter  tv  —  2  ns.  The 
geometrical  parameters  for  parallel  planes  are  as  follows:  p  —  20 m,  <j>  -  8/1 —  5/6x,  a  =  2m, 

h  can  be  seen  from  Fig,  3,8  that  the  waveforms  of  the  received  pulses  are  not  second  order  derivative  Gaussian 
anymore.  It  follows  from  Section  3.3.3  that  they  are  actually  “semi-integral”  forms  of  the  second  order  derivative 
Gaussian  pulse.  Again,  the  agreement  between  the  waveform  obtained  by  time  domain  convolution  (labeled  ‘TD”) 
and  frequency  domain  result  (labeled  “FD+IFFT”)  shows  our  time  domain  derivation  is  correct. 

Although  Eq,  (3.22)  shows  infinite  terms,  we  can  see  from  Fig.  3.8  that  these  sums  converge  very  fast. 


3.54  Optimum  Transmit  Waveform  Design 

Fig.  3.9  shows  the  optimum  transmit  waveforms  for  different  channels.  These  waveforms  have  been  obtained 
by  substituting  /*{*)  of  corresponding  channels  into  (3.38)  and  numerically  computing  associated  with  the 
maximum  eigenvalue.  For  the  convenience  of  waveform  comparison,  the  energy  of  the  computed  optimum  transmit 
waveforms  (<£(£))  for  different  channels  have  been  normalized.  The  parameters  used  in  the  simulation  are  as  follows: 
T  —  2ns,  fi  —  3 GHz7  /2  =  bCHz.  Moreover,  in  each  model,  the  environmental  parameters  we  have  used  here 
are  the  same  as  those  we  have  used  for  the  computation  of  distortion  waveforms  in  this  section.  It  has  been  shown 
in  Section  3,4.1  that  the  optimum  transmit  waveform  for  the  two  parallel  half  planes  channel  is  approximately  the 
same  as  wedge  channel,  under  some  circumstances. 

For  the  purpose  of  maximizing  the  SNR  at  the  receiver  side,  all  the  energy  should  be  sent  through  the  optimum 
transmit  waveform,  which  is  the  eigenfunction  associated  with  the  maximum  eigenvalue  of  *(£),  For  the  purpose 
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Figure  3.9:  Optimum  transmit  waveform  design  for  bandlimited  UWB  signal  passing  through  different  channels 

of  flexibility,  we  are  interested  in  transmission  of  a  linear  combination  of  the  eigenfunctions  whose  eigenvalues 
are  close  to  the  maximum  one.  Fig.  3.10  shows  the  first  six  eigen  waveforms  for  the  passband  channel.  Their 
corresponding  eigenvalues  are  1,0000,  0.9999,  0.9979,  0.9976,  0.9621,  and  0.9599.  We  can  also  transmit  these 
waveforms  parallel  to  increase  the  data  rate.  Note  that  the  eigenfunctions  for  the  low  pass  channel  ([0,  /])  is  knows 
as  prolate  spheroidal  wave  functions  [20],  and  has  been  proposed  as  a  novel  UWB  pulse  shaping  technique  in  [21  ]. 


3.6  Conclusion 

For  wideband  signals  such  as  UWB,  by  principle,  time  domain  frame  work  appears  to  be  a  better  frame  work  to  study 
its  propagation.  The  time  domain  impulse  responses  for  a  series  of  scattering  environments  have  been  derived  for 
the  first  time.  Specifically,  several  electromagnetic  communication  fundamental  channels,  including  perfect  conduct 
electromagnetic  wedge,  large  rectangular  plates,  and  parallel  planes  waveguide,  are  studied.  The  receiving  signal, 
calculated  through  convolution  in  the  time  domain  impulse  response,  can  provide  waveform  information,  as  well 
as  the  amplitude  information  (Path  loss).  The  pulse  waveforms  passing  through  these  channels  arc  compared  with 
their  original  waveforms.  The  comparison  show  that  for  wideband  signals,  a  transmit  pulse  could  be  extremely 
distorted  by  diffraction  mechanism.  The  results  converted  from  frequency  domain  by  applying  a  inverse  fast  Fourier 
Transform  (IFFT)  serves  as  sanity  checks  for  our  time  domain  derivation.  Given  the  lime  domain  channel  impulse 
response,  we  also  study  the  optimum  transmit  waveform  such  that  the  signal  to  noise  ratio  (SNR)  at  the  receiver  side 
is  maximized,  under  the  constraint  of  fixed  transmitted  power. 
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Figure  3,10:  The  first  six  eigen  waveforms  for  a  passband  channel 
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Chapter  4 


Optimum  Transmission  Waveform  in  RF 
Harsh  Environments 

4.1  Introduction 


Communications  and  networking  in  confined  metal  environments  are  critical  to  applications  such  as  intra-ship,  intra- 
vehicle,  intra-engine,  manufacturing  plants,  assembly  lines,  nuclear  plants,  body  area  network  sensors  suiroundcd 
by  vehicles  and  tanks,  etc.  Narrow  band  wireless  technologies  have  proved  ineffective  in  these  environments,  due 
to  resonance  caused  by  the  metal  walls.  The  physical  mechanisms  of  short  pulse  propagation  in  a  confined  metal 
environment  are  studied  by  Felsen[l],  from  a  transient  radar  cross  section.  An  Ultra-wideband  (UWB)  system  [2]- 
[5]  is  proposed  here  to  resolve  the  resonance  into  many  time-resolvable  pulses — corresponding  to  extremely  rich 
multipath.  The  energy  of  these  pulses  is  difficult  to  collect  using  traditional  RAKE  or  multicarrier  technologies.  As 
an  alternative,  a  powerful  time  reversal  channel-matching  technique  is  suggested  for  this  purpose. 

Periodically  modulated  radar  returns  from  jet  engines  can  be  modeled  as  a  chirp  [6],  This  view  provides  some  insight 
into  our  problem,  with  measured  channel  pulse  response  from  a  pulse  train  as  long  as  800ns. 

Chirp  UWB  systems  [7] — accepted  into  IEEE  802. 1 5.4a  in  Dec.  2006— have  some  favorable  advantages  such  as 
low  cost,  low  transmission  power  (long  range),  low  battery  power,  anti-jamming,  easy  synchronization,  robustness 
to  Doppler  frequency  shift  (mobility),  etc.  Analog  SAW  filters  can  be  employed  to  perform  receiver  matched  filter. 

This  paper  is  motivated  by  the  UWB  applications  in  the  harsh  environment,  by  adapting  the  existing  standards. 
Time  reversal[9][10) — see  [8]  for  a  review — has  been  used  to  compress  the  channel  response  into  a  sharp  main 
lobe.  Similarly,  chirp  matched  filter  results  in  pulse  compression.  We  design  the  time  reversal  chirp  system  from 
the  classical  modal  modulation  and  give  some  relative  performances  in  different  situations.  Finally,  the  spectral 
efficiencies  of  different  transmission  schemes  are  compared  to  benchmark  the  time  reversal  chirp  scheme. 
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Figure  4. 1 :  Rectangular  metal  cavity  used  for  propagation  experiments. 


4.2  Pulse  Propagation  in  Rectangular  Metal  Cavity 


The  measurement  is  performed  using  frequency  domain  technique  to  analyze  the  characteristics  of  the  UWB  channel 
in  the  confined  metal  environment.  The  rectangular  confined  metal  cavity  in  Fig.  4. 1  measures  1 6  feet  by  8  feet  by 
8  feet.  The  materials  of  the  walls  are  aluminum. 

The  sounding  is  carried  out  by  sweeping  a  set  of  narrowband  sinusoid  signals  (tones)  through  a  wide  frequency 
band.  The  data  of  channel  frequency  response  is  collected  using  a  Vector  Network  Analyzer  (VNA) — Agilent 
N5230A  (300kHz  -  13.5GHz).  The  VNA  sweeps  from  3GHz  to  10GHz  using  7001  points  with  frequency  step 
of  1MHz  .  The  power  of  each  tone  is  lOdBm,  and  the  same  for  all  the  measurements  reported  in  this  paper.  The 
number  of  averaging  is  128,  which  is  used  to  enhance  the  accuracy  of  the  measurement.  Meanwhile,  the  cables  and 
the  connectors  are  calibrated  before  measurement  to  compensate  for  the  frequency  dependent  losses. 

The  measurement  is  done  for  line  of  sight  (LOS)  situation.  The  distance  between  transmitter  antenna  and  receiver 
antenna  varies  from  0.5m  to  4m  with  distance  step  of  0.5m  and  both  antennas  are  placed  at  1 .35m  above  the  ground. 
The  transmitter  antenna  is  fixed,  and  the  receiver  antenna  is  moved  to  different  locations  along  the  middle  line  of 
one  side  of  the  rectangular  metal  cavity. 

After  the  channel  response  is  measured  in  the  frequency  domain  (Fig.  4.2),  zero  padding  is  done  between  0Hz  and 
3GHz,  conjugate  reflection  between  -10GHz  and  0Hz  and  inverse  fast  Fourier  transform  (IFFT)  are  used  to  obtain 
the  channel  pulse  response  ft(f)  in  Fig.  4.3  and  Fig.  2.1(d).  Channel  pulse  response  in  the  office  environment 
(Fig.  2.1(d))  is  very  different  from  that  in  the  rectangular  metal  cavity  (Fig.  4.3),  for  the  transceiver  distance  of  4m. 
Channel  pulse  response  in  the  metal  cavity  is  long  and  strong  for  the  transceiver  distance  of  4m  (Fig.  4.3). 

The  energy  of  the  channel  is 

Eh=  [  *  h2(t)dt  (4.1) 

Jo 

where  T*  is  the  time  duration  of  the  channel.  The  autocorrelation  of  the  channel  is  shown  in  Fig.  4.4,  defined  as 


Rhh(t)=k(t)*h(J'h-t). 


(4.2) 
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Figure  4*2:  Measuring  frequency  response  H(f)  in  a  rectangular  metal  cavity,  A  total  of  7001  equal  frequency  tones 
arc  sent  from  3GHz  to  10GHz  (with  a  step  of  I  MHz). 


x  10  4  4  meters  (nwtal  cavity) 


Figure  4.3:  Channel  pulse  propagation  h(t)  in  a  rectangular  metal  cavity.  The  pulse  has  7001  equal  frequency 
components  (with  a  step  of  1  MHz)  from  3GHz  to  10GHz. 

It  is  interesting  to  notice  that  the  total  energy  of  the  channel  pulse  response  in  a  confined  metal  cavity  drops  very 
slowly  as  compared  to  office  and  hallway  environments,  as  shown  in  Fig.  4.5. 


4.3  Time  Reversal  Transmission  with  Chirp  Signaling 

4.3.1  Modal  Modulation 

The  starting  point  of  our  system  formulation  is  from  the  well  known  result  of  Holsinger  (1964)  [13],  regarding 
optimum  detectors  for  a  linear  time-invariant  (LTI)  channel  with  memory,  as  shown  in  Fig.  4.6.  Ifx(f)  is  the  input  to 
the  LTI  channel  with  channel  pulse  response  h(t)  mentioned  before.  The  symbol  duration  is  '1$,  and  the  waveform 
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Figure  4.4:  Autocorrelation  of  channel  pulse  response  —  h(f)  *  E{Th  ~  t)/Eh  in  a  rectangular  metal  cavity. 


Figure  4.5:  The  total  energy  of  the  channel  pulse  response  in  a  confined  metal  cavity  decays  very  slowly,  as  compared 
with  office  and  hallway. 

duration  of  r(i)  is  T.  The  output  waveform  is 

r  (t)  =  f  x  (r)h  ( i  —  r)  dr,  0  <  t  <  Ti  (4.3) 

Jo 

for  the  observation  time  [0,7)j.  White  Gaussian  noise  is  assumed.  According  to  the  channel  pulse  response,  the 
observation  lime  and  waveform  duration,  a  symmetric  function  K ((,  s)  can  be  defined  as 

K  (t,  s)  =  [  h(o  —  t)h(cr  —  s)  da,  0  <  t,s  <  7'  (4.4) 

Jo 

Furthermore,  we  can  define  a  possibly  infinite-size  set  of  orthonormal  eigenfunctions,  #,(£),  that  satisfy  the  relation 

I13j 

An  ■  <Pn{t)  =  f  K(t,T)ipn{r)dT,  0  <t<T 

Jo 


(4.5) 
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Figure  4.6:  Modal  modulation  is  optimum  for  the  LTT  channel  with  memory.  When  a  self-reproducing  eigenfunction 
ipn (f )  is  transmitted,  after  the  channel  matched  filter  h(mm(Th,T\)  -  t),  the  same  waveform  v?n(0  (within  a  gain 
factor  A„)  is  obtained.  The  optimal  detection  waveform  is  matched  to  the  <#,(£). 


T;  tT*T. 

Figure  4.7:  Time  reversal  modal  modulation  is  proposed  to  simplify  the  transceiver  complexity.  The  channel  matched 
filter  h(Th  —  t)  is  moved  to  the  transmitter  side.  The  self-reproducing  4»(0  can  be  defined,  similarly  to  the  modal 
modulation  illustrated  in  Fig.  4.6. 


A  special  form  of  (4.5)  is  given  by  [14].  That  is  a  function  that  when  convolved  with  K  (t,  s )  over  the  interval 
[0,7  ]  reproduces  itself,  scaled  by  a  constant  called  eigenvalue. 

^n(/)  n  =  1, 2, . . . ,  oc  can  be  seen  as  a  set  of  basis  functions,  so  x(t)  can  be  written  as 


OO 


x{t)  =  o  <  t  <  T 

(4.6) 

n=0 

and 

f T 

xn=  x  ( t )  ipn  ( t)dt ,  n  =s  1, 2, . . . ,  00 

J  o 

(4.7) 

so  the  energy  of  the  transmitted  signal  is 

fT 

Et  =  /  |x  (f)|2dt 

Jo 

(4.8) 

OO 

=  £4 

n=  1 

(4.9) 

From  detection  point  of  view,  after  matched  filter,  the  output  signal  to  noise  ratio  (SNR)  is  2 Ej JN6»  where 

n=  1 


(4.10) 
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and  No/ 2  is  the  double  sided  noise  spectral  density.  So  the  probability  of  error  for  the  detector  is  given  by 


1  /■* 

\Z2tt  J ^2e/n„  rK  2 


exp(— ]-x2)dx 


(4.  II) 


In  order  to  minimize  Pe,  SNR  should  be  maximized.  Since  by  conversion  A|  >  A2  >  A3  >  *  •  ■ ,  if  the  energy  of 
x(t)  is  fixed,  maximizing  E  means  x(f)  =  xi<p\(t)  which  is  the  optimal  waveform  for  detection.  So,  the  maximum 
SNR  can  be  expressed  as. 


SNR.max 


2  A 

“nT 


2Amax  Et 


No 


(4.12) 

(4.13) 


From  capacity  point  of  view,  the  modal  modulation  system  have  many  parallel  orthogonal  channels  and  each  channel 
corresponds  to  one  eigenfunction  #i(£).  We  can  use  the  classical  water  filling  scheme  to  calculation  coefficients 
xn  n  =  1,2, . . .  ,00  based  on  the  energy  constraint  of  transmitted  signal  x[t). 

The  observation  time  T\  and  the  symbol  duration  Tg  are  very  important  for  the  system  design.  If  7J  <  7’  +  7),,  then 
Ts  should  be  equal  to  or  greater  than  max  (2i  +  nun(7i,7/l),7'  +  7),).  And  if  71  >  7  +  Th,  then  K(t,s)  reduces 
to  K(t  —  s)  and  Tg  should  be  equal  to  or  greater  than  7'  +  27A.  In  these  situations,  the  generalized  Nyquisl  criterion 
can  be  satisfied,  so  we  do  not  need  to  consider  ISI  when  using  modal  modulation. 


4.3.2  Time  Reversal  Modal  Modulation 


Time  reversal  modal  modulation  is  proposed  to  simplify  the  transceiver  complexity.  The  h(%  —  f),  time  reversal 
filter  is  moved  to  the  transmitter  side.  Similarly  to  the  modal  modulation  illustrated  in  Fig.  4.6,  the  self-reproducing 
eigenfunctions  <&,(()  with  eigenvalues  illustrated  in  Fig.  4.7  can  be  defined  as 


where 


«n 


Khh(t  ~  r)(f>n{T)dT,  0  <  l  <  T 


(4.14) 


Rhh{t)=k(t)*h(-t)  (4.15) 

Notice  that  when  T\  >  7’  +  7),  and  Tg  >  T  +  27h,  moving  the  channel  matched  filter  h('l'h  -  1)  in  Fig.  4.6  to  the 
transmitter  side  will  not  affect  the  optimality  of  the  system.  In  other  words,  the  systems  in  Fig.  4.6  and  Fig.  4.7 
arc  equivalent  in  the  above  situation.  In  Appendix  B,  we  give  the  proof  that  when  the  optimal  waveform  is  used  for 
detection,  the  maximum  SNR  after  matched  filter  is  equal  to  2tirnttxEt/N()  similarly  as  shown  in  Eq.  4. 1 3. 


4.3  J  Time  Reversal  Chirp  System 

Time  reversal  modal  modulation  gives  us  a  lot  of  insights  to  design  UWB  system.  The  main  work  shifts  from 
how  to  handle  ISI  to  how  to  design  the  transmitted  waveform.  We  can  calculate  the  orthonormal  eigenfunctions 
<i>n(t)  n  =  1,2, . . .  ,00  according  to  the  autocorrelation  of  the  memory  channel.  Each  eigenfunction  can  be  seen  as 
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Figure  4.8:  Output  waveform  of  the  chirp  matched  filter  without  time  reversal,  defined  as  in  (4,17). 


Figure  4.9:  Output  waveform  of  the  chirp  matched  filter  with  time  reversal. 


a  waveform  and  any  linear  combinations  of  these  waveforms  can  be  used  in  different  situations  and  requirements  as 
well. 

It  is  well  known  that  chirp  signal  has  the  sharp  autocorrelation  characteristic  and  it  is  widely  used  in  the  radar  field. 
We  can  easily  use  chirp  signal  as  the  waveform  instead  of  the  eigenfunction  to  get  the  time  reversal  chirp  system. 


Chirp  signal  is  defined  as 


ch  (t)  =  cos  (21 T  ^ fat  ±  ,  0  <  i  < 


(4.16) 


where  /o,  7  ",  /i  arc  center  frequency,  signal  duration  and  sweep  rate,  respectively.  The  amplitude  spectrum  of  a  chirp 
signal  approaches  to  a  band  pass  filter  as  the  time  bandwidth  product  of  it  becomes  large. 


The  frequency  of  a  chirp  sweeps  continuously,  and  its  duration  can  be  changed  independently.  The  bandwidth  B  of 
a  chirp  signal  is  1 i  —  f.t  *  7',  Signal  with  +fi  is  an  up  chirp  while  —ft  is  a  down  chirp.  The  pair  is  matched  to  each 
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Figure  4.10:  Simulation  result  for  time  reversal  chirp  system  in  the  rectangular  metal  cavity. 

other.  The  output  ofthe  matched  filtering  is  Red t)  —  Ch(t)  *  ch{'T  -  t) — implemented  using  surface  acoustic  wave 
(SAW)  devices — expressed  by  [19] 

Rcc(t)  »  v^cos  (27T/0(t  -  T))  Sm  >  0  <  t  <  21'  (4.17) 

which  is  illustrated  by  Fig.  5.3. 

if  time  reversal  is  not  considered,  after  the  chirp  matched  fitter,  a  compressed  pulse  that  has  a  short  duration  of  2/ B 
is  generated.  The  peak  in  Fig.  5.3  is  >/TB.  As  a  result  of  matched  filter,  a  fiat  envelope  pulse  with  duration  7’  is 
compressed  into  a  much  shorter  waveform  of  duration  2/S.  If  time  reversal  is  used,  the  output  waveform  of  the 
chirp  matched  filler  is  shown  in  Fig.  4.9. 


4.3,4  System  Performance  Comparison 

The  system  performance  comparison  between  time  reversal  modal  modulation  system  and  time  reversal  chirp  system 
will  be  considered  here  when  Ts  >T  +  27\. 


From  Detection  Point  of  View 

From  detection  point  of  view,  SNR  after  matched  filter  in  the  receiver  side  is  the  metric  to  determine  the  system 
performance. 

The  maximum  SNR  for  time  reversal  modal  modulation  system  can  be  shown  in  Eq.  4.1 08. 

For  time  reversal  chirp  system,  if  chirp  signal  q,(f)  is  used,  then, 

00 

Ch  {t)  =  J2cn<i>n(t),0<t<T 

n—  I 


(4.18) 
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where 


Cn=  f  Ch{t)4>n(t)dt,  n  =  1,2,... , 


oo 


and  <j>n  (0  is  from  Eq.  4J4. 

The  energy  of  the  transmitted  signal  is 


Et  = 


=  J  jf  ^^c„0n(T)  j  /i{r  —  t)  rfT-  jT  ^^c„^n(s)j /i(s -i) 

-  J  J  J  h(r  -t)h{s~t)dt  c«<An  («)  MS  Cn<^n  (T)j 

=  jfjf  flft/i  (T  —  a)  (s)^  cn<An 

rT  /  <*  fT  \  (  00 

-  jo  ( H  ^  ^  fl«>  tT  _  *)  ^  (*) ds J  ( £ 

<®  -r 

=  Jo  \<j>n  (r)|2  dr 

=  £*n<£ 


dsdt 


dsdr 


,  ^n<^n  (5)  j  ^ „  cn0n  (r)  j  dsdr 
Khh  (T  -  s)  <(>n  (s)  ds'j  ( y  Cn<t>n  (r)  ]  dr 
I  dr 


In  the  receiver  side,  after  matched  filter,  SNR  is 


SNR^hjrp  = 


[j„r  ( j:  to)  (g,**-  w) dl] 

f /oI'(£j^w)  * 

(£q  («"<£)) 


OO 

& 


n=€ 


2£? 


oo 

n^O 


(4.19) 


(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 


(4.28) 


(4.29) 


(4.30) 
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So,  the  maximum  SNBgap  between  time  reversal  modal  modulation  system  and  lime  reversal  chirp  system  is 


SSNK  — 


*ma x&t  £ 
n=0 

w 

00 

^max 

oo 

T. 

n— 1 


(4.31) 


(4.32) 


If  T  goes  to  infinity,  Cn  can  be  treated  equally  and  Eq.  4.32  is  reduced  to 


ffSNR  = 


(max{|  H  (/)1})2 

hrl£\nw?« 


(4.33) 


which  can  be  treated  as  the  peak  to  average  ratio  from  perspective  of  frequency  domain. 

When  T  is  equal  to  Ins,  10ns,  100ns  and  goes  to  infinity,  the  maximum  SNR  gap  will  become  2.5dB,  4.5dB,  8.4dB 
and  !2.7dB,  respectively.  The  result  is  based  on  the  measurement  data  shown  in  Fig.  4.2. 


From  Data  Rate  Point  of  View 

If  L'levcl  pulse  amplitude  modulation  (PAM)  is  considered  and  there  are  L  equally  likely  amplitude  values  placed 

at  level  ±d,  ±3 d . ±(L  -  I)d.  In  lime  reversal  system,  the  waveform  filter  in  the  transmitter  side  is  Or{t)  and 

the  corresponding  frequency  response  is  Gr(f):  the  time  reversal  filter  is  h(-t)  and  the  corresponding  frequency 
response  is  //*(/).  So,  the  average  transmitted  power  is 

n  =  i4M) 


In  the  receiver  side,  after  matched  filter,  the  minimum  distance  of  the  signal  is 

and  the  noise  power  is 

°2  =  t  /_JGt(/)|2^ 

where  jV0/2  is  the  double  sided  PSD  of  the  white  Gaussian  noise.  Thus,  the  PAM  probability 


(4.35) 

(4.36) 

of  symbol  error  is 


P, 


3  zrartJZ\GTU)\2\inf)\24f 

U- 1  *0  \GT(f)?df 


(4.37) 

(4.38) 
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6! 


where  Q{x)  is  defined  as  Q{x)  =  e  t2^dt  and  the  corresponding  bit  enor  rate  (BER)  1%  can  be  determined 

as  follows  if  L  is  given, 

n  =  r-^r  (4.39) 

log2  L 

Meanwhile  we  assume  that  the  time  duration  of  Gr(f.)  is  T  and  T$  =T  +  2 7},.  So,  the  number  of  bits  that  can  be 
transmitted  in  each  Ts  is 

b  =  log2  L  (4.40) 

and  L  can  be  equal  to  2, 2s,  23, _ 

In  the  time  reversal  chirp  system,  Gr(t)  is  c*(£)  defined  in  Eq.  4.18.  If  T,  I\  and  Lc  are  given,  we  can  get  2T„I\/ No 
from  Eq.  4.37  and  Eq.  4.39  and  the  data  rate  of  the  time  reversal  chirp  system  is 


logs  Lc 
TS 


(4.41) 


In  the  time  reversal  modal  modulation  system,  there  are  many  parallel  orthogonal  channels  and  each  channel  corre¬ 
sponds  to  the  orthonormal  eigenfunction  (t)  and  eigenvalue  from  Eq.  4.14.  We  assume  >  k2  >  >  •  ■  • . 

Meanwhile,  the  bits  transmitted  in  these  different  channels  are  assumed  to  be  totally  independent. 

For  the  nth  channel,  the  PAM  level  is  Ln‘,  Gjn(t)  is  xn<fi„  (t)  and  x„  is  the  value  to  be  determined  later;  the 
transmitted  power  is  Ptn ;  BER  is  Ptm',  the  PAM  probability  of  symbol  CTtor  is  Pen\  the  number  of  bits  transmitted  in 
each  T$  is  bn.  So,  from  Eq.  4.34,  it  follows  that, 

f_Ll-  1  /■“ 

Ts  “  j  -00 

kn^(/)i2l//(/)|2^/ 


Ptn  = 


PjGTn(f)\2\H  (/) Is 


(if 


Ll- 1 


Ts  3 
(fi  (Lj  ~  l)KnX3 
Ts  3 


/oo 

l&»  (/}|! 

OO 


df 


and  from  Eq.  4,38,  it  follows  that. 
Pen 


-  2('-£)q{ 
■  2('-^)Q(' 


3  2 T'Ptn  J 

r_“  |GTn(/)!2|W(/)|2d/\ 

LI-  !  No 

fZo\GTn(f)\2df  ) 

3  2T3Ptnroo\^4>n(f)\i\H(f)\2df^ 

£ 

1 

fZe\*n4>n(f)\2df  j 

3  2'/,Ptn/°°00|^(/)|2|W(/)|2d/N\ 

3K> 

1 

£ 

1 

SZolMDfdf  ) 

(4.42) 

(4.43) 

(4.44) 

(4.45) 

(4.46) 

(4.47) 

(4.48) 

(4.49) 

(4.50) 
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and 

u  _  Ren 
bn‘  log  2Ln 

(4.51) 

and 

bn  — 

(4.52) 

Thus,  the  total  power  Pt  used  for  time  reversal  modal  modulation  is 

pt 

oc 

=  ^Ptn 

—  _  1 

(4.53) 

fl=l 

(P  y,  (L*  ~  1  )  «„Ij^ 

-  5*2-  3 

(4.54) 

and  the  total  bits  transmitted  in  each  1$  is 

*  = 

(4.55) 

n—l 

oo 

=  5]log2Ln 

n—l 

(4.56) 

So,  the  data  rate  of  lime  reversal  modal  modulation  system  is 

—  FfT 

(4.57) 

From  the  data  rate  point  of  view,  we  have  the  optimization  problem  shown  as  below, 


<30 


max 

s.t 


b= 

n= 1 

Pt=E  Pn 


(4.58) 


Pbn  —  constant 

The  data  rate  of  time  reversal  modal  modulation  system  is  b/'lg. 

The  iterative  loading  algorithm  shown  in  Appendix  C  is  used  here  to  solve  the  above  the  optimization  problem.  So, 
the  data  rate  gap  is  defined  as, 


Rm 

Rc 

OG 

£  log2  Rn 

n— 1 _ _ 

log2  Lc. 


(4.59) 

(4.60) 


When  BER  Pb  is  set  to  10“4  and  the  same  pulse  amplitude  modulation  (PAM)  constellation  is  considered,  the  data 
rate  gap  between  time  reversal  modal  modulation  system  and  time  reversal  chirp  system  is  calculated  as  follows: 
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When  Lc  -  2,  grate  =  0,3.0, 7.0  (dB>  for  T  =  1, 10, 100  (ns);  when  Lc  =  4,  gTOie  =  5.4, 7.4, 10.2  (dB)  for 
T=  1,  10, 100  (ns);  when  Lc  =  8,  grale  =  7.5, 11.1, 13.9  (dB)  for  7'  =  1,10, 100  (ns). 

From  the  comparison,  we  know  the  performances  of  time  reversal  modal  modulation  system  are  better  than  those 
of  time  reversal  chirp  system  when  Ts  >  T  +  27V  If  we  want  to  fully  employ  the  advantage  of  chirp  signal,  the 
requirement  of  Ts  should  be  relaxed.  In  this  case,  the  introduction  of  ISI  is  inevitable.  In  Section  4.3.5,  we  will 
consider  the  system  performance  from  BER  point  of  view  when  Ts  is  less  than  T. 

4.3.5  BER  Performance  for  Time  Reversal  Chirp  System 

In  this  section  we  consider  the  time  reversal  chirp  system  with  the  existence  of  ISI.  In  this  case,  duration  of  symbol  £ 
is  far  less  than  7  ’.  Supposing  PAM  signals  are  used  in  a  single-user  system,  the  transmitted  signal  can  be  expressed 
as 

OO 

$(t)  =  51 ~  nTs)  (4.61) 

n— O 

where 

p(t)  =  Ch(t)  *  h{-t )  (4.62) 

l„  represents  the  nth  discrete  information  symbol  with  duration  'Is,  c^(t)  is  the  chirp  signal  with  the  bandwidth  W 
and  the  time  duration  T,  and  h(—t)  is  the  time  reversed  channel  impulse  response  h(£).  The  received  signal  is 

OO 

r(t)  =  J2  -  nTs)  +  n(i)  (4.63) 

T*=0 


where 

n(t)  is  AWGN. 


y{t)  =  p(f)  *  h{t) 


(4.64) 


The  received  signal  r(f)  passes  the  matched  filter  followed  by  a  sampler  with  a  sampling  rate  of  1  /7s ■  The  impulse 
response  of  the  matched  filter  is  <},(-£).  So  the  output  of  the  matched  filter  is  expressed  as 


OO 

9(0  =  E  URpp{t  -  nTs)  +  v(t)  (4.65) 

n=0 

where  /ipp{£)  is  the  autocorrelation  of  p(£)  and  u(i)  is  the  response  of  the  matched  filter  for  AWGN  noise  n(f). 
Denote 

f^k-n  —  Rpp(t  “  ™7s)|t=fcTs  (4.66) 

Vk  =  (4>67) 

wc  get 

OO 

9Jt  =  5^  InRk-n  +  Vk  (4.68) 

n={J 

OO 

=  Ftoh  +  IfiRk-n  + 

n=flTn/fe 


(4,69) 


64 


CHAPTER  4.  OPTIMUM  TRANSMISSION  WAVEFORM  IN  RF  HARSH  ENVIRONMENT 


Figure  4.1!:  Spectral  efficiency  in  a  rectangular  metal  cavity. 

where  lk  term  represents  the  expected  information  symbol  of  the  k(h  sampling  period,  the  2^  term  is  ISI  and  vk  is 
the  additive  Gaussian  variable  at  the  kth  sampling  point. 

Here  we  performed  the  time  reversal  chirp  simulation  using  the  channel  depicted  in  Fig.  4.3.  The  time  reversal 
chirp  system  has  2-level  PAM  and  500Mbps  data  symbol  rate.  The  frequency  range  of  the  chirp  signal  is  from  3GHz 
to  10GHz  and  the  duration  is  100ns.  Perfect  synchronization  was  assumed.  The  simulation  result  is  depicted  in 
Fig.  4.10.  From  the  BER  curve  we  can  see  that  the  Eb/Ny  penalty  for  10~3  is  approximately  l7dB.  As  fy/N0 
reaches  25dB,  the  trend  of  the  slope  becomes  flat  and  ki,/N0  approaches  between  10“ 4  and  10“5.  For  a  large  class 
of  UWB  applications,  it’s  sufficient  for  BER  at  10-3  with  data  rate  up  to  500MHz  in  harsh  environment  with  I  us 
delay  spread. 


4.4  Capacity  in  Rectangular  Metal  Cavity 

The  channel  transfer  function  is  H(f)  with  bandwidth  W  -  J\-  f0  where  /()  is  the  starting  frequency  and  f\  the  end 
frequency.  Here,  11(f)  contains  path  loss  and  multipath.  The  functional  relation  between  path  loss  and  multipath  is 
complicated  to  parameterize.  Instead,  we  use  the  measured  H(f)  for  calculation.  No  fading  is  observed  for  a  fixed 
transceiver  pair. 

In  the  transmitter  side,  different  precod ing  schemes  are  employed  to  shape  the  spectrum  of  the  signal.  Here,  water 
filling,  time  reversal,  channel  inverse,  and  constant  power  spectrum  density  (PSD)  of  signal  arc  considered.  See 
Appendix  A  for  details.  Their  relative  spectral  efficiencies  are  compared  in  Fig.  4.1 1.  By  simulation  it  is  easy  to 
find  that  the  curve  of  constant  PSD  can  represents  that  of  ideal  chirp  signal,  while  spectral  efficiency  of  time  reversal 
is  almost  the  same  as  that  of  ideal  time  reversal  chirp  signal.  The  channel  inverse  scheme  is  good  for  detection,  but 
performs  poorly  in  terms  of  capacity. 

The  capacity  is  limited  by  the  PSD  of  the  transmission  signal .  In  the  US,  the  FCC  requires  that  the  maximum  PSD 
be  — 41.3dBm/MHz.  If  the  PSD  of  the  AWGN  noise  is  -114dBm/MHz,  then  the  ratio  of  the  transmitted  signal 
power  and  the  received  noise  power,  called  TX  SNR  (defined  in  Appendix  A)  in  Fig.  4,1  J,  is  around  70dB. 
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Appendix  A  Spectral  Efficiency 

The  PSD  of  the  signal  in  the  transmitter  side  and  noise  the  receiver  side  are,  respectively,  denoted  as  S(f)  and^f*. 
The  power  of  the  signal  is  defined  as 

r  =  [  1  s(/)df 

The  power  of  noise  in  the  receiver  side  is  expressed  as 


(4.70) 


N  =  f(h-M  =  ^w 

with  ratio  of  the  transmitted  signal  power  to  received  noise  power  (TX  SNR) 


P  N 


(4.71) 


(4.72) 


The  capacity  is 

H 

fij  2 

(4.73) 

and  the  spectral  efficiency  is 

c 

^logj(l  +  ^)<f/ 

(4.74) 

w  ’ 

h  ~  /.» 

If  wateT  filling  is  used,  then 


where  (jt)  +  —  max[0,  x],  the  constant  /i 


<?(/)  =  („- 


l)  + 


|H(/)P 

is  the  water  level  chosen  to  satisfy  the  power  constraint  with  equality 


(4.75) 


and  the  spectral  efficiency  in  this  case  is 

c  /£(' 
w~  /,-/<> 

When  lime  reversal  is  used,  it  follows  that 

the  constant  a  is  the  factor  chosen  to  satisfy  the  poweT  constraint  with  equality 


(4.76) 

(4.77) 

(4.78) 

(4.79) 

(4.80) 

(4.81) 
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and 


so 


The  spectra]  efficiency  in  this  case  is 


C_ 

W 


_  P 

“  ft  \nu)M 

s(f)  ==  fMDl 

ft  \H(f)\2df 


l  jfa 

fi  ~  fo 


emmpp 


w 


ft -fa 


For  channel  inverse  transmission  scheme,  it  follows  that 

S{/)  =  |H(7jF 

where  the  constant  a  is  the  factor  chosen  to  satisfy  the  power  constraint  with  equality 


with 


Thus, 


S(f)  =  —f - - - - 

K  mbw«  l»('>l2 


and  the  spectral  efficiency  in  this  case  is 


C_ 

W 


ft  log2  (1  +  w,  r  |  —  )df 

2  Jfa  \Hin\2^ 


ft  log2(l  + 


fl-fo 

pW 


w 


log2  (  l  + 


/i-/o 

pW 

ft  wuwdf 


(4.82) 

(4.83) 

(4.84) 

(4.85) 

(4.86) 

(4.87) 

(4.88) 

(4.89) 

(4.90) 

(4.91) 

(4.92) 

(4.93) 


(4.94) 
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If  the  transmitted  signal  is  of  the  constant  PSD  from  Jo  to  f\. 


-  w 


(4.95) 


so  the  spectral  efficiency  is 


w 


h  -  h 


//*  l°g2  ((  +  p\H(f)f2)df 
fi  ~  fo 


(4.96) 

(4.97) 


Appendix  B  Maximum  SNR  after  Matched  Filter  for  Time  Reversal  Modal  Modulation 

From  Eq.  4.14,  we  assume  the  maximum  «„  is  itm>.  and  the  corresponding  orthonormal  eigenfunction  is  0((),  0  < 

t  <  T. 

If  the  transmitted  signal  is  x<p  (t)  *  h  (— f),  then  the  energy  of  the  transmitted  signal  is 


Et  ~  j  [ar^(f)  *h(-t)]2df 

—  x2  J  J  0(r)  h(r  —  t)dT  J  <p(3)h(s-t)dadt 

—  x2  f  f  f  h(r  —  t)h(s  —  t)  dt<f>  (,<?)  0  (r)dsdr 

J  o  Ju  J 

=  X2  [  f  Rhh  (t  -  s)<f>  ( s )  0  (r)  dsdr 

Jo  J 0 

=  X2  /  Kmax<p(T)<f>(T)dT 
JO 

'  [T\4>(T)fdT 
JO 


—  %jax  £ 


(4.98) 

(4.99) 

(4.100) 

(4.101) 

(4.102) 

(4.103) 

(4.104) 


In  the  receiver  side,  after  matched  filter,  the  maximum  SNR  is 


SNR.n 


j/(,T  Kmaxx4>  ( t )  {x4>  (())  dfj 

fjfMON 

2  (Km**xf 
No 

No 

2^man  N 
No 


(4.105) 

(4.106) 

(4.107) 

(4.108) 
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Appendix  C  Iterative  Loading  Algorithm 

The  iterative  loading  algorithm  is  employed  here  to  solve  the  optimization  problem.  The  algorithm  can  be  expressed 

as: 

1  Initialization:  bn  =  0,  n  =  1,2, ...;Ln  =  I,  n  =  1,2, Ptn  ~  0,  n  =  1,2, P^n  =  constant,  n  =  t,2, _ 

2  For  each  channel  n,  if  one  more  bit  is  loaded,  (he  incremental  power  is  calculated  as  follows, 

A  Ptn  =  Ptn  (6„  +  1)  -  Ptn  { bn )  (4. 1 09) 

where  Ptn  (6„)  satisfies 


Pbn  (logs  Ln)  =  2  1- 


( I  1  ^  £>  (  I  3  2T„Ptn  (6„)  Kn  \ 


(4.1 10) 


and  Pm  (6„  +  1 )  satisfies 


„  ,  ,  r  s  nft  '  I  3  7lsPln  (bn  +  \)^\ 


3  We  get  the  channel  index  j  according  to 


j  —  arg  min  A  Ptn 


(4.111) 


(4.112) 


If  Ptn  +  A Ptj  <  Pt,  bj  +  1  —*  bj,  2Lj  —*  Lj,  Ptj  +  A/ tj  — ►  Ptj  and  go  to  step  2;  otherwise,  go  to  step  4. 


n=l 


4  The  number  of  total  bits  transmitted  is  b  =  bn. 

n-  3 
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Chapter  5 


Chirp  UWB  Systems 


Chirp  UWB — accepted  into  IEEE  802.15.4a  in  December  2006 — is  a  promising  technology.  Our  interest  mainly  lies 
in  solving  the  problems  faced  with  impulse  radio — narrowband  interference  and  synchronization.  Range  extension 
is  our  goal. 


5.1  Chirp  UWB 

5.1.1  Challenges  for  IR-UWB 

Though  IR-UWB  has  the  advantages  of  low  complexity,  low  power,  etc.,  it  is  facing  two  challenges,  limited  range 
and  Narrow  Band  Interference  (NBI). 

The  range  for  IR-UWB  communications  is  limited  by  the  transmit  power.  At  present,  typical  range  is  less  than 
30m.  In  order  to  extend  the  range,  an  IR-UWB  system  has  to  increase  its  transmit  power.  However,  UWB  system 
is  a  power-limited  system.  Since  IR-UWB  pulse  has  ultra-short  period  while  ultra-low  duty  cycle,  increasing  the 
transmit  power  will  increase  the  pulse’s  Peak -to- Average  Power  Ratio  (PAPR)  largely.  Such  high  PAPR  makes  it 
difficult  to  implement  in  hardware. 

The  second  challenge  for  UWB  is  narrow  band  interference  (NBI).  Since  UWB  covers  a  large  bandwidth,  it  is  easy 
to  suffer  NBI  with  a  Signal  to  Interference  power  Ratio  (SIR)  less  than  -20dB.  NBI  suppression  algorithms  have  to 
be  used  but  the  complexity  of  such  a  system  will  be  very  high. 

Timing  and  synchronization  is  the  classical  problem  in  an  impulse  radio,  even  if  non-coherent  transmitted  reference 
and  energy  detection  is  used. 


5.1.2  Range  Extension 

In  this  report  we  would  like  to  propose  another  UWB  system  that  has  the  potential  to  overcome  the  above  challenges. 
It’s  Chirp  UWB,  a  chirp  spread  spectrum  (CSS)  based  UWB  system.  CSS  is  originated  from  the  RADAR  system’s 
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pulse  compression  theory.  Pulse  compression  is  developed  to  solve  the  RADAR  system’s  paradox:  maximizing  the 
detection  range  while  minimizing  the  resolution  of  the  detected  object.  Pulse  compression  requires  signals  that  can 
be  simultaneously  spread  in  time  and  frequency.  Spreading  in  time  keeps  the  transmit  power  low  while  accumulates 
high  SNR  after  pulse  compression  at  the  receiver.  Spreading  in  frequency  ensures  high  resolution  for  the  detected 
object.  Chirp  signal  inherently  has  this  dual  spreading  nature.  A  typical  chirp  signal  has  the  form  of  Equation  5.1. 

Ch  (t)  =  COS  ^2jr  (jot  ,  -y  <  t  <  y  (5.1) 

fo,  Tc,  P  are  center  frequency,  signal  duration  and  sweep  rate,  respectively.  From  Figure  5.1,  wc  can  sec  that  chirp 
signal  is  a  linear  frequency  modulation  (LFM)  signal.  Figure  5.2  is  its  amplitude  spectrum.  Its  frequency  'sweeps’ 
continuously  and  its  duration  can  be  changed  independently.  Signal  with  ’+(*’  is  an  up  chirp  while  '-ft'  is  a  down 
chirp,  They  are  matched  to  each  other.  Equation  5.2  and  Equation  5.3  are  the  match  filter’s  impulse  response  and 
matched  output  for  Equation  5.1,  respectively. 


MO  =  v^cos  (j*  .  ~y 


27 Tfit 


2  “  “  2 

(5.2) 

~TC  <  t  <  Tc, 

(5.3) 

Figure  5.3  illustrates  the  compressed  pulse.  After  match  filtering,  a  compressed  pulse  is  generated,  which  has  a  short 
duration  of  2  j B  and  a  peak  of  \/TB.  As  a  result,  a  flat  envelope  pulse  with  duration  %  is  compressed  after  matched 
filtering.  Note  that  the  TB  product  is  the  processing  gain,  which  increases  the  SNR  after  match  filtering. 

This  feature  can  solve  the  challenge  of  limited  range  for  UWB,  because  higher  SNR  in  turn  compensates  the  path 
loss  costs  by  distance.  SNR  gain  equals  to  TB  product  of  a  chirp,  which  can  be  expressed  in  Equation  5.4. 

SNR=^--TB,  indB  (5.4) 

No 

As  a  result,  the  receiver  sensitivity  can  be  as  low  as  (EJ,//V0  -  TB)  in  dB.  Since  the  entire  chirp  signal  is  spread 
in  time  and  frequency,  an  increased  signal  power  can  be  achieved  by  increasing  %  or  bandwidth,  while  maintaining 
the  same  transmit  power  and  PAPR.  This  is  fundamentally  different  from  IR-UWB.  Take  a  CSS  system  for  example. 
A  Chirp  UWB  pulse  has  T  —  IfiS  and  B  =  I  GHz,  and  the  processing  gain  is  TB  =  30  in  dB.  BPSK.  modulation 
is  used.  Under  the  FCC  power  spectrum  limit  and  -1  !4dBm  noise  level,  the  transmit  power  is  -4  l.3dBm  and  the 
transmit  SNR  is  72.7dB.  The  path  loss  model  of  LOS  channel  [1 J  is  depicted  in  Figure  5.4.  With  30dB  processing 
gain,  the  receiver  sensitivity  is  about  -20dB,  which  means  that  the  communication  range  can  be  1000m  with  BER  of 
10-3.  In  this  case,  when  BPSK  modulation  is  used,  the  data  rate  is  I  Mbps. 


5.1.3  NBI  Suppression 

Chirp  based  system  also  has  a  simple  solution  for  NBI  suppression.  Chirp  signal’s  spectrum  can  be  controlled  in  the 
lime  domain.  If  NBI  appears  at  a  certain  frequency,  that  frequency  can  be  filtered  out  by  ’gating’  in  the  time  domain. 
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Time 


Figure  5.1:  Waveform  of  an  up  chirp 

The  block  diagram  is  dipictcd  in  Figure  5.5.  Figure  5.6  is  a  ’gated’  chirp.  Figure  5.7  is  its  amplitude  spectrum. 
Equation  5.5  and  Figure  5.5  illustrate  the  ’gating’  process. 


M*)  = 

t  = 

i=l  ' "  ' 

It  can  be  noticed  that  a  “gated”  chirp  is  summed  up  with  N  continuous  chirp  signals.  Each  rr*  represents  a  frequency 
component.  A  certain  frequency  can  be  “filtered”  by  setting  corresponding  it*  to  0.  If  m,  —  ]  for  all  i,  the  summed 
signal  becomes  a  normal  chirp. 

From  the  amplitude  response  we  can  conclude  that  such  ’gated’  chirp  acts  like  a  notch  filter.  However,  this  notch 
filter  is  easy  to  implement.  It  has  been  reported  that  this  gating  process  can  be  done  by  the  surface  acoustic  wave 
(SAW)  device,  which  is  analog  [2],  Figure  5.8  and  Figure  5.9  shows  the  NBI  suppression  result.  It  can  be  seen 
that  when  SIR  is  -20dB,  compared  with  the  normal  match  filter,  the  NBI  is  suppressed  successfully  after  the  ’gated' 
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Frequency 


Figure  5.2:  Amplitude  spectrum  of  chirp 

SAW.  It  should  be  noted  that  this  gating  process  is  implemented  by  a  fast  adaptive  switching  network,  which  is 
digital.  Such  analog-digital  hybrid  system  makes  it  easy  to  control  the  spcctmm. 


5.1.4  Other  Advantages 

In  addition  to  the  range  extension  and  NBI  suppression  ability.  Chirp  UWB  has  the  advantage  of  fast  synchronization 
because  of  its  analog  match  filter  and  high  processing  gain.  Match  filter  can  be  implemented  by  the  SAW  device  in 
the  analog  domain.  SAW  is  a  passive  device.  A  compressed  pulse  with  processing  gain  will  automatically  appear 
after  match  filtering.  Timing  acquisition  becomes  a  task  of  peak  searching,  which  is  not  difficult  to  implement. 
Also,  chirp  signal  is  ’immune’  to  frequency  offset.  A  frequency  offset  Af  will  result  a  timeshift  St  in  the  match 
filter  output.  Suppose  a  chirp  signal  with  frequency  offset  A  f  has  the  form  of  Equation  5.6. 

s  (()  =  cos  ^2tt  ((/„  +  Af)t  +  *0)  (5.6) 


Equation  5.7  is  the  output  of  the  match  filter. 


5.1.  CHIRP  UWB 


75 


> 

CL) 

1  JT 

a 

a 

< 

> 

- 1 

( 

B 

1  . *  , 

r 

f  T 

-Tc  0  Tc 

Time 

Figure  5.3:  Waveform  of  a  compressed  pulse 
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s  (f)  h  (t  —  f)  dr 


rr~  fn  f  r  ,  a  rw .  sin  ()r  (A/ + /if)  (/e  |f|))  ^  ^ 

=  V4Mcos(2jt(/o  +  A/)i) - 2tt  (A/  +  pt) - 


(5.7) 


Equation  5.7  shows  that  the  output  is  still  a  compressed  pulse.  However,  comparing  it  with  Equation  5.3,  there  is  a 
timeshift  determined  by  A f.  The  relationship  between  Af  and  St  can  be  shown  in  Equation  5.8. 


(5.8) 


Compared  with  the  sweep  rate  mu,  frequency  offset  is  very  small,  and  the  influence  can  be  ignored.  For  example, 
if  we  add  Af  =  1  MHz  to  a  chirp  with  p  =  500  MHz//ts,  the  timeshift  is  only  2  ns.  In  all,  the  tasks  of  timing  and 
frequency  synchronization  can  be  accomplished  together  by  simple  peak  searching. 

Since  Chirp  UWB  satisfies  the  definition  of  UWB,  it  has  all  the  merits  of  traditional  UWB.  It  has  the  same  multipath 
resolution  as  IR-UWB.  Ranging  is  not  a  problem  for  Chirp  UWB,  since  pulse  compression  with  CSS  has  already 


SNR  at  Receiver  (dB) 
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Figure  5.4:  Path  loss  model  for  indoor  LOS  channel 


Received 

Signal 


Figure  5.5:  Block  diagram  for  gating  process 


been  popular  in  RADAR  applications.  Besides,  Chirp  UWB  can  be  implemented  by  the  analog  SAW  device,  its 
power  consumption  will  be  low  and  it  will  have  long  battery  life.  As  a  result.  Chirp  UWB  provides  the  ability  to 
enhance  current  UWB  performance. 


5.2.  ADVANTAGES  OF  CSS  OVER  DSSS 
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Figure  5.6:  Waveform  of  a  ’gated’  chirp 

5.2  Advantages  of  CSS  over  DSSS 

CSS  is  similar  to  Direct  Sequence  Spread  Spectrum  (DSSS).  Instead  of  spreading  the  spectrum  using  PN  code 
like  DSSS,  CSS  uses  chirp  signals.  However,  CSS  is  featured  with  Doppler  immunity,  fast  synchronization  and 
controllable  spectrum,  which  are  impossible  for  DSSS.  In  RADAR  applications,  DSSS  is  also  a  choice  for  pulse 
compression,  But  it  is  well  known  that  BAKER  code  is  very  sensitive  to  Doppler.  DSSS  with  large  bandwidth 
and  high  processing  gain  requires  very  high  speed  Digital -to- Analog  Converter  (DAC)  and  Digital  Signal  Processor 
(DSP).  Both  DAC  and  DSP’s  performances  are  limited  by  the  state-of-art  and  they  consume  a  lot  of  dissipation. 
For  synchronization,  DSSS  needs  to  do  the  digital  domain  correlation  to  generate  a  compressed  pulse  and  detect  its 
autocorrelation  peak.  This  process  cost  a  lot  of  hardware  resource  and  increases  the  system  complexity.  It  is  also 
impossible  for  DSSS  to  filter  out  NB1  as  easily  as  CSS,  because  the  pulse  for  DSSS  is  fixed. 


5.3  Implementation 

The  key  device  in  the  Chirp  UWB  system  is  SAW.  SAW  is  widely  used  as  bandpass  filter,  delay  line  and  pulse 
compression  Dispersive  Delay  Line  (DDL).  It  is  estimated  that  the  current  world-wide  production  of  SAW  is  well  in 
excess  of  1000  million  devices  per  year  [3].  For  CSS  purpose,  DDL  SAW  is  exclusively  used.  A  DDL  based  solution 
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Frequency  (Hz) 


Figure  5.7:  Amplitude  spectrum  of  the  ’gated’  chirp 

was  approved  by  IEEE  802. 1 5.4a  WPAN  task  group  in  March,  2007.  In  the  system,  a  pair  of  SAW  devices  arc  used 
as  chirp  generator  and  match  filter.  Each  SAW  generates  an  up/down  chirp  with  80MHz  bandwidth  and  l/ts  duration 
[4],  It  is  reported  that  under  1CT3  Bit  Error  Rate  (BER),  such  communication  link  can  support  transmission  distance 
up  to  9800m  with  I W  transmit  power  [5].  However,  the  acceptance  of  IEEE  standart  only  justifies  CSS's  perfor¬ 
mance  for  wideband.  CSS’s  UWB  application.  Chirp  UWB,  should  be  of  interest.  The  state-of-art  performances  of 
SAW  devices  can  generate  a  chirp  signal  with  bandwidth  over  500MHz  [6,  7,  8]  or  fraction  bandwidth  over  25%  [9], 
satisfying  both  UWB  definitions.  In  addition,  small  size,  low  cost  and  mass  production  are  the  features  of  SAW  [6], 
which  will  be  perfect  for  sensor  network  applications.  Combining  CSS  with  UWB,  Chirp  UWB  system  will  inherit 
the  advantages  of  both  CSS  and  UWB. 


5.4  The  History  of  CSS 

The  idea  of  pulse  compression  was  proposed  by  different  researchers  in  different  countries  after  World  War  II 
[10,  II,  12,  13].  In  the  late  fifties,  several  major  labs  published  their  research  results  on  chirp  RADAR  [14,  15,  16], 
The  first  paper  applying  chirp  for  communications  other  than  RADAR  was  published  in  1962  by  M.  Winkler  [17], 
where  he  suggested  using  up  and  down  chirps  to  represent  binary  ’0’  and  *1%  respectively.  This  is  the  basics 
of  a  Binary  Orthogonal  Keyed  (BOK)  system,  utilizing  the  quasi-orthogonal  property  of  up  and  down  chirps.  In 
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Figure  5.8:  Normal  match  filter  output  with  NBI,  SIR  =  -20dB 


1966,  Hala  [18]  and  Gott  [19]  independently  proposed  a  BOK  system  for  the  High  Frequency  (HF)  band.  In  their 
proposals,  chirp  signal’s  Doppler  Effect  immunity  was  utilized.  After  Gott’s  work,  Forgan  [20]  and  Edge  [21]  had 
finished  their  research  on  chirp  signals  and  they  all  reference  some  Royal  Aircraft  Establishment  and  Ministry  of 
Defense  reports.  It  seems  that  their  research  was  supported  by  the  British  military  and  they  had  shown  great  interests 
in  CSS.  The  paper  of  Zaytsev  [22]  also  showed  some  interests  from  the  Russian  military.  In  the  early  seventies,  Gott 
and  his  group  published  several  papers  on  Low  Earth  Orbit  (LEO)  satellites  and  ground-to-air  communications  in 
the  HF  band  [23,  24,  25,  26,  27].  However,  the  idea  of  using  SAW  device  as  a  chirp  generator  was  first  proposed 
by  Bush  [28]  in  1973.  Before  that,  the  only  practical  chirp  signal  generator  was  lumped  element  filters.  Interest  in 
CSS  became  less  since  late  seventies.  At  that  time  DSSS  and  Frequency  Hopping  Spread  Spectrum  (FHSS)  was  the 
more  popular  choice  because  of  the  development  of  Digital  Signal  Process  (DSP)  algorithms.  Only  a  few  papers 
mentioned  CSS  by  applying  it  to  DSSS  against  the  Doppler  Effect  [29, 30, 3 1  ].  At  that  time,  SAW  device  was  mainly 
developed  for  RADAR  pulse  compression  [32,  33,  34],  satellite  signal  processing  [35,  36, 37]  and  channel  sounders 
[38,  39].  Interest  in  CSS  fell  can  be  explained  by  the  target  frequencies  and  the  SAW  implementation  techniques  of 
that  time.  The  excess  delay  for  HF  band  is  in  the  order  of  milliseconds,  which  means  that  the  chirp  signal’s  duration 
must  be  no  less  than  that  level.  Since  the  velocities  in  SAW  arc  thousands  of  meters  per  second,  the  corresponding 
physical  size  of  SAW  would  be  over  1  meter,  which  is  impractical. 

However,  when  the  frequencies  are  not  limited  in  HF  band,  SAW  devices  found  its  applications.  For  indoor  appli¬ 
cations,  the  excess  delay  is  in  the  order  of  microseconds,  which  is  practical  for  physical  size  of  SAW.  SAW  devices 


so 
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were  first  applied  in  the  indoor  channels  by  Tsai  and  Chang  in  1994,  using  BOK  scheme  [40],  After  that,  several 
groups  took  interests  in  CSS.  Hucmer  [41]  in  1998  published  another  BOK  based  system.  Gugler's  group  [42]  and 
Pinkney  [43]  focused  their  CSS  research  on  Direct  Modulation  (DM),  another  CSS  modulation  scheme  which  is 
more  efficient  than  BOK.  In  their  study,  the  data  rate  could  be  as  high  as  70Mbps  with  80MHz  bandwidth.  In  March, 
2007,  IEEE  approved  Nanotron’s  CSS  physical  layer  (PHY)  in  its  new  wireless  standard  802. 15.4a.  A  pair  of  SAWs 
are  used  in  its  BOK  scheme.  The  target  applications  for  such  system  are  Real  Time  Location  Systems  (RTLS),  in¬ 
dustrial  control,  sensor  networking  and  medical  devices  in  the  ISM  band,  taking  the  advantages  of  CSS:  long  range, 
high  robustness  and  low  power  consumption. 


5.5  A  Time  Reversed  Chirp  UWB — Forward  Looking 


Chnp  UWB  radio  can  be  regarded  as  a  compromise  between  the  optimum  waveform  and  ease  of  implementation — 
analog  or  digital,  by  selecting  chirp  waveform  for  signal  modulation.  The  research  proposed  here  has  a  significance 
of  practical  value,  and  the  chirp  system  is  used  as  a  benchmark  to  be  compared  against  the  optimum  modal  modula¬ 
tion. 

A  RF  harsh  environment  shown  in  Fig.  2.1(a)  has  been  investigated.  One  simulation  example  is  included  here.  A 


5.5.  A  TIME  REVERSED  CHIRP  UWB— FORWARD  LOOKING 


81 


Figure  5.10:  The  photo  of  a  pair  of  SAW  filters.  The  size  is  about  1  x  2  inches. 

data  rate  of  500Mbps  can  be  supported  using  time  reversal  chirp  waveform  with  a  bandwidth  7GHz  and  waveform 
duration  T  =  100  ns — the  time-bandwidth  product  is  TB  =  700.  BPSK  modulation  is  used.  For  a  bit  error  rate 
{BER)  of  10-4,  the  ratio  of  bit  energy  to  the  noise,  Et/No,  is  required  to  be  larger  than  23  dB. 

One  PhD  student  (Peng  Zhang) — who  did  his  MS  thesis  on  chirp  UWB  for  three  years — joined  the  Pi’s  team  in  2006. 
He  is  currently  working  on  the  testbed  for  chirp  UWB,  using  a  SAW  chirp  device.  On  this  topic,  Mr.  Zhang  is  work¬ 
ing  with  one  R&D  engineer  Dr.  N.  Guo  and  another  MS  student  Amanpreet  Singh  Saini.  A  pair  of  SAW  filters  shown 
in  Fig.  5.10  arc  available  at  the  Pi’s  lab.  The  vendor  is  Sawcom  Tech,  Inc.  (http  :  // www.8awr.omtech.cotn/). 

This  is  the  first  phase  of  the  chirp  UWB  radio.  The  goal  is  to  prove  the  concept  as  quickly  as  possible.  The  second 
phase  of  this  project  is  to  cope  with  the  cognitive  radio  concept,  which  is  the  proposed  topic  for  this  proposal.  The 
SAW  filters  can  be  electronically-controlled.  Narrowband  interference  suppression  is  an  important  research  task, 
due  to  the  co-existence  of  UWB  and  narrowband  systems  such  as  Wi-Fi  and  Wi-MAX. 

A  chirp  UWB  system  combines  the  concepts  of  CSS  and  UWB  with  the  following  features: 

1 .  Long  range  transmission 

2.  NBI  suppression 

3.  Fast  synchronization 

4.  Low  transmit  power,  low  probability  of  intercept 

5.  Low  power  consumption 

6.  Low  system  complexity,  low  cost 

In  summary,  chirp  UWB  is  practical  and  useful.  Its  combination  with  time  reversal  is  a  promising  transmission 
scheme.  Time  reversal  deals  with  the  extremely  quasi-continuous  multipath — unique  to  UWB — while  a  pair  of 
chirp  pulses  is  used  for  signaling.  Time-bandwidth-product  (TB)  solely  determines  the  system  performance.  The 
range  and  data  rate  can  be  trade-offed,  by  conveniently  controlling  the  time-bandwidlh-producl. 
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Chapter  6 


UWB  Channel  Sounding 


In  this  chapter,  UWB  channel  sounding  is  discussed.  In  general,  our  Lab  has  the  capabilities  of  sounding  the  channel 
in  both  time  domain  and  frequency  domain,  with  and  without  multiple  antennas.  The  goal  of  these  measurements 
is  to  evaluate  the  non-fading  wireless  communications  enabled  by  time  reversal  transmission.  Both  space  and  time 
are  investigated  as  a  function  of  signal  bandwidth.  A  critical  question  is:  What  is  the  optimum  signal  bandwidth? 
This  question  depends  on  the  specific  transmission  scheme — we  are  interested  in  two  schemes:  time  reversal  modal 
modulation,  and  time  reversal  chirp  UWB. 


6.1  Measurement  Setup 

The  knowledge  about  the  UWB  radio  channel  can  be  gained  by  using  either  time  or  frequency  domain  measurement 
technique.  We  will  discuss  both  techniques  in  detail  in  this  section. 


6,1.1  Time  Domain  Measurement  Technique 

We  will  discuss  the  time  domain  measurement  technique  based  on  impulse  transmission.  In  this  technique  a  narrow 
pulse  is  sent  through  the  propagation  channel  and  the  channel  impulse  response  is  recorded  by  using  a  Tektronix 
Digital  sampling  oscilloscope  (DSO),  The  pulse  sent  through  the  channel  is  modulated  by  a  special  waveform  like 
sine,  square  or  ramp  which  enables  us  to  analyze  the  effects  of  different  paths  on  the  received  signal[l  ].  The  data 
stored  in  DSO  can  be  processed  by  using  a  computer.  The  bandwidth  of  received  signal  depends  on  the  shape  and 
width  of  transmitted  pulse.  The  time  domain  measurement  setup  consists  of  a  pulse  generator,  transmitter  antenna 
and  receiver  antenna,  a  triggering  signal  generator.  Low  noise  amplifier  (LNA)  and  the  DSO.  The  setup  for  time 
domain  measurement  is  shown  in  Fig.  6. 1 

The  whole  setup  consist  of  two  section,  he  transmit  and  receive  part.  The  signal  generator  and  pulse  generator 
constitute  the  transmitter  part  and  DSO  along  with  LNA  constitutes  the  receiver  part.  The  signal  generator  is  used 
to  trigger  the  pulse  generator  and  pulse  generator  generates  the  pulse  that  is  transmitted  through  the  channel.  On 
the  receiver  side  the  signal  is  passed  through  LNA  and  amplified.  The  final  results  are  displayed  on  the  DSO. 
A  triggering  signal  from  signal  generator  is  used  to  synchronize  the  DSO  to  record  the  measurements.  By  using 
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Triggo-  &gnal 


Figure  6. 1 :  Time  domain  measurement  setup  block  diagram 
averaging  the  signal  to  noise  ratio  (SNR)  is  improved. 

The  main  advantages  of  this  technique  are  less  complexity,  lower  cost  and  channel  responses  is  readily  available  in 
time  domain.  With  multiple  antennas,  our  system  can  sound  the  UWB  MIMO  channel  in  the  time  domain,  by  using 
two  sample  modules  from  Tektronix. 

6.1.2  Frequency  Domain  Measurement  Technique 

In  the  frequency  domain  measurement  technique,  the  channel  sounding  is  carried  out  by  sweeping  a  set  of  narrow- 
band  sinusoid  signals  (tones)  through  a  wide  frequency  band.  Each  tone  has  the  same  power  (typically  3dBm  is  used 
in  our  experiments).  The  channel  frequency  response  is  recorded  by  using  a  vector  network  analyzer  (VNA).  The 
VNA  is  operated  in  transfer  function  mode  where  one  of  its  port  serves  as  transmitting  port  and  the  other  as  receiving 
port.  S-parameters  are  used  to  measure  the  complex  frequency  domain  transfer  function.  The  setup  for  frequency 
domain  measurements  is  shown  in  Fig.  6.2 


Figure  6,2:  Frequency  domain  measurement  setup  block  diagram 

Generally  S21  parameter  is  used  where  port  1  and  2  are  acting  as  transmitting  and  receiving  port  respectively. 
VNA  sends  a  frequency  tone  /  through  the  channel  and  the  frequency  response  of  the  channel  is  represented  as 
$21  (/)  corresponding  to  frequency  tone  f.  By  sweeping  a  frequency  range  from  to  /i  we  can  get  the  frequency 
response  of  that  particular  band.  If  N  is  the  number  of  frequency  points  per  sweep  with  frequency  step  fc  MHz,  their 
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relationship  with  bandwidth  B  MHz  can  be  represented  as 


where  B  is 


k  = 


B 

N  -  I 


B  =  A  -  A 


(6-1) 

(6.2) 


The  maximum  detectable  delay  of  the  channel  can  be  calculated  as  [2] 

_  N  -  1  ,,,, 

Tmax  “  g  (b-3) 

The  main  advantage  of  the  frequency  domain  measurement  technique  is  that  the  cables  and  connectors  can  he 
calibrated  before  measurement  to  compensate  for  different  kind  of  errors  and  frequency  dependent  losses  that  can 
occur  during  the  measurement.  There  are  three  different  kind  of  measurement  errors  [3] 


•  Systematic  errors 

•  Random  errors 

•  Drift  errors 


Systematic  errors  are  caused  by  imperfections  in  the  test  equipment  and  test  setup.  These  errors  can  be  removed  by 
calibration.  Systematic  errors  are  related  to  signal  leakage,  signal  reflections,  and  frequency  response. 

Random  errors  are  random  in  nature  and  are  unpredictable.  As  a  result  they  cannot  be  removed  by  calibrations. 
Random  errors  are  present  mainly  due  to  instrument  noise.  So  they  can  be  removed  by  increasing  the  source  power 
or  by  narrowing  the  IF  bandwidth. 

Drift  errors  occurs  when  the  system  performance  changes  afteT  the  completion  of  calibration  process.  Drift  er¬ 
rors  occurs  mainly  due  to  change  in  temperature  of  measurement  environment  and  can  be  removed  by  additional 
calibrations. 

A  mechanical  calibration  kit  consists  of  a  set  of  physical  devices  called  standards  [4].  Each  standard  has  a  precisely 
known  magnitude  and  phase  response  as  a  function  of  frequency.  These  standards  are 

•  Gpen,short/offset  short, load/sliding  load 

•  Adapters 

•  Precision  offsets  waveguide  or  coaxial 


There  arc  various  calibration  method.The  four  most  widely  used  are 


#  short,  open,  load,  and  through  (SOLT) 

•  LRM 
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•  LRRM 

•  TRL 


Out  of  the  above  mentioned  methods  the  most  dominant  one  is  SOLT.  With  SOLT,  the  four  known  references 
are:  short,  open,  load,  and  through.  Once  the  calibration  is  done,  the  calibration  state  can  be  saved  and  recall 
while  performing  the  experiment.  By  using  frequency  domain  measurement  technique  we  can  get  larger  dynamic 
range  which  improves  the  measurement  precision.  The  channel  impulse  response  (CIR)  which  yields  the  required 
information  to  characterize  the  UWB  channel,  is  obtained  by  taking  the  inverse  fast  Fourier  transform  (IFFT)  of  the 
recorded  signal.  So  this  technique  requires  extra  signal  processing. 

One  of  the  main  limitation  of  frequency  domain  measurement  technique  is  that  it  need  static  environment  through  out 
the  measurement.  But  Time  domain  measurement  technique  can  support  non-stationary  channels.  Another  problem 
is  the  high  attenuation  of  the  cables  so  frequency-domain  measurement  is  good  for  short  range. 


6.2  Measurement  Examples 

In  this  section,  we  will  discuss  some  of  the  experiments  to  gain  more  insight  about  the  channel  sounding. 


6.2.1  Measurement  for  SISO  Case 

The  measurement  is  performed  is  for  SISO  case  in  three  different  environments  i.e  rectangular  metal  cavity,  hallway 
and  office  to  analyze  the  characteristics  of  the  UWB  channel  in  these  environments.  The  VNA  used  is  Agilent 
N5230A  (300kTb.-I3.5GHz),  The  calibration  method  used  is  SOLT.  Table  6.1  list  the  main  parameters  for  the 
measurement  This  set  of  parameters  is  used  for  measurement  in  all  the  three  cases.  The  measurement  is  done  for 


Table  6. 1 :  Measurement  Parameters  setup 


Parameter 

value 

Frequency  Band 

3GHz-10GHz 

Bandwidth 

7GHz 

Number  of  points 

7001 

Transmission  power 

lOdBm 

Frequency  step 

1MHz 

Antenna  Polarization 

vertical 

Averaging  number 

128 

Line  of  sight  (LOS)  situation.  The  distance  between  transmitter  antenna  and  receiver  antenna  varies  from  0.5m  to 
4m  in  steps  of  0.5m  and  both  the  antennas  are  placed  at  L35  meters  above  the  ground  for  all  the  three  measurement 
environments. 
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Measurement  in  Rectangular  Metal  Cavity 

The  rectangular  metal  cavity  used  in  this  measurement  measures  1 6  feet  by  8  feet  by  8  feet.  The  material  of  inside 
walls  of  the  cavity  is  aluminum.  The  setup  for  measurement  in  rectangular  metal  cavity  is  shown  in  Fig.  6.3 


Figure  6.3:  Setup  for  measurement  in  rectangular  metal  cavity 


The  transmitter  antenna  is  fixed,  and  the  receiver  antenna  is  moved  from  distance  0.5m  to  4m  in  steps  of  0.5m  w.r.t 
transmitter  antenna  along  the  middle  line  of  the  rectangular  metal  cavity.  The  frequency  domain  response  at  distance 
Im,  2m,  3m  and  4m  between  transmitter  antenna  and  receiver  antenna  is  shown  in  Fig.  6.4 
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Figure  6.4:  Frequency  domain  response  H(f)  at  distance  lm,2m,3m  and  4m  between  transmitter  antenna  and 
receiver  antenna  in  rectangular  metal  cavity 


After  the  channel  response  is  measured,  zero  padding  is  done  between  QHz  and  3  GHz,  conjugate  reflection  between 
-10  GHz  and  QHz  and  inverse  fast  Fourier  transform  (IFFT)  is  used  to  get  channel  impulse  response  in  time  domain. 
The  time  domain  response  at  distance  1m,  2m,  3m  and  4m  between  transmitter  antenna  and  receiver  antenna  is 
shown  in  Fig.  6.5 
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Figure  6.5:  Channel  impulse  response  h(t )  at  distance  Im,  2m,  3m  and  4m  between  transmitter  antenna  and  receiver 
antenna  in  rectangular  metal  cavity 

The  channel  impulse  response  forms  a  channel  autocorrelation  function  according  to 

Rkh  =  MO  *  (6.4) 

where  h’(-t)  is  the  time  reversed  conjugated  version  of  h(t).  The  autocorrelation  of  channel  pulse  response  at 
distance  I  m,  2m,  3m  and  4m  between  transmitter  antenna  and  receiver  antenna  is  shown  in  Fig.  6.6 
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Figure  6.6:  Autocorrelation  function  at  distance  1m,  2m,  3m  and  4m  between  transmitter  antenna  and 

receiver  antenna  in  rectangular  metal  cavity 

Measurement  in  Office  Environment 

The  same  procedure  of  measurement  is  repeated  in  office  environment.  The  same  set  of  parameters  as  shown  in 
table  6.1  is  used.  The  frequency  domain  response  at  distance  1m,  2  m,  3m  and  4m  between  transmitter  antenna  and 
receiver  antenna  is  shown  in  Fig.  6.7 
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Figure  6.7:  Frequency  domain  response  H(f)  at  distance  Im,2m,3m  and  4m  between  transmitter  antenna  and 
receiver  antenna  in  office 


The  time  domain  response  at  distance  Im,  2m,  3m  and  4m  between  transmitter  antenna  and  receiver  antenna  is 
shown  in  Fig.  6.8 
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Figure  6,8:  Channel  impulse  response  h{t)  at  distance  lm,2m,3m  and  4m  between  transmitter  antenna  and  receiver 
antenna  in  office 

The  autocorrelation  of  channel  pulse  response  at  distance  1m,  2m,  3m  and  4m  between  transmitter  antenna  and 
receiver  antenna  is  shown  in  Fig.  6,9, 
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Figure  6.9:  Autocorrelation  function  Hhh(t)  at  distance  lm,  2m,  3m  and  4m  between  transmitter  antenna  and 
receiver  antenna  in  office 

Measurement  in  Hallway  Environment 

The  same  procedure  of  measurement  is  repeated  in  hallway  environment.  The  same  set  of  parameters  as  shown  in 
table  6.1  is  used.  The  frequency  domain  response  at  distance  lm,  2m,  3m  and  4m  between  transmitter  antenna  and 
receiver  antenna  is  shown  in  Fig.  6.10 
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Figure  6.10:  Frequency  domain  response  H(f)  ai  distance  lm>2in,3m  and  4m  between  transmitter  antenna  and 
receiver  antenna  in  hallway 


The  lime  domain  response  al  distance  lms  2m,  3m  and  4m  between  transmitter  antenna  and  receiver  antenna  is 
shown  in  Fig.  6.!  I 
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Figure  6. !  I :  Channel  impulse  response  h(t)  at  distance  Im,2m,3m  and  4m  between  transmitter  antenna  and  receiver 
antenna  in  hallway 


The  autocorrelation  of  channel  pulse  response  at  distance  Im,  2m,  3m  and  4m  between  transmitter  antenna  and 
receiver  antenna  is  shown  in  Fig,  6. 1 2 
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Figure  6,12:  Autocorrelation  function  Hkh(t)  at  distance  Im,2m,3ni  and  4m  between  transmitter  antenna  and  re¬ 
ceiver  antenna  in  hallway 

Comparison  of  Energy  for  Rectangular  Metal  Cavity,  Office  and  Hallway  Environments 


The  energy  of  channel  impulse  response  for  the  three  cases  can  be  calculated  by  using 

Eh  =  fTh  h\t)dt 
J  o 


(6.5) 


where  'J),  is  ihc  time  duration  of  the  channel.  The  energy  for  the  three  cases  is  shown  in  Fig.  6.13 

From  the  Fig.  6.13,  it  can  be  seen  clearly  that  the  total  energy  of  the  channel  impulse  response  in  a  confined  metal 
cavity  decays  very  slowly,  as  compared  to  those  in  office  and  hallway. 
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Figure  6. 1 3:  Energy  of  channel  impulse  response  for  rectangular  metal  cavity,office  and  hallway  environments 


6.2.2  Measurement  for  MIMO  Case 

The  measurement  is  performed  for  a  3x3  MIMO  case  using  virtual  array  inside  the  rectangular  metal  cavity.  The 
measurement  parameters  used  are  the  same  as  described  in  table  6. 1 .  The  distance  between  the  antennas  in  transmit¬ 
ter  array  and  receiver  array  is  2  feet.  The  setup  for  the  measurement  is  shown  in  Fig.  6.14 
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Figure  6. 14:  Setup  for  3x3  virtual  MIMO  case  in  rectangular  metal  cavity 
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Figure  6.15:  Frequency  domain  response  //(/)  from  tansmittcr  antenna  I  to  different  receiver  antennas 


Figure  6. 16:  Channel  impulse  response  h(t)  from  transmitter  antenna  I  to  different  receiver  antennas 


The  frequency  domain  response  and  time  domain  response  from  transmitter  antenna  1  to  different  receiver  antennas 
is  shown  in  Fig.  6.15  and  Fig.  6.16 
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Figure  6. 1 7:  Frequency  domain  response  //(/)  from  transmitter  antenna  2  to  different  receiver  antennas 


Figure  6*  18:  Channel  impulse  response  h(t)  from  transmitter  antenna  2  to  different  receiver  antennas 


The  frequency  domain  response  and  time  domain  response  from  transmitter  antenna  2  to  different  receiver  antennas 
is  shown  in  Fig*  6J7  and  Fig*  6.18 

The  frequency  domain  response  and  lime  domain  response  from  transmitter  antenna  3  to  different  receiver  antennas 
is  shown  in  Fig.  6.19  and  Fig*  6.20 
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Figure  6,20:  Channel  impulse  response  h(t.)  from  transmitter  antenna  3  to  different  receiver  antennas 


63  Conclusion 


In  this  chapter,  channel  sounding  for  UWB  channel  is  discussed.  Time  domain  measurement  technique  and  fre¬ 
quency  domain  measurement  technique  are  discussed  in  detail.  Then  two  measurement  examples  are  considered. 
Firstly  we  compare  the  characteristics  of  UWB  SISO  channel  in  rectangular  metal  cavity,  in  office  and  in  hallway, 
and  then  the  measurement  case  for  a  3x3  UWB  M1MO  channel  in  rectangular  metal  cavity  is  also  considered. 
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Chapter  7 


Overview  of  Testbed 


This  report  documents  some  of  our  recent  development  work  with  focus  on  design  and  implementation  aspect. 

The  ultimate  goal  is  to  achieve  range  extension  by  optimizing  system  design.  Specifically,  the  overall  goal  is  to 
implement  a  UWB  radio  test-bed  with  waveform  optimization.  The  work  of  design  is  to  provide  a  feasible  solution 
under  some  theoretical  and  practical  criteria.  We  proposed  a  system  architecture  relying  on  eontinuous-time  (or 
waveform  level)  lime  reversal  technique,  though  there  are  a  few  other  types  of  optimization  rules.  One  of  the  key 
components  in  the  UWB  time  reversal  system  is  the  programmable  waveform  generator  based  on  advanced  FPGA 
and  D/A  conversion  technologies,  A  XilJinx  Virtex  5  chip,  one  of  the  fastest  FPGA  chips  to  date,  is  chosen  to 
perform  major  digital  back-end  functions  at  the  transmitter.  The  ultra- wide  bandwidth  and  ultra-high  sampling  rate 
pose  many  design  and  implementation  challenges.  Major  difficulties  are  associated  with  the  interfaces  between  the 
FPGA  and  mixed-signal  devices,  timing  closure  issues  in  FPGA  routing,  fixed-point  implementation  with  limited 
hardware  resources,  etc. 

Management  wise,  in  additional  to  frequently  documenting  and  periodically  reporting,  we  have  come  up  with  a 
manageable  mechanism  including  detailed  plans  and  goals  for  different  periods  of  times,  weekly  project  meetings, 
monthly  project  reviews,  and  work  assignments  on  daily  basis. 

Currently  the  system  architecture  has  been  stabilized  and  we  arc  working  toward  the  next  milestone  (see  milestone 
list  in  Table  l). 

The  rest  part  of  this  report  is  organized  as  follows. 


_ Table  1,  Milestones _ _ _ 

Stable  Architecture  8/2007 

Time  Reversal  System  with  Cable  1/2008 

Time  Reversal  System  over  the  Air  (Some  Field  Test)  6/2008 
Fully  Functional  UWB  System  with  Time  reversal  1/2009 
System  Field  Test  6/2009 

Potential  Technology  Transfer  to  US  Navy_ 1 2/2009 
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Chapter  8 


System  Design  and  Architecture 


With  the  overall  goal  and  constraints  from  both  theoretical  and  practical  perspectives,  a  completed  top-down  solution 
has  been  given.  A  very  broad  range  of  issues  need  to  be  considered.  Here  just  name  a  few:  frequency  and  band,  data 
rate,  modulation,  synchronization,  coexistence,  commercial  availability  of  devices,  etc.  Simply  speaking,  a  feasible 
and  advanced  design  must  be  made. 


8.1  Bandpass  Time  Reversal 

From  implementation  point  of  view,  it  is  too  difficult  to  build  a  system  working  from  DC  to  500  MHz  or  higher. 
Instead,  a  pass-band  system  is  more  feasible  and  flexible.  To  reduce  the  requirement  for  sampling  rate,  a  typical 
technique  used  with  digital  signal  processing  is  up-  and  down-  conversion.  The  carrier  phase  difference  between  the 
received  signal  and  the  output  of  local  oscillator  is  usually  unknown.  We  will  show  that  quadrature  conversion  can 
take  care  of  this  unknow  phase  difference. 

For  convenient  description  in  the  following,  complex  representation  is  accepted  to  represent  bandpass  signals.  As¬ 
sume  a  bandpass  sounding  signal  is  transmitted  over  the  channel,  where  uj  is  the  carrier  frequency  and 

#ti  is  an  unknown  initial  phase.  All  filtering  effect  introduced  over  the  radio  link  is  taken  into  account  by  p(t).  Con¬ 
sider  a  channel  with  its  channel  impulse  response  (CIR)  A(fc)  modeled  as  tap  delay  line  (implying  no  consideration 
of  frequency  selective  pulse  distortion) 


(8!) 


t 


Up  on  receiving  the  noisy  sounding  signal 


(8-2) 


the  receiver’s  down-conversion  stage  outputs  a  baseband  signal 


(83) 
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where  0\  is  an  unknown  phase  introduced  in  down-converting.  There  are  a  variety  of  estimation  techniques  to 
reduce  the  impact  of  noise  rt'(t)  and  multiple  sounding  pulse  may  be  necessary.  To  describe  the  central  idea  more 
conveniently,  we  assume  perfect  estimation  here.  Based  on  this  baseband  signal  y(t),  there  can  be  a  few  waveform 
for  performing  time  reversal  at  the  transmitter.  Consider  the  following  two  optional  waveforms: 
option  1 

h(t)  =  ^ai<S(f-Tj). 
i 

option  2 

p(t)  =  ^  ctip(t  —  Ti)e>(WTi_flo+91*.  (8.4) 

The  option  I  needs  a  great  deal  of  effort  in  the  estimation  stage,  while  the  option  2  is  extremely  simple  since  it 
does  not  require  any  parametric  estimation.  In  particular,  the  option  2  can  work  directly  with  “arbitrary”  waveform 
generator  with  need  for  pulse  shaping. 

Let  us  consider  a  simple  case  of  one  pulse  per  symbol  and  use 

p(7i,  - =53a*P(^b  -  *-'n)ei(M*+T,,_9o+S,+fal  (8.5) 

i 

as  the  symbol  waveform,  where  Tq  is  a  constant  to  guarantee  the  transmitted  signal  being  causal.  The  received 
symbol  waveform  can  be  expressed  as 

ia*p(7o  -  t  -  n  + 

i  k 

=  Y, afr(T<>  -  t)e>^t-9o+9i+e^)  +  Y  <*i<*kP(Ta  -t-Ti  +  Tk)e* M*+* -*>>-* +*i  +*1,  (8.6) 

i  i^k 

which  leads  to  a  baseband  version  after  down-conversion: 

Y  «?PC*0  -  W*  +  E  a*aM -t-Ti  +  Tk)^-^+e\  (8.7) 

i  i/fc 


where 

9  =  —Oq  +  6?i  +  $2  +  (8,8) 

has  been  applied.  In  (8,6)  and  (8,7)  the  summations^  and  'H&k  correspond  to  main  lobe  and  side  lobes,  respec¬ 
tively,  For  coherent  demodulation,  there  has  to  be  some  mechanism  to  estimate  the  unknow  phase  0,  but  this  is  not 
necessary  for  non-coherent  demodulation*  Also,  down-conversion  is  not  required  for  energy  detector  based  scheme. 


8.2  System  Architecture 

The  transmitter  and  receiver  architectures  are  shown  in  Fig.  1  and  2,  where  either  one  can  be  divided  into  four  major 
functional  parts:  R.F  front-end,  mixed-signal,  digital  back-end  and  control/intcrface  via  PC. 

8.3  Link  Budget 


The  link  budget  estimation  shown  below  is  based  on  OOK  modulation  and  free  space  propagation. 


8J.  LINK  BUDGET 
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Transmitter  Architecture 


Figure  8. 1 :  Transmitter  Architecture, 


Table  II.  A  Link  Budget  Example. 


Bit  rate 

Minimum  required  SNR  per  bit 
Center  frequency 
Transmitter  antenna  gain 
Receiver  antenna  gain  2 
Distance  between  the  two  antennas 
Receiver  noise  figure 
Implementation  loss 
Transmitted  power  to  antenna 
Bandwidth  (approximately) 

J  m  EIRP  spectral  density  over  the  whole  bandwidth 
Path  loss 

Received  power  the  antenna 
Link  margin 

Proposed  minimal  receiver  sensitivity 


5  Mb/s 

—  16 

4  GHz 
2  dBi 
dBi 

40  meters 
7dB 
6dB 

0,51  dBm  (1.13  mW) 
500  MHz 
-68,96  dBm/MHz 
76,52  dB 
-72,01  dBm 
6.00  dB 
-78.01  dBm 
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Receiver  Architecture 


Figure  8.2:  Architecture  of  a  reference  receiver  design. 


Chapter  9 


RF  Front-Ends  and  Mixed-Signal  Stage 


We  mainly  relay  on  off-the-shelf  products  in  our  test-bed  development.  Reported  in  this  section  include  selection  of 
components/modules  and  some  test  results.  Major  components  and  modules  listed  in  Table  III. 

_ Table  111.  A  list  of  Selected  Components/Modules. _ 

Local  oscillator  center  frequency  3.5  GHz  or  4  GHz 

Modulator  (up-converter)  bandwidth  500  MHz 

Dual-channel  D/A  sampling  rate  lGS/s  or  higher,  with  LVDS  ports 

A/D  sampling  rate  up  to  1 .5  GS/s 


9.1  Local  Oscillator 


PSA3550C  produced  by  Z-Communication,  Inc.  is  used  as  the  local  oscillator  (LO)  whose  frequency  range  is  from 
3545MHz  to  3555MHz.  Fig.  9.1  shows  the  picture  of  evaluation  board  ofPSA3550C. 

PSA3550C  includes  ADF4106  made  by  Analogy  Device,  which  is  used  as  the  frequency  synthesizer.  ADF4I06  is 
programmable  through  P5  port  on  the  evaluation  board.  When  we  do  the  test,  the  power  supply  is  set  to  6Vdc.  The 
amplified  RF  output  P2  is  connected  to  Spectrum  Analyzer.  P3  is  connected  to  the  termination.  No  attenuations  are 
used.  The  onboard  12MHz  reference  oscillator  is  used,  so  PI  is  not  connected  to  the  external  reference  oscillator 
temporarily.  The  REF  IN  frequency  is  set  to  12MHz.  The  desired  pfd  reference  frequency  is  set  to  lOOOKHz.  If  the 
desired  VCO  output  frequency  is  set  to  3550MHz,  the  spectrum  of  the  output  is  shown  in  Fig.  9.2. 


9.2  Modulator 


TRF3703-15  produced  by  Texas  Instruments  is  used  as  the  quadrature  modulator.  Fig.  9.3  shows  the  picture  of 
evaluation  board  ofTRF3703-I5. 

The  TRF  3703  is  a  very  low-noise  direct  quadrature  modulator,  capable  of  converting  complex  modulated  signals 
from  baseband  directly  up  to  RF.  The  TRF3703  is  ideal  for  high  performance  direct  RF  modulation  from  400 MHz 
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Figure  9.1:  The  picture  of  evaluation  board  of  PSA3550C  used  as  the  local  oscillator. 

up  to  4GHz.  The  RF  output  block  consists  of  a  differential  to  single-end  converter  and  an  RF  amplifier  capable  of 
driving  a  single-end  5011  load  without  any  need  of  external  components  [  I  ]. 

When  we  do  the  test,  we  connect  the  5Vdc  supply  to  headers  Wl  and  W2;  we  use  PSA3550C  to  supply  the  LO 
signal  at  3555MHz  to  J I  and  terminate  J2  with  50 fl  to  ground;  we  use  a  waveform  generator  to  provide  the  pulse 
waveform  whose  width  is  20ns  as  the  I  input  signal  and  set  the  common  mode  voltage  on  the  input  signals;  we 
connect  a  Spectrum  Analyzer  to  the  SMA  connector  marked  RFOUT  (J7)  and  monitor  the  output  The  spectrum  of 
the  output  is  shown  in  Fig.  9.4. 


9J  Digital  to  Analogy  Converter 


MB86064  produced  by  Fujitsu  is  used  as  the  Digital  to  analogy  converter  (DAC).  Fig.  9.5  is  the  picture  of  evaluation 
board  of  MB86064.  Fig.  9.6  shows  the  GUI  of  control  software  for  MB86064. 

The  Fujitsu  MB86064  is  a  dual  14-bit  lGSa/s  DAC,  delivering  exceptional  dynamic  performance.  Each  high 
performance  DAC  core  is  capable  of  generating  multi-stand,  multi-carrier  communication  transmitted  signals.  DAC 
data  is  input  via  two  high-speed  LVDS  ports.  These  operate  in  a  pseudo  double  data  rate  (DDR)  mode,  with  data 
latched  on  both  rising  and  falling  edges  of  the  clock.  To  simplify  system  integration  the  DAC  operates  from  a  clock 
running  at  half  the  DAC  conversion  rate  [2].  A  function  block  diagram  is  shown  in  Fig.  9.7.  The  reasons  why  we 
choose  MB86064  as  DAC  in  our  test-bed  can  be  summarized  below, 

From  clock  point  of  view,  MB86064  requires  a  RF  input  clock  at  half  the  DAC  conversion  rate,  with  sufficient 
spectral  purity  to  not  impact  the  target  analog  output  performance.  Because  the  DAC  cores  arc  clocked  on  both 
rising  and  falling  edges  of  the  input  clock,  the  requirement  of  the  external  clock  source  is  reduced.  The  sampling 
rate  of  DAC  cores  is  lGaa/s  and  the  frequency  of  the  RF  input  clock  is  500MHz  which  is  easier  to  be  generated  by 
the  current  FPGA  devices.  At  first,  we  use  AD9734  as  DAC  in  our  test-bed.  Though  the  sampling  rate  of  AD9734 
is  1.2Gsa/s,  the  DAC  clock  with  frequency  of  1.2GHz  is  needed. 
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Figure  9.2:  The  spectrum  of  the  PSA3550C  output  when  the  desired  VCO  output  frequency  is  set  3550MIiz. 


Figure  9.3:  The  picture  of  evaluation  board  ofTRF3703-l5  used  as  the  quadrature  modulator. 

From  the  number  of  DAC  cores  point  of  view,  MB86064  has  two  DAC  cores  which  can  generate  two  signals  for  1/Q 
simultaneously.  So,  we  do  not  need  to  buy  and  integrate  two  DACs  with  only  one  core  in  the  test-bed,  which  will 
greatly  reduce  the  difficulty  of  the  system  integration. 

From  the  interface  point  of  view,  The  evaluation  board  of  MB86064  is  designed  to  accommodate  up  to  six  SMA 
conncctors-two  for  the  differential  clock  input  and  four  for  the  two  differential  analog  outputs.  Boards  arc  supplied 
with  transformers  on-board  to  perform  single  ended -to-differential  and  differ cntial-to-single  ended  conversions  [3], 
So  using  MB86064,  we  can  easily  make  one  simple  adapter  to  connect  the  evaluation  board  of  MB86064  and  the 
evaluation  board  ofTRF3703-15. 

From  the  waveform  memory  module  point  of  view,  MB  86064  incorporates  a  waveform  memory  module  featuring 
two  16k  point  on-chip  waveform  memories  [4],  These  allow  the  DAC  cores  to  be  driven  with  user  programmed 
waveforms  without  the  need  for  external  high  speed,  pattern  generators  which  can  increase  the  flexibility  of  system 
design. 
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Figure  9.4:  The  spectrum  of  TRF3703-I5  output, 


Figure  9.5:  The  picture  of  evaluation  board  of  mb86064  used  as  DAC. 

Now,  we  just  test  the  function  of  the  waveform  memory  module.  According  to  the  user  manual,  we  follow  the  steps 
to  power  up  and  start  to  use  the  evaluation  board  of  MB86064.  A  single-ended  clock  source  is  connected  to  SMA 
connector  CLK;  the  single-ended  analog  outputs  on  SMAs  IOUTA  and  IOUTB  are  connected  to  the  test  equipment; 
the  jumper  links  LKI  and  LK2  are  set  as  required;  the  PC  USB  programming  cable  is  connected  to  the  CONTROL 
I/F  header;  the  power  supply  cable  is  connected  to  the  POWER  header  with  correct  orientation;  the  power  supplies 
are  turned  on;  the  PC  control  software  are  installed  and  operated. 


9.4  MAX108  A/D 


An  8-bit  monolithic  bipolar  A/D  converter  with  sampling  rate  up  to  1.5Gbps  is  employed  to  convert  the  baseband 
signal  (analog  signal)  into  digital  domain.  A  baseband  amplifier  between  the  low  pass  fitter  following  the  square  law 
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Figure  9.6:  The  GUI  of  control  software  for  MB86064. 


Figure  9.7:  The  function  block  diagram  of  MB86064. 

detector  and  ADC  is  necessary  for  signal  conditioning.  To  facilitate  digital  interface,  the  A/D  converter  features  an 
on-chip,  selectable  8:16  output  demultiplexing  that  slows  the  LSGsps  data  rate  to  750  mega-words/second  ported 
to  two  parallel,  differentia!  8-bit,  low- voltage  (PECL)  outputs.  A  programmable  high-frequency  clock  synthesizer 
is  used  to  generate  the  sampling  clock  for  the  A/D  converter  The  output  frequency  of  the  clock  synthesizer  is 
controlled  by  FPGA,  so  variable  ADC  sampling  rate  is  achieved.  The  number  of  bits  of  the  A/D  converter  is 
dynamically  adapted  to  the  conditions  of  the  environment.  It  can  downscale  the  digital  baseband  complexity  and 
total  power  assumption.  A  double-data-rate  (DDR)  interface  with  FIFO  is  implemented  in  FPGA  to  connect  the 
A/D  converter  to  FPGA.  For  the  design  of  the  high-speed  interface  between  ADC  and  FPGA,  signal  integrity  has 
become  a  critical  issue.  Many  signal  integrity  problems  are  electromagnetic  phenomena  in  nature  and  hence  related 
to  the  EMI/EMC.  There  are  two  concerns  for  signal  integrity  -  the  timing  and  the  quality  of  the  signal.  Signal 
timing  mainly  depends  on  the  delay  caused  by  the  physical  length  that  the  signal  must  propagate.  Signal  waveform 
distortions  can  be  caused  by  reflection,  cross  talk  and  power/ground  noise.  An  interface  board  has  been  carefully 
designed  to  solve  the  signal  integrity  issue. 
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Figure  9.8:  MAX  1 08  function  diagram. 


9.5  SPI  Control 

In  this  system.  There  are  many  devices  need  to  configure  and  control,  such  as  LO,  AGC,  DAC.  Each  one  has  its  own 
way  of  initialization  and  configuration.  In  order  to  reduce  the  system’s  complexity,  SPI  Bus  is  chosen  to  configure 
all  these  devices.  Also  it  could  be  implemented  in  FPGA, which  is  easy  and  efficient. 

SPI  standards  for  Serial  Peripheral  Interface,  It  is  a  simple  interface  that  allows  one  chip  to  communicate  with  one 
or  more  other  chips.  It  can  easily  achieve  a  few  Mbps  ( mega-bits-per-seconds).  Basically:  (I)  It  is  synchronous.  (2) 
It  is  serial.  (3)  It  is  full-duplex.  (4)  (l  is  not  plug-and-play.  (5)  There  is  one  (and  only  one)  master,  and  one  (or  more) 
slaves.  In  more  details:  (I)  A  clock  is  generated  by  the  master,  and  one  bit  of  data  is  transferred  each  time  the  clock 
toggles.  Because  SPI  is  synchronous  and  full-duplex,  every  time  the  clock  toggles,  two  bits  are  actually  transmitted 
(one  in  each  direction).  (2)  Data  is  serialized  before  being  transmitted,  so  that  it  fits  on  a  single  wire.  There  are  two 
wires  for  data,  one  for  each  direction.  (3)  The  master  and  slave  know  beforehand  the  details  of  the  communication 
(bit  order,  length  of  data  words  exchanged,  etc.)  (4)  The  master  is  always  the  one  who  initiates  communication. 
Only  one  slave  is  active  at  a  time. 

A  SPI  master  can  communicate  with  single  stave  as  well  as  multiples  slaves  by  connecting  most  signals  in  parallel 
and  adding  SSEL  lines.  SPI  requires  4  wires  to  be  used  in  between  the  two  chips. As  you  can  see,  the  wires  are  called 
SCK,  MOSI,  MI  SO  and  SSEL,  and  one  of  the  chip  is  called  the  SPI  master,  while  the  other  the  SPI  slave.  Fig.  9. 1 0 
shows  A  SPI  master  communicate  with  one  slave  and  three  slaves. 

(1)  MOSI  -  Master  device  data  out.  Slave  device  data  in  (2)  MISO  -  Master  device  data  in.  Slave  device  data  out  (3) 
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Figure  9.9:  MAX  108  development  board. 


SCLK  -  Clock,  Generated  by  Master  Device  (4)  SSEl  -Slave  device  enable,  controlled  by  Master  device 

When  there  are  three  slaves.  The  master  must  have  three  SSEL  lines,  and  activate  only  one  SSEL  line  at  a  time. 
Since  all  three  slaves  are  connected  to  the  MISQ  line,  slaves  that  are  not  activated  must  not  drive  the  MISO  line. 

Let's  assume  that  the  master  and  slave  expect  Sybils  data  transfers,  with  MSB  transmitted  first.  Fig.  9.1 1  shows  how 
would  look  a  single  8-bits  data  transfer  and  its  the  internal  interface. 

( 1 )  The  master  pulls  SSEL  down  to  indicate  to  the  slave  that  communication  is  starting  (SSEL  is  active  low).  (2)  The 
master  toggles  the  clock  eight  limes  and  sends  eight  data  bits  on  its  MOSI  line.  At  the  same  time  it  receives  eight 
data  bits  from  the  slave  on  the  MISO  line.  (3)  The  master  pulls  SSEL  up  to  indicate  that  the  transfer  is  over. 

If  the  master  had  more  than  one  8-bits  data  to  send/reccive,  it  could  keep  sending  or  receiving  and  rc-asscrt  SSEL 
only  when  it  is  done. 

In  our  system,  SPI  bus  is  used  to  implement  system  control. 

For  LO,  The  evaluation  board  is  supposed  to  come  with  a  cable  for  connecting  to  the  printer  port  of  a  PC.  The  silk 
screen  and  cable  diagram  for  the  evaluation  board  are  shown  below  as  Fig.  9.12.  We  intend  to  replace  the  PC  by 
FPGA,  Using  SPI  protocol  to  write  ADF4106  registers  through  the  DB-9  interface.  The  ADF4106  has  a  simple 
SPI-compatible  serial  interface  for  writing  to  the  device.  CLK,  DATA,  and  LE  control  the  data  transfer. 

For  DAC,  Its  configuration  and  control  are  expected  to  through  Fujitsu’s  PC  USB  Programming  cable  and  accom¬ 
panying  PC  software,  as  well  as  download  test  vectors  to  the  waveform  memory  module.  The  programming  cable 
interfaces  between  a  host  PC’s  USB  port  and  the  4-wire  serial  interface  implemented  on  the  device,  The  connection 
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Figure  9.10:  The  block  diagram  of  SPI  Bus. 
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Figure  9. 1 1 :  The  block  diagram  of  SPI  Transfer, 

is  shown  in  Fig.  9.13. 

Similarly  we  also  intend  to  configure  the  device  by  SPI  protocol.  To  facilitate  this.  An  appropriate  connector  should 
be  fitted  to  the  target  board,  providing  connections  to  the  signals  as:  ( I )  SERIAL  IN.  (2)  SERIAL  OUT,  (3)  SERIAL 
CLK.  (4)  SERI  AL  EN.  Fig.  9.14  shows  this  connection. 


Other  devices  could  be  configured  with  SPI  protocol  similarly. 
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Figure  9, 12:  The  diagram  of  LO  interface. 


Figure  9.13:  The  diagram  of  DAC  Kit, 
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Figure  9.14:  The  diagram  of  DAC  interface. 
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Chapter  10 


Digital  Back-Ends 


Major  function  played  by  the  digital  back-ends  include  system  control/managemcm,  parameter  setting,  AGC,  syn¬ 
chronization,  interfacing  with  PC,  as  well  as  modulation  and  demodulation. 


10.1  System  Control  Flow 

10.2  Filter  Design 


In  order  to  realize  time  reversal  in  the  test-bed,  we  need  to  build  the  time  reversal  waveform  generator  in  the  transmit¬ 
ter  side.  According  to  the  design  of  the  test-bed,  the  FIR  filter  built  in  FPGA  will  be  used  to  generate  the  waveform 
of  time  reversal  signaling.  Meanwhile,  FIR  filter  will  be  used  in  the  receiver  side  to  accomplish  the  functions  of 
matched  filter,  correlation  and  so  on  in  the  digital  domain. 


10.2.1  The  Typical  Methods  for  Filter  Design 

The  typical  methods  for  the  FIR  filter  design  in  FPGA  are  intellectual  property  (IP)  core  based  method  and  user- 
defined  coding  based  method.  IP  core  based  method  means  we  use  IP  core  in  FPGA  to  realize  the  function  of  digital 
filter.  There  arc  two  kinds  of  IP  cores  in  FPGA  can  be  employed  for  this  purpose.  One  kind  is  filter  based  IP  core 
and  the  other  kind  is  memory  based  IP  core.  The  user-defined  coding  based  method  means  wc  directly  write  Vcrilog 
code  to  build  the  FIR  filter  according  to  the  design  specification. 

Multiply  accumulate  (MAC)  FIR  filter  IP  core  and  distributed  arithmetic  (DA)  FIR  filter  IP  core  are  the  frequently 
used  filler  based  IP  cores.  The  MAC  FIR  IP  core  implements  a  highly  configurable,  high-performance,  and  area 
efficient  FIR  filter.  All  internal  data-paths  provide  full-precision  arithmetic  to  avoid  the  possibility  of  overflow.  The 
full-prec  ision  sum-of-products  is  presented  on  the  data  output  port.  A  set  of  three  handshake  control  signals  provides 
an  easy-to-use  user  interface.  The  MAC  FIR  filter  IP  core  uses  one  or  more  time-shared  multiply  accumulate 
functional  units  to  service  the  N  sum-of-product  calculations  in  the  filler  [1].  The  DA  FIR  filler  IP  core  is  also  a 
highly  paramcterizable,  area-efficient  high-performance  FIR  filter.  The  mam  difference  between  the  MAC  FIR  filter 
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and  ihe  DA  FIR  filter  is  the  latter  employs  no  explicit  multipliers  in  the  design,  only  look-up  tables  (LUTs),  shift 
registers,  and  a  scaling  accumulator. 

In  its  most  obvious  and  direct  form,  DA-based  computations  are  bit-serial  in  naturcserial  distributed  arithmetic 
(SDA)  FIR.  The  advantage  of  a  distributed  arithmetic  approach  is  its  efficiency  of  mechanization.  The  basic  opera¬ 
tions  required  are  a  sequence  of  table  look-ups,  additions,  subtractions  and  shifts  of  the  input  data  sequence  [  I  ].  All 
of  these  functions  efficiently  map  to  FPGAs.  Input  samples  arc  presented  to  the  input  parallel-to-serial  shift  register 
(PSC)  at  the  input  signal  sample  rate.  As  the  new  sample  is  serialized,  the  bit-wide  output  is  presented  to  a  bit-serial 
shift  register  or  time-skew  buffer  (TSB).  The  TSB  stores  the  input  sample  history  in  a  bit-serial  format  and  is  used  in 
forming  the  required  inner-product  computation.  The  nodes  in  the  cascade  connection  of  TSBs  arc  used  as  address 
inputs  to  a  look-up  table.  This  LUT  stores  all  possible  partial  products  over  the  filter  coefficient  space. 

In  a  conventional  MAC  FIR  filter  realization,  the  sample  throughput  is  coupled  to  the  filter  length.  As  the  filter  length 
is  increased,  the  system  sample  rate  is  proportionately  decreased.  With  a  DA  architecture,  the  system  sample  rate 
is  related  to  the  bit  precision  of  the  input  data  samples.  The  trade  off  introduced  here  is  one  of  silicon  area  (FPGA 
logic  resources)  for  time.  As  the  filter  length  is  increased  in  a  DA  FIR  filter,  more  logic  resources  are  consumed, 
but  throughput  is  maintained  [I],  Meanwhile,  one  approach  to  increase  the  speed  of  FIR  filter  using  DA  architecture 
is  to  partition  the  input  words  into  M  subwords  and  process  these  subwords  in  parallel  [I],  This  method  requires 
M-times  as  many  memory  look-up  tables  and  so  comes  at  a  cost  of  increased  storage  requirements.  Maximum  speed 
is  achieved  by  factoring  the  input  variables  into  single-bit  subwords.  The  resulting  structure  is  a  fully  parallel  DA 
(PDA)  FIR  filter.  With  this  factoring  a  new  output  sample  is  computed  on  each  clock  cycle.  Another  advantage  of 
the  DA  FIR  filter  IP  core  is  the  on-line  coefficient  reload  which  provides  additional  flexibility  to  the  system  design. 
While  the  new  coefficient  values  are  being  loaded,  and  some  internal  data  structures  are  subsequently  initialized, 
the  filter  ceases  to  process  input  samples.  The  coefficient  reload  time  is  a  function  of  the  filter  length  and  type  [I], 
When  we  use  the  DA  FIR  filter  IP  core  to  build  the  single  rate  FIR  filter,  the  reload  delay  is  210  clock  cycles  if  the 
number  of  taps  is  10,  the  coefficient  width  is  10  and  the  input  width  is  I.  Fig.  10.3  shows  the  simulation  result  of 
the  coefficient  timing  using  DA  FIR  filter  IP  core  in  FPGA. 

If  wc  want  to  use  memory  based  IP  core  to  realize  the  FIR  filter,  FIFO  Generator  IP  core  and  Block  Memory 
Generator  IP  core  can  be  employed.  The  basic  thought  of  this  method  is  that  the  data  of  the  candidate  waveform  are 
put  into  the  memory  firstly  and  then  the  data  are  read  from  the  memory  according  to  some  control  signals  later.  IP 
core  in  FPGA  has  powerful  functions.  Even  though  IP  core  can  work  efficiently,  it  will  occupy  more  resources  than 
we  expect.  So  if  our  design  is  not  complex,  we  can  directly  write  Verilog  code  to  build  FIR  filter.  For  example,  wc 
can  build  the  shift  register  and  adder  by  coding  to  generate  the  waveform  of  time  reversal  signaling. 


10.2.2  Quantization  Consideration 

Quantization  issue  is  very  important  for  the  FIR  filter  design  in  FPGA,  because  any  number  in  FPGA  should  be  finite 
word  length.  Quantization  will  introduce  error  in  the  filter  implementation.  There  are  three  sources  of  quantization 
error  -  coefficient  quantization,  input  quantization  and  the  internal  rounding  that  is  the  quantization  of  the  results 
of  arithmetic  operations  within  the  filter  [2].  According  to  [2],  if  each  coefficient  is  quantized  to  m  +  1  bits  and 
the  input  is  also  quantized  to  n  +  1  bits  where  n  <  m,  the  output  quantization  noise  power  is  nearly  optimal  for 
internal  rounding  precision  set  to  less  than  m  4-  I  and  greater  than  n  +  I  bits  if  there  is  no  restriction  on  the  output 
quantization.  However,  in  our  test-bed,  the  output  quantization  is  fixed  and  set  to  1 4  bits,  because  the  input  resolution 
of  DAC  is  14  bits.  Meanwhile  if  we  use  shift  register  to  build  time  reversal  waveform  generator,  we  will  only  have 
coefficient  quantization  and  internal  rounding  problems.  If  the  time  reversal  waveform  has  no  overlapping  in  the 
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transmitter  side,  the  coefficient  quantization  should  be  14  bits  which  is  the  same  as  the  output  quantization.  Based 
on  Dr  Guo’s  research  results,  nearly  half  of  the  ideal  energy  of  the  channel  impulse  response  can  be  attained  in  the 
receiver  side  using  I  bit  coefficient  if  inter-symbol  interference  (ISI)  is  not  present  Because  mono-bit  resolution  will 
cause  the  big  side  lobe  in  the  received  signal.  So,  if  waveform  overlapping  is  required  in  the  transmitter  side  and  ISI 
is  present  in  the  receiver  side,  the  resolution  of  the  coefficient  and  the  internal  rounding  should  be  increased  in  order 
to  suppress  the  side  lobe  and  increase  the  energy  of  the  main  lobe  in  the  received  signal.  Because  the  channel  state 
varies  in  different  scenarios,  we  should  make  the  tradeoffs  between  the  system  performance  and  system  complexity 
to  determine  the  resolution  of  the  coefficient  and  the  internal  rounding  when  the  output  resolution  is  fixed  and  given* 
This  research  issue  will  be  further  studied. 

In  sum,  no  matter  what  method  we  use,  we  should  follow  some  steps  to  design  the  FIR  filter  in  FPGA,  (1)  Use 
Mat  lab  to  design  a  floating-point  FIR  filler,  (2)  Convert  the  floating-point  FIR  filter  into  fixed-point  FIR  filter  using 
Matlab.  (3)  Build  a  hardware  description  in  FPGA  and  verify  that  the  fixed-point  design  functions  correctly.  (4) 
Evaluate  the  fixed-point  FIR  filter  in  terms  of  its  performance  and  hardware  efficiency  [3].  (5)  Program  FPGA  and 
verify  the  FIR  filler’s  performance  experimentally  on  the  workbench. 
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Figure  10. 1 :  System  control  flow  for  transmitter. 
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Figure  !  0.2:  System  control  flow  for  receiver. 
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Figure  10.3:  The  simulation  result  of  coefficient  reload  timing  using  DA  FIR  filter  IP  core  in  FPGA. 
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Chapter  11 


Graphic  User  Interface  (GUI)  between  FPGA 
and  PC 

In  ihis  section  wc’II  introduce  the  PC  interface  in  our  test-bed  system.  The  PC  interface  is  responsible  for  in 
put/output  communications  between  host  PC  and  FPGA. 


11.1  Architecture 

Fig.  M .  I  shows  the  architecture  of  the  PC/FPGA  connection  with  PC  interface.  Here  we  use  laptop  as  the  host  PC, 
since  it’s  powerful  enough  and  easy  to  carry.  The  data  acquisition  (DAQ)  board  acts  as  the  PC  interface.  For  digital 
data  output,  DAQ  gets  data  from  PC  then  send  it  to  FPGA.  For  digital  data  input,  DAQ  samples  data  from  FPGA 
then  sends  it  to  PC.  We  use  USB  cable  for  the  interface  between  PC  and  DAQ  while  flying  lead  wire  for  DAQ/FPGA 
connection  (Fig.  11.2). 


1 1.2  Communications  between  FPGA  and  host  PC 


FPGA  is  extensively  used  in  our  test-bed  system.  However,  if  we  want  to  change  the  parameters  in  FPGA  for  testing, 
there  will  be  some  difficulties.  We  can  let  the  FPGA  change  the  parameters  by  itself,  but  this  will  take  up  the  valuable 
FPGA  resource.  Also,  FPGA  does  not  have  graphical  user  interface  (GUI).  We  must  change  the  codes  and  download 
them  to  FPGA  each  lime,  which  is  time  consuming  and  not  intuitive.  With  a  PC  interface  between  host  PC  and 
FPGA,  we  can  bypass  these  difficulties.  We  can  set  up  the  parameters  in  GUI  then  transmit  them  to  FPGA  via  PC 


Figure  1 1 . 1 :  The  architecture  of  the  PC/FPGA  connection  with  PC  interface 
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Figure  1 1 .2:  The  flying  lead  wire  interface. 


Figure  1 1 .3:  Byte  Paradigm,  GP22050  DAQ  board 


interface  without  making  any  changes  in  FPGA.  This  is  much  more  efficient  and  intuitive. 

We  can  also  use  PC  interface  to  read  data  from  FPGA.  Though  this  function  can  be  done  by  logic  analyzer,  it  can 
only  display  static  signals,  not  real  time  signals.  Also,  it  does  not  support  real  time  signal  processing.  With  PC 
interface  bridged  between  FPGA  and  PC,  real  time  data  can  be  sampled  and  stored  to  PC.  After  capturing,  the  data 
can  be  processed  and  viewed  in  the  monitor  in  real  time. 


U.3  Implementation 


We  choose  Byte  Paradigm’s  GP22050  as  our  DAQ  (Fig.  1 1.3).  It  is  a  powerful  DAQ  supporting  up  to  1 6  bidirectional 
data  and  6  bidirectional  control  signals.  The  maximum  processing  rate  is  50MHz.  Therefore,  the  maximum  burst 
throughput  can  be  up  to  lOOMBytc/s.  GP22050  is  also  provided  with  C  language  interface.  The  C  functions  can  be 
called  to  control  DAQ  by  Visual  C  (VC),  Labview  and  other  programs  [1,2],  Thus,  a  user-friendly  GUI  can  be  made 
to  control  the  PC  interface  between  DAQ  and  FPGA. 
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Figure  1 1 .4:  Data  stream  structure  for  PC/DAQ  connection. 
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Figure  1 1 .5:  The  relationship  of  the  clock,  trigger  and  data  signals. 


11.4  Application:  BER  Measurement 


Now  we  are  developing  the  bit  error  rate  (BER)  measurement  application.  To  get  the  BER,  we  need  to  compare  the 
received  data  with  transmitted  data.  To  simplify  the  process,  the  transmitted  data  periodic  and  known  at  the  receiver. 
Also,  at  the  beginning  of  each  period,  there  is  a  trigger  signal  to  inform  the  program  when  to  start  comparing.  In 
our  PC  interface  based  system,  the  received  data,  trigger  signal  and  clock  signal  are  the  output  of  FPGA.  DAQ  uses 
three  data  streams  to  get  those  signals.  The  data  stream  structure  is  illustrated  in  Fig.  1 1.4.  DAQ  will  not  start 
sampling  until  it  receives  the  trigger  signal.  After  5  clock  cycles  of  the  trigger  signal  [3],  DAQ  samples  data  from 
FPGA  at  each  clock  rising  edge.  The  relationship  of  the  clock,  trigger  and  data  signals  is  illustrated  in  Fig.  1 1 .5. 
DAQ  samples  at  rising  edge  of  a  clock  and  it  only  starts  sampling  after  trigger  signal  is  detected  [3].  The  maximum 
sampling  rate  is  the  same  as  DAQ’s  maximum  sampling  rate,  which  is  50MIIz  (4],  Since  DAQ  can  sample  16  data 
lines  in  parallel,  the  maximum  throughput  can  be  lOOMbyte/s  [4]. 

Fig.  1 1,6  is  the  BER  measurement  software  GUI.  Before  BER  measurement,  we  need  to  initialize  the  measurement 
environment  in  the  “Settings”  dialog  first  (Fig.  1 1.7).  If  “Stop  Manually”  is  chosen,  the  software  will  keep  running 
until  the  “Stop”  button  in  Fig.  1 1.6  is  clicked.  Otherwise,  it  stops  when  it  reaches  “Number  of  bits”  or  “Errors”  set 
by  user.  BER  values  will  be  updated  every  100ms  in  the  “BER”  edit  box.  A  new  BER  value  will  be  plotted  onto  the 
“Real  time  display”  box  every  1  second.  BER  VS  Distance  (or  SNR)  curves  will  be  displayed  on  “Statistic  display” 
box  for  each  environment.  They  can  be  saved  and  loaded  any  time. 

This  software  is  still  under  construction.  It  is  expected  to  be  finished  in  one  week. 
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Figure  I !  ,6:  The  BER  software  GUI. 
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Figure  11.7:  The  Settings  dialog.  This  dialog  is  used  to  initialize  the  measurement  environment. 
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11.5  Application:  FPGA  Test 

Wc  arc  going  to  develop  an  FPGA  test  software.  In  order  to  test  FPGA  codes,  we  need  to  change  the  parameters 
frequently,  such  as  FIR  filter’s  coefficients,  system  data  rate,  transmitted  data,  etc..  Since  FPGA  docs  not  have  GUI, 
we  can  either  change  the  code  each  time  or  let  FPGA  change  the  parameters  by  itself.  Both  of  these  methods  are  lime 
consuming  and  not  intuitive.  With  DAQ,  however,  FPGA  will  be  able  to  read  parameters  from  PC’s  GUI  without 
changing  the  codes.  This  is  a  much  more  intuitive  way  and  saves  time. 

GP22050’s  "arbitrary  digital  waveform  generator”  (ADWG)  function  [4]  can  be  applied  to  our  FPGA  test  applica¬ 
tion.  The  ADWG  function  is  able  to  generate  arbitrary  16  bit-depth  output  at  a  rate  of  50MHz.  As  a  result,  GP22050 
supports  setting  16  coefficients  simultaneously  in  FPGA  at  a  high  speed.  With  the  ADWG  Library  [2]  provided  by 
the  device,  we  can  customize  the  GUI  design  and  make  the  test  much  more  intuitive. 


11.6  Limitations 

Although  PC  interface  provides  a  way  for  communications  between  FPGA  and  PC,  the  processing  speed  is  limited. 
Since  PC  utilizes  a  lot  of  resources  to  support  the  operation  system  and  other  necessary  softwares,  the  signal  pro¬ 
cessing  speed  is  far  lower  than  the  speed  of  the  DAQ.  According  to  result  of  a  test  with  a  Pentium  4  3  GHz  PC,  the 
maximum  continuous  throughput  is  only  1  IMByte/s,  which  is  far  less  than  the  burst  throughput  of  lOOMBylc/s. 


142 


CHAPTER  1 ! .  GRAPHIC  USER  INTERFACE  (GUI)  BETWEEN  FPGA  AND  PC 


Bibliography 


[1]  C/C++  library  for  Analyser  application,  Byle  Paradigm,  Belgium,  2006. 

[2]  C/C++  library  for  ADWG  application.  Byte  Paradigm,  Belgium,  2006. 

[3]  HP!  Control  Panel  Users  Guide  for  G P-22050,  Byte  Paradigm,  Belgium,  2006. 

[4]  G P-2 2050  complete  data  sheet,  Byte  Paradigm,  Belgium,  2006. 


143 


